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ABSTRACT 


A  nondivergent,  barotropic  analytical  model  to  predict  steady  tropical  cyclone  (TC) 
propagation  relative  to  the  large-scale  environment  is  developed  in  terms  of  a  "£elf- 
advectiori"  process  in  which  the  TC  is  advected  by  an  azimuthal  wavenumber  one  gyre 
flow  that  results  from  TC-environment  interaction.  The  model  is  comprehensive  in  that 
it  includes  the  first-order  effects  of  all  of  the  dynamic^7 influences  that  are  presently  un¬ 
derstood  to  be  important  to  barotropic  propagation:  gradients  of  planetary  and  envi¬ 
ronmental  vorticity,  changes  in  TC  wind  structure,  and  environmental  windshear.  An 
unforced  version  of  the  model  is  used  to  show  that  angular  windshear  in  the  symmetric 
TC  circulation  acts  to  damp  perturbations  from  axisymmetry  by  tilting  the  perturbations 
downshear.  The  resultant  transfer  of  kinetic  energy  from  perturbation  to  symmetric 
circulation  thus  tends  to  restore  axisymmetry.  Thus,  steady  propagation  of  TC-like 
barotropic  vortices  is  a  manifestation  of  a  stable  response  to  asymmetric  forcing.  To 
predict  both  the  asymmetric  gyre  flow  and  the  propagation  it  induces,  the  forced 
Barotropic  Self-Advection  Model  (BSAM)  is  closed  by  seeking  a  particular  pattern  in  the 
vorticity  tendencies  of  the  TC-environmental  interaction  flow.  For  realistic  combina¬ 
tions  of  environmental  vorticity  gradients  and  linear  windshear,  the  BSAM  predicts 
propagation  speeds  and  directions  that  are  consistent  with  TC  propagation  character¬ 
istics  observed  in  composite  data.  The  capability  of  the  BSAM  to  account  for  variable 
TC  structure  is  used  to  show  that  errors  in  determining  TC  outer  wind  strength  of  ±1 
m/s  can  result  in  an  85  km  forecast  error  at  48  h.  Finally,  and  most  importantly,  the 
capability  of  the  BSAM  to  initialize  a  barotropic  numerical  model  so  that  quasi-steadv 
TC  propagation  occurs  almost  immediately  is  demonstrated  for  several  simple  dynamical 
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I.  INTRODUCTION 


It  is  well  known  that  tropical  cyclone  (TC)  motion  persistently  differs  from  the  mean 
advective  effect  (i.e.,  the  "steering")  of  the  large-scale  environment  (e.g.,  George  and 
Gray  1976;  Chan  and  Gray  1982).  The  precise  speed  and  direction  of  this  velocity  dif¬ 
ference  (hereafter  referred  to  as  "propagation")  depends  to  some  extent  on  the  particular 
steering  flow  definition  employed,  but  is  generally  westward  and  poleward  in  direction 
and  1  to  2.5  m/s  in  magnitude  (Carr  and  Elsberry  1989).  Although  this  propagation 
velocity  is  typically  smaller  than  the  advection  by  environmental  steering,  it  can  induce 
significant  TC  track  errors  (2  m  s  *  170  km  d)  if  not  properly  accounted  for  in  dynam¬ 
ical  models.  Interestingly,  it  is  well  known  that  numerical  model  forecasts  of  TC  tracks 
generally  have  less  skill  than  CLIPER  (CLlmatologv  and  PERsistance)  over  the  first 
2-1-36  h,  and  that  the  error  is  primarily  an  underestimation  of  TC  speed  (e.g.,  Neumann 
and  Pelissicr  1981;  their  Figs.  4,8  respectively).  Since  TC's  track  generally  westward  and 
poleward  throughout  much  of  their  lifetimes,  such  statistics  provide  circumstantial  evi¬ 
dence  that  dynamical  forecast  models  do  not  as  yet  properly  account  for  TC  propagation 
until  about  1  day  into  the  numerical  integration.  It  will  be  shown  that  the  initial 
slowness  of  numerical  models  may  result  from  initializing  the  model  with  a  symmetric 
TC  structure  that  lacks  any  propagation-inducing  asymmetries  caused  by 
TC-environment  interaction.  Although  past  modeling  experience  indicates  that  both 
barotropic  and  baroclinic  processes  can  cause  TC  propagation  asymmetries,  the  general 
goal  of  this  research  is  to  develop  an  improved  theoretical  model  for  understanding  and 
predicting  the  barotropic  contributions  to  TC  propagation  asymmetries. 

A.  MODELING  REVIEW 

1.  Status  of  basic  understanding 
a.  /?- induced  propagation 

The  propagation  of  an  initially  symmetric,  nondivergent  barotropic  (NDBT) 
vortex  in  a  quiescent  environment  due  to  the  influence  of  planetary  rotation,  as  ap¬ 
proximated  by  a  />-plane,  provides  the  most  basic  dynamical  scenario  for  beginning  a 
review  of  TC  motion  theory.  The  northwestward  propagation  at  2-3  m  s  noted  by  Chan 
and  William  (1987;  hereafter  CW)  in  a  recent  numerical  study  of  this  problem  is  gener¬ 
ally  consistent  with  the  propagation  observed  in  previous  numerical  studies  (Anthcs  and 
Iloke  1975;  Kitadc  1981;  Dc. Maria  1985,  hereafter  DM).  The  authors  of  the  previous 
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numerical  studies  attributed  their  results  to  the  traditional  linear  theories  of  barotropic 
vortex  propagation  advocated  in  various  forms  ( e.g .,  Rossby  1939, 1948;  Kasahara  1957; 
Kasahara  and  Platzmann  1963;  Holland  1983  ).  An  interesting  aspect  of  these  theoretical 
predictions  is  that  the  westward  and  poleward  components  of  the  propagation  are  asso¬ 
ciated  with  different  linear  mechanisms,  e.g.,  Rossby  wave  propagation  by  the  TC  to  the 
west  as  opposed  to  poleward  acceleration  due  to  net  Coriolis  force  on  the  TC.  In  con¬ 
trast,  CW  concluded  that  only  the  nonlinear  "self-advection"  (i.e.,  the  TC  advecting  its 
own  vorticity)  process  associated  with  the  development  of  a  horizontal  asymmetry  in  the 
TC  windfield  was  responsible  for  both  the  westward  and  poleward  components  of  the 
propagation. 

The  basis  for  this  important  distinction  can  be  demonstrated  by  partitioning 
the  vortex  flow  into  a  steady  symmetric  (S)  component  and  an  asymmetric  (A)  compo¬ 
nent  that  is  generally  small  by  comparison.  In  a  coordinate  system  centered  at  any  in¬ 
stant  on  the  symmetric  component,  the  propagation  process  in  a  quiescent  environment 
is  approximately  described  by 

4f-  =  -  v<.v:5  -  vs.vo  -  vs.v/,  (i. i) 

C  l 

(a)  (b)  (c)  (d) 

where  term  (a)  represents  the  vorticity  tendency  associated  with  motion  and  or  dis¬ 
tortion  of  the  TC  due  to  the  combined  effect  of  the  two  dominant  self-advection  terms 
(b  and  c)  and  the  advection  of  planetary  vorticity  (d).  Linear  theories  regard  the  west¬ 
ward  component  of  propagation  as  a  Rossby  wave-like  propagation  associated  with  the 
east-west  dipole  of  vorticity  tendency  generated  by  term  (d),  and  either  eliminate  self- 
advection  by  requiring  the  TC  to  be  symmetric  at  all  times  (e.g.,  Kasahara  and  Platzman 
1963),  or  assume  that  self-advection  is  implicit  in  computations  of  observed  steering 
flows  (e.g.,  Holland  1983,  1984).  When  the  nonlinear  self-advection  term  was  omitted 
from  their  model.  CW  noted  that  advection  of  planetary  vorticity  by  the  TC  is  a  strong 
dispersive  distorting  influence,  but  by  itself  generates  negligible  motion  of  the  TC  center 
(5:0.3  m  s).  Thus,  they  concluded  that  the  linear  process  generates  the  asymmetric  flow 
essential  to  the  existence  of  the  self-advection  process,  but  the  self-advcctive  interaction 
of  the  symmetric  and  asymmetric  TC  components  actually  causes  the  propagation. 

I'iorino  (1087)  and  Tiorino  and  Llsberry  (1989;  hereafter  FT)  clarified  the 
individual  roles  of  the  two  primary  self-advection  processes  (terms  b  and  c)  via  a  diag- 


nostic  technique  that  extracts  \s  and  \'A  from  a  numerical  model  similar  to  that  used 
by  CW.  Their  analysis  showed  that  the  structure  of  \A  is  predominantly  an  azimuthal 
wavenumber  one  asymmetry  that  consists  of  two  counterrotating  gyres  (hereafter  re¬ 
ferred  to  as  "wavenumber  one  gyres")  that  produce  a  near-uniform  flow  across  the  cen¬ 
tral  region  of  the  TC.  The  average  velocity  of  the  uniform  flow  is  very  nearly  equal  to 
the  TC  propagation  vector,  which  indicates  that  NDBT  /J-induced  propagation  is  es¬ 
sentially  an  advcction  of  symmetric  TC  vorticity  by  the  asymmetric  flow  component  of 
the  TC,  i.e.,  term  (b)  of  (1.1).  This  is  somewhat  surprising  since  simple  scale  arguments 
suggest  that  the  advection  of  asymmetric  vorticity  by  the  symmetric  tangential  wind 
(term  c)  should  be  equally  important  to  the  propagation  process.  However,  FE  demon¬ 
strated  that  streamfunction  tendencies  associated  with  term  (c)  nearly  cancel  those  as¬ 
sociated  with  term  (d)  after  a  quasi-steady  state  is  achieved.  Thus,  term  (c)  acts 
primarily  to  limit  the  growth  of  the  wavenumber  one  gyres  and  the  resultant  propagation 
to  nearly  steady  values  after  an  initial  adjustment  period  of  approximately  24-48  h. 

A  theoretical  model  that  accurately  predicts  ^-induced  gyre  structure  and 
the  associated  propagation  based  on  this  recently  identified  self-advection  process  has 
yet  to  be  developed,  and  this  is  the  principal  objective  of  this  research.  The  successful 
development  of  such  a  model  that  confirms  the  numerical  work  of  CW  and  FE  will  be 
an  important  step  in  establishing  the  nonlinear  self-advection  principle  as  the  mechanism 
responsible  for  barotropic  TC  propagation  as  opposed  to  the  traditional  linear  Rossby 
wave  arguments. 

Since  advection  of  planetary  vorticity  by  a  vortical  flow  represents  essen¬ 
tially  steady  asymmetric  forcing,  FE's  results  suggest  that  advection  by  the  symmetric 
tangential  wind  tends  to  stabilize  a  barotropic  approximation  of  a  TC  to  such  forcing. 
Nonlinear  advection  is  necessary  to  permit  dispersion-resistant  movement  of  model 
vortices  in  general  (e.g.,  McWilliams  and  Flierl  1979;  Mied  and  Lindemann  1979),  or 
nondispersive  movement  in  the  case  of  very  specialized  entities  known  as  either  modons 
or  solitary  eddy  solutions  (see  Flierl  a  al.  19S0).  However,  the  precise  mechanism  by 
which  nonlinear  advection  stabilizes  a  vortex  of  arbitrary  radial  structure  has  not  been 
identified.  Because  steady  barotropic  vortex  propagation  is  essentially  a  manifestation 
of  this  stability  to  asymmetric  forcing,  a  thorough  analysis  of  this  process  is  important 
to  improved  understanding  of  barotropic  TC  propagation,  and  thus  is  another  objective 
of  this  research. 

It  is  important  to  briefly  address  the  question  whether  barotropic  vortex 
propagation  stability  is  a  sufficiently  important  process  compared  to  potential  baroclinic 


influences  to  warrant  the  present  barotropic  modeling  effort.  Although  Holland  (1983) 
has  suggested  that  inertial  stability  accounts  for  near  axisymmetrv  in  the  inner  300  km 
of  a  TC,  inertial  stabilization  cannot  occur  in  NDBT  or  quasi-geostrophic  models.  The 
qualitative  similarity  of  TC  propagation  tracks  in  simple  baroclinic  models  (e.g., 
Madala  and  Piacsek  1975;  Kitade  1980)  to  barotropic  results  provides  at  least 
circumstantial  evidence  that  baroclinic  TC  stability/ propagation  is  more  a  modification 
to,  as  opposed  to  being  fundamentally  different  from,  the  barotropic  phenomenon. 
Thus,  it  is  reasonable  to  assume  that  TC  propagation,' stability  is  fundamentally  a 
vorticity  advection  process  that  should  be  fully  understood  within  a  barotropic  context 
before  additional  baroclinic  complexities  are  added. 
b.  Influence  of  TC  structure 

By  using  an  initial  vortex  specification  that  permits  independent  variations 
of  the  inner  and  outer  tangential  wind  structure,  FE  clearly  showed  TC  propagation  is 
strongly  dependent  on  outer  wind  changes  (strength),  but  is  virtually  independent  of  in¬ 
ner  wind  changes  (intensity)  as  predicted  by  earlier  researchers  (Holland  1983;  DM).  In 
particular,  an  increase  of  only  a  few  m.'s  in  the  initial  tangential  wind  between  300  and 
800  km  results  in  a  faster  and  more  westward  TC  propagation  that  can  significantly  alter 
the  /?-induced  TC  track.  Although  the  initial  value  of  TC  total  relative  angular  mo¬ 
mentum  (RAM)  has  been  cited  as  having  a  potentially  important  influence  on  TC  tracks 
(DeMaria  1987),  Shapiro  and  Ooyama  (1989)  have  recently  shown  that  /hinduced  mo¬ 
tion  does  not  depend  on  initial  TC  RAM  in  any  well-defined  way.  Since  RAM  is  es¬ 
sentially  an  integral  (and  therefore  not  unique)  measure  of  TC  strength,  TC  propagation 
variations  due  to  changes  in  RAM  are  likely  a  manifestation  of  the  dependence  on  TC 
strength  described  by  FE. 

In  view  of  the  above,  it  is  clear  that  properly  initializing  dynamical  forecast 
models  for  TC  strength,  and  ensuring  that  the  model  dynamics  physics  maintain  an  ac¬ 
curate  TC  outer  wind  profile  throughout  the  integration  are  important  areas  for  ongoing 
research.  Without  aircraft,  improved  measurements  of  TC  structure  for  initializing  op¬ 
erational  forecast  models  can  only  be  achieved  by  advancement  in  remote  sensing  ca¬ 
pabilities.  Designing  utilizing  remote  sensors  to  better  detect  TC  strength  and  choosing 
model  parametcrizations  to  best  maintain  a  desired  TC  structure  both  depend  on  im¬ 
proved  understanding  of  the  nature  and  sensitivity  of  TC  propagation  dependence  on 
outer  wind  strength.  Thus,  an  additional  objective  of  this  research  is  to  identify  the  dy¬ 
namical  basis  for  this  dependence  and  provide  at  least  a  preliminary  assessment  of  the 
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wind  measurement  accuracy  required  to  adequately  account  for  the  influence  of  TC 
outer  wind  strength  on  the  barotropic  propagation  process. 

c.  Influence  of  divergence 

Anthes  and  Hoke  (1975)  noted  significant  differences  in  the  ^-induced  TC 
motion  tracks  predicted  by  nondivergent  and  divergent  (5  km  depth)  barotropic  numer¬ 
ical  models,  which  they  attributed  to  divergence-induced  slowing  of  Rossby  wave  phase 
speeds.  Such  a  result  would  suggest  that  nondivergent  forecast  models  such  as 
SAXBAR  (Sanders  ei  al.  1975)  fail  to  account  for  a  significant  propagation  process. 
In  a  model  comparison  similar  to  Anthes  and  Hoke,  Shapiro  and  Ooyama  (1989)  found 
nearly  identical  propagation  tracks,  even  though  substantially  more  divergence  was  in¬ 
cluded  ( 1  km  depth).  Thus,  Shapiro  and  Ooyama  concluded  that  the  earlier  result  was 
due  to  dissimilar  initial  TC  structures  that  W’ere  employed  in  the  divergent  and  nondi¬ 
vergent  models.  This  result  is  important  for  clarifying  the  minimal  role  divergence  plays 
in  barotropic  TC  propagation,  and  also  supports  the  choice  of  a  NDBT  dynamical 
framework  for  the  theoretical  model  employed  in  this  study.  In  addition,  their  result  re¬ 
emphasizes  the  important  influence  of  initial  TC  specification  on  the  accuracy  of  TC 
tracks  predicted  by  a  barotropic  model. 

d.  Environmentally-induced  propagation 

Barotropic  models  also  have  been  used  to  study  the  influence  of  horizontally 
nonuniform  environmental  winds  on  TC  propagation  relative  to  steering.  The  additional 
influences  on  the  propagation  process  may  be  described  by  rewriting  (1.1 )  as 
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in  which  the  environmental  (F)  processes  V£ •  V(/  +  £f)  and  d^Jdt  have  been  omitted 
since  they  relate  to  the  evolution  of  the  steering  flow.  Terms  (a),  (b)  and  (c)  of  (1.2)  are 
the  same  as  in  (1.1).  However,  term  (d)  of  (1.2)  now  includes  advection  of  environmental 
vorticity  by  the  TC,  and  a  new  term  (e)  accounts  for  both  steering  and  shearing  of  the 
symmetric  TC  and  wavenumber  one  gyres  by  the  environmental  wind. 

With  the  addition  of  variable  environmental  winds,  it  becomes  essentially  a 
matter  of  convention  whether  V„  should  be  associated  with  the  symmetric  TC  (Vs)  in 
which  case  the  term  "self-advection"  is  still  meaningful,  or  whether  V*  should  be  re¬ 
garded  as  a  modification  to  V£,  and  thus  a  contribution  to  steering.  The  comparison  of 
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observational  data  with  theory  presented  in  Section  C  will  provide  strong  evidence  to 
suggest  that  V*  will  be  manifested  as  propagation  rather  than  a  contribution  to  con¬ 
ventionally  computed  environmental  steering.  Thus,  V,  will  continue  to  be  referred  to 
as  a  self-advection  flow  associated  with  TC-environment  interaction,  and  will  be  re¬ 
garded  as  distinct  from  the  steering  associated  with  V£. 

In  a  NDBT  numerical  study  with  steady  zonal  winds  that  varied  only  with 
latitude,  DM  noted  that  the  poleward  component  of  propagation  associated  with  a  TC 
in  a  poleward  (equatonvard)  gradient  of  environmental  vorticity  increased  (decreased) 
the  northward  propagation  component  relative  to  p  alone.  DM  attributed  his  results  to 
the  linear  theory  of  Kasahara  and  Platzman  (1963),  which  predicts  that  horizontal  vari¬ 
ability  of  the  environmental  winds  will  induce  TC  propagation  in  the  direction  of  and 
90°  to  the  left  of  the  large-scale  gradient  of  absolute  vorticity.  This  theory  is  essentially 
a  direct  extension  of  Rossby's  westward  propagation  (1939)  and  poleward  acceleration 
(1948)  of  TC's  due  to  p  alone.  The  presence  of  in  term  (d)  of  (1.2)  would  make 
Kasahara  and  Platzman's  hypothesis  intuitively  plausible.  However,  the  recent  demon¬ 
stration  by  C\V  that  linear  advection  of  planetary  vorticity  only  distorts  the  TC  must 
also  apply  to  the  advection  of  environmental  vorticity  by  the  TC.  DM's  results  clearly 
indicate  that  environmental  vorticity  gradients  alter  barotropic  TC  propagation.  This 
suggests  that  Kasahara  and  Platzmann's  theory  of  a  linear  absolute  vorticity-induced 
propagation  due  to  the  combined  effects  of  p  and  can  be  recast  within  the  context  of 
the  nonlinear  self-advection  mechanism  recently  identified.  Thus,  another  objective  of 
this  research  will  be  to  determine  the  extent  to  which  a  self-advection  model  for 
/j-induced  propagation  can  also  accomodate  the  presence  of  a  horizontally  variable 
windfield. 

Any  analogy  between  the  propagation-inducing  effects  of  /?  and  is  com¬ 
plicated  by  the  additional  distorting  influence  of  environmental  wind  shear  in  term  (e) 
of  (1.2).  In  a  barotropic  model  that  includes  P  and  an  environmental  current  with  only 
linear  shear,  Chan  and  Williams  (1989)  found  that  cyclonic  (anticyclonic)  shear  induces 
a  westward  (eastward)  curving  TC  propagation  relative  to  the  tracks  induced  by  /?  alone. 
They  also  noted  that  the  orientation  of  the  wavenumber  one  gyres  rotates  consistently 
with  the  propagation  changes,  which  suggests  that  this  process  is  due  to  shearing  of  gyre 
vorticity  by  the  rotational  part  of  the  environmental  current.  Although  curving  TC 
tracks  are  also  evident  in  the  results  of  DM  who  used  a  sinusoidal  windfield.  the  effect 
is  much  smaller  than  that  seen  by  Chan  and  Williams.  Because  of  the  potentially  im- 
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portant  impact  of  environmental  shear  on  TC  propagation  direction,  modeling  of  this 
linear  shear  process  will  be  included  in  this  research. 

2.  Status  of  propagation/gyre  prediction  models 
a.  Recent  results 

The  preceding  discussion  of  the  self-advection  mechanism  for  barotropic 
TC  propagation  has  been  based  on  diagnostic  analyses  of  numerical  modeling  results. 
Theoretical  models  of  TC  motion  that  explicitly  predict  wave  number  one  gyre  structure 
and  associated  propagation  based  on  the  self-advection  concept  are  also  being  devel¬ 
oped.  The  limited  success  of  such  models  as  discussed  below  has  provided  motivation 
for  this  research. 

Given  the  intense  advective  influence  of  the  TC's  tangential  winds,  the 
quasi-steady  nature  of  the  /?-induced  wavenumber  one  gyre  observed  by  FE  suggests  that 
it  is  a  standing  azimuthal  wave.  Such  a  wave  would  propagate  clockwise  in  a  Rossby 
wave-like  manner  on  the  basic  state  vorticity  gradient  represented  by  the  symmetric 
component  of  the  TC  via  term  (b)  of  (1.1).  On  the  basis  of  such  a  hypothesis  and  an 
assumption  as  in  (1.1)  and  (1.2)  that  <  <  Vs  over  a  large  area,  Willoughby  (1988) 
formulated  a  linear  shallow  model  that  explicitly  predicts  the  structure  of  such  "vorticity 
waves"  as  well  as  the  propagation  they  induce.  He  solved  the  potential  vorticity  and 
divergence  equations  in  a  polar  coordinate  system  moving  with  the  TC  for  a  range  of 
possible  propagation  velocities,  and  assumed  that  the  correct  propagation  velocity  would 
minimize  the  Lagrangian  of  the  system  as  required  by  Hamilton's  Principle  from  vari¬ 
ational  calculus.  By  incorporating  a  rotating  mass  source-sink  pair  to  approximate 
asymmetric  eyewall  convection,  he  obtained  a  vorticity  wave  angular  frequency  that 
agrees  well  with  that  of  small-scale  trochoidal  oscillations  observed  in  TC  tracks.  This 
result  suggests  that  TC  track  oscillations  are  also  explainable  in  terms  of  a  barotropic 
self-advection  process  associated  with  asymmetric  convection,  as  opposed  to  the  tradi¬ 
tional  Magnus  force  theory  (Kuo  1950;  Yell  1950). 

Willoughby's  model  failed  to  predict  a  correct  propagation  velocity  on  a 
/?-plane  (quiescent  environment),  presumably  due  to  a  barotropic  instability  studied  by 
Peng  and  Williams  (1989)  in  a  nondivergent  version  of  the  model.  Peng  and  Williams 
showed  that  existence  of  the  instability  depends  on  the  vorticity  gradient  sign  reversal 
of  the  TC  (necessary  condition)  and  sufficient  model  resolution  to  resolve  the  small  scale 
of  the  instability  (radius  of  maximum  winds).  Thus,  they  suggest  that  this  linear  insta¬ 
bility  is  dynamically  based,  but  is  highly  damped  in  full  numerical  models  due  to  a 
combination  of  nonlinear  interaction  and  coarse  horizontal  resolution!!.  However,  ana- 
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lytical  gyre/propagation  models  such  as  the  one  discussed  in  the  next  paragraph  and  the 
one  developed  herein  do  not  experience  the  instability.  Further  analysis  of  this  process 
is  necessary  to  resolve  the  apparent  paradox,  but  this  issue  is  not  considered  within  the 
scope  of  the  present  research. 

Smith  et  al.  (1989)  have  devised  an  approximate  analytical  model  for  the 
/^-induced  gyre  structure  by  ignoring  vortex  propagation  and  computing  the 
streamfunction  associated  with  a  redistribution  of  the  initial  absolute  vorticity  by  the 
symmetric  TC.  The  inner  gyre  uniform  flow  predicted  by  the  model  closely  approxi¬ 
mates  the  propagation  speeds  predicted  by  nonlinear  numerical  models  up  to  about  24 
h.  Smith  (personal  communication,  1989)  is  presently  modifying  the  model  to  extend  the 
period  of  accurate  propagation  prediction  by  empirically  accounting  for  the  movement 
of  the  TC. 

b.  Application 

A  natural  application  for  a  barotropic  TC  propagation  model  would  be  as 
part  of  a  track  prediction  scheme  that  computes  total  TC  motion  at  any  time  by 
adding:  i)  a  steering  velocity  obtained  from  a  numerical  model  for  the  evolution  of  the 
large-scale  environmental  winds;  and  ii)  a  propagation  velocity  derived  from  the  theo¬ 
retical  model  based  on  the  gradient  of  absolute  (or  potential)  vorticity  of  the  numer¬ 
ically  generated  environmental  winds.  Early  efforts  based  on  linear  propagation  theory 
(Kasahara  1957;  Kasahara  and  Platzman  1963)  included  the  cross- vorticity  gradient 
propagation  effect,  but  ignored  the  hypothesized  up-vorticity  gradient  acceleration  of  the 
TC.  In  fact,  this  supposed  acceleration  due  to  a  net  force  on  "the  vortex"  has  never  been 
transformed  into  a  calculated  propagation  velocity.  A  recent  track  prediction  model 
developed  by  Holland  (19S3;  1984)  computes  westward  propagation  due  to  /?  similar  to 
the  earlier  efforts,  but  also  includes  a  northward  propagation  effect  by  empirically  ac¬ 
counting  for  net  convergence  into  the  TC.  That  these  models  have  shown  only  limited 
forecast  skill  is  not  surprising  in  light  of  the  recent  self-advection  theory  of  propagation, 
and  the  skill  that  they  do  possess  is  probably  due  to  some  portion  of  the  actual  nonlinear 
self-advection  phenomenon  being  parameterized  by  the  linear  models. 

The  present  unavailability  of  an  accurate  TC  propagation  model  based  on 
nonlinear  self-advection  has  precluded  any  implementation  into  an  "advection  +  prop¬ 
agation"  track  prediction  model  analogous  to  the  linear  models.  Although  the  present 
model  may  prove  suitable  for  this  type  of  endeavor,  development  and  testing  of  such  a 
model  is  outside  the  scope  of  this  research.  Unlike  the  linear  propagation  models  that 
arc  insensitive  to  small  deviations  of  the  TC  from  axisymmetry.  self-advection  models 


8 


depend  fundamentally  on  such  asymmetries,  and  more  importantly,  actually  predict  their 
structure  ( e.g .,  Willoughby  1988;  Smith  et  at.  1989).  Accurate  predictions  of  such  gyre 
structures  would  be  ideally  suited  for  initializing  barotropic  numerical  models  such  as 
SAXBAR  (Sanders  et  at.  1975)  so  that  quasi-steady  propagation  occurs  immediately  at 
the  start  of  the  integration  rather  than  after  the  24-48  h  typically  observed  for  an  initially 
symmetric  TC  (e.g.,  DM;  CW;  FE).  This  dissertation  will  conclude  with  a  proof-of- 
concept  demonstration  of  how  the  gyre/propagation  model  developed  here  can  be  used 
to  improve  numerical  model  initialization. 

B.  OBSERVATIONAL  EVIDENCE  FOR  PROPAGATION 

Given  the  many  barotropic  models  of  TC  propagation  in  the  literature,  it  is  natural 
to  ask  whether  observational  evidence  exists  to  support  the  predictions  of  the  theoretical 
or  numerical  models.  Unfortunately,  data  deficiencies  have  made  computation  of  prop¬ 
agation  for  an  individual  TC  difficult.  However,  studies  that  average  the  data  from 
many  TC's  to  produce  a  "composite"  TC  are  available  for  this  purpose.  The  propa¬ 
gation  vectors  in  Fig.  1.1  from  Carr  and  Elsberry  (1989)  are  based  on  data  from  two  such 
studies  (George  and  Gray  1976;  Chan  and  Gray  1982).  A  description  of  the  process  by 
which  the  vectors  were  obtained  is  given  in  Appendix  A. 

1.  Comparison  with  modeling  results 

The  propagation  vectors  in  Fig.  1.1  exhibit  a  number  of  interesting  properties 
that  strongly  resemble  the  TC  propagation  in  the  numerical  models  cited  above.  Except 
for  the  anomalous  "after  recurvaturc"  vector  (Fig.  1.10,  the  vectors  have  magnitudes 
ranging  from  1.0  to  2.5  m  s  and  directions  that  tend  to  be  westward  and  poleward  in 
both  hemispheres,  which  is  consistent  with  the  numerical  and  analytical  results  previ¬ 
ously  cited.  In  addition,  the  rotation  of  the  propagation  vector  direction  from  west- 
soutlnvestward  for  westward  moving  TC's  to  northwestward  for  eastward  moving  TC's 
in  the  direction  stratification  (Fig.  1.1b)  is  consistent  with  DM's  numerical  results.  DM 
showed  that  the  change  in  the  direction  of  the  environmental  vorticity  gradient  from 
poleward  on  the  poleward  side  of  the  subtropical  ridge  to  equatorward  on  the 
equatorward  side  caused  a  decrease  in  the  meridional  component  of  TC  propagation 
similar  to  that  in  Fig.  1.1b.  Finally,  the  propagation  vectors  in  the  intensity 
stratification  (Fig.  1  .Id)  have  a  direct  dependence  on  TC  intensity.  As  discussed  earlier, 
the  modeling  studies  of  DM  and  FE  demonstrate  that  a  nondivergent  barotropic  pre¬ 
diction  of  TC  propagation  due  to  /?  is  independent  of  TC  intensity,  but  is  well  correlated 
with  outer  wind  strength  (sec  Merrill  1984  for  typical  definitions  of  strength  and  inten- 
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Fig.  1.1  Vector  difTerences  of  composite  TC  motion  minus  steering  for  the  (a)  latitude, 
(b)  direction,  (c)  speed  and  (d)  intensity  stratifications  of  Chan  and  Gray  (1982),  and  the 
(e)  direction  and  (0  recurvature  stratifications  of  Holland  (1984).  The  vector  identifica¬ 
tion  labels  correspond  to  those  in  column  1  of  Tables  A.l  and  A. 2. 


sity).  Since  a  weak  correlation  exists  between  the  intensity  and  strength  of  TC's 
(Weatherford  and  Gray  1988),  the  increase  in  propagation  vector  magnitude  in  Fig.  l.ld 
may  be  a  manifestation  of  the  numerically-predicted  dependence  of  p- induced  propa¬ 
gation  on  TC  strength. 

The  results  in  Fig.  1.1  also  contain  some  apparent  inconsistencies.  Examples 
are:  i)  the  significantly  larger  meridional  components  of  the  Northern  Hemisphere  vec¬ 
tors  compared  to  that  of  the  Southern  Hemisphere  vectors;  and  ii)  the  presence  of 
equatorward  components  in  some  of  the  propagation  vectors.  Such  properties  may  be 
associated  with  boundary  layer  or  baroclinic  processes  not  considered  in  the  barotropic 
theories.  However,  a  possible  barotropic  explanation  might  be  the  presence  of  east-west 
vorticity  gradients  in  the  TC  environment  that  were  excluded  by  DM.  For  example,  a 
large-scale  westward  relative  vorticity  gradient  is  present  during  the  summer  in  the 
troposphere  between  the  anticyclone  over  the  western  North  Pacific  and  the  heat  low 
over  southeastern  Asia.  Based  on  DM's  results,  a  TC  vortex  embedded  in  such  a 
vorticity  gradient  should  have  a  westward,  and  more  importantly,  a  southward  compo¬ 
nent  of  propagation.  Since  the  meridional  gradient  of  environmental  relative  vorticity 
and  P  are  in  opposite  directions  south  of  the  Northern  Hemisphere  subtropical  ridge,  the 
zonal  gradient  of  environmental  vorticity  might  tend  to  dominate,  and  thus  explain  the 
southward  component  of  vector  W  in  Fig.  1.1b.  In  contrast,  the  meridional  environ¬ 
mental  vorticity  gradient  and  p  are  in  the  same  direction  north  of  the  Northern  Hemi¬ 
sphere  subtropical  ridge,  and  thus  might  dominate  over  the  influence  of  a  zonal  vorticity 
gradient  in  the  cases  of  vectors  N  and  E  in  Fig.  1.1b.  Differences  in  the  direction  of  the 
large-scale  absolute  vorticity  gradients  in  the  Northern  and  Southern  Hemispheres  thus 
may  contribute  to  the  hemispheric  variability  of  the  data  shown  here.  The  present 
propagation  model  will  permit  testing  of  such  a  hypothesis  within  a  modeling  context 
since  the  influence  of  an  east-west  environmental  vorticity  gradient  will  be  included. 

Statistical  influences  in  Fig.  1.1  also  must  be  considered,  such  as:  i)  ambiguities 
introduced  by  composite  stratifications  that  may  incorporate  multiple  propagation- 
inducing  influences;  and  ii)  possible  random  or  systematic  errors  in  the  composite  data 
that  may  be  significant  relative  to  the  small  size  of  the  propagation  vectors  being  ana¬ 
lyzed  (e.g.,  accuracy  of  rawinsondc  wind  measurements  and  error  in  locating  the  TC 
center  position).  Random  and  systematic  errors  should  be  reduced  as  sample  sizes  are 
increased  and  observational  accuracies  arc  improved  respectively. 
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2.  Interpretation  of  observed  propagation  vectors 

The  fundamentally  important  issue  that  determines  the  mechanism  responsible 
for  observed  TC  propagation  is  whether  or  not  a  TC-interaction  flow  (V„)  will  be  man¬ 
ifested  as  an  essentially  inseparable  contribution  to  conventionally  computed 
environnmental  steering.  If  VA  is  included  within  the  steering,  then  the  observed  TC 
propagation  vectors  must  be  due  to  some  mechanism  other  than  self-advection,  at  least 
within  the  context  of  the  barotropic  theory.  If  a  self-advection  flow  is  not  accounted  for 
in  the  computed  environmental  steering  flows,  then  the  observed  propagation  vectors 
may  indeed  represent  the  self-advection  process. 

Holland  (1983)  assumed  that  is  included  in  the  steering,  and  he  has  proposed 
a  linear  propagation  model  that  is  consistent  with  the  near-zonal  orientation  of  the 
Australian-Southwest  Pacific  difference  vectors  (Fig.  1.1  e-f).  Holland  assumes  that 
inertial  stability  constrains  the  inner  core  of  the  TC  to  move  with  the  outer  envelope, 
which  is  assumed  to  propagate  westward  as  a  Rossby  wave  with  a  phase  speed  appro¬ 
priate  to  an  "effective  radius."  In  practice,  the  "effective  radius"  parameter  for  a  partic¬ 
ular  storm  and  time  is  chosen  to  give  a  barotropic  Rossby  wave  propagation  speed  that 
equals  the  observed  westward  component  of  TC  propagation  over  a  preceding  time  in¬ 
terval.  Although  Holland  proposed  low-level  convergence  to  account  for  small  devi¬ 
ations  from  pure  westward  motion,  the  above  discussion  of  environmental  shear  in 
addition  to  />  is  an  alternative  explanation. 

In  contrast  to  Holland's  assumptions,  FE's  explicit  illustration  of  the 
wavenumber  one  structure  of  the  self-advection  flow  provides  strong  evidence  that  \A 
is,  for  the  most  part,  not  included  in  the  steering.  Key  aspects  of  this  structure  are  that 
the  central  uniform  flow  portion  of  the  gyres  that  account  for  the  propagation  is  con¬ 
fined  to  radii  less  than  300  km  from  the  TC  center,  and  that  the  gyre  flow  is  considerably 
weaker  outside  that  radius.  Even  allowing  for  somewhat  larger  gyres  in  actual  TC's,  this 
suggests  that  the  uniform  flow  region  of  the  interaction  flow  would  be  largely  unac¬ 
counted  for  in  a  steering  flow  calculated  over  an  annulus  of  5°-7°  lat.  radius  from  the 
TC  center  as  is  the  case  for  the  data  used  to  produce  Fig.  1.1.  Thus,  "self-advection" 
should  be  manifested  primarily  as  propagation,  rather  than  as  a  contribution  to  con¬ 
ventionally  calculated  steering  flows.  Thus,  Holland  s  linear  model  may  actually  include 
nonlinear  motion-inducing  processes  since  the  selection  of  an  effective  radius  is  based 
on  the  observed  difference  between  TC  motion  and  steering  over  a  preceding  time  in¬ 
terval. 
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Accurate  measurement  of  steering  near  the  TC  center  will  continue  to  be  prob¬ 
lematic  for  the  foreseeable  future.  Thus,  it  seems  advisable  as  a  practical  matter  to 
continue  to  compute  steering  at  a  large  scale  (*  1000  km)  as  is  now  done,  and  to  regard 
the  self-advection  flow  as  a  propagation-inducing  process  that  is  distinct  from  steering 
by  the  large-scale  environment.  In  light  of  the  subsynoptic  scale  region  in  which  the 
self-advective  flow  is  influential,  such  a  flow  partitioning  may  also  be  a  good  approxi¬ 
mation  to  the  long-sought  scheme  to  uniquely  separate  the  TC  from  its  environment. 
Thus,  this  research  will  use  the  three-part  partitioning  scheme  proposed  in  Elsberry 
(1986)  that  consists  of:  i)  a  symmetric  TC;  ii)  a  propagation-inducing  TC-environment 
interaction  flow;  and  iii)  a  large-scale  environmental  flow  component  that  accounts  for 
TC  steering. 

C.  DISSERTATION  OVERVIEW 

In  summary,  the  specific  objectives  of  this  research  are: 

1.  identification  and  analysis  of  the  specific  mechanism  responsible  for  barotropic 
vortex  stability  to  asymmetric  forcing,  and  how  this  stabilizing  influence  is  related 
to  TC  motion; 

2.  development  and  analysis  of  a  theoretical  model  for  barotropic  TC  propagation  due 
to  both  planetary  and  environmental  vorticity  effects  that  is  based  on  nonlinear 
self-advection,  and  that  accurately  predicts  gyre  structure  and  associated  propa¬ 
gation  velocity  relative  to  equivalent  numerical  model  solutions; 

3.  identification  of  the  dynamical  basis  for  the  dependence  of  barotropic  TC  propa¬ 
gation  of  outer  wind  strength,  and  an  assessment  of  the  wind  measurements  re¬ 
quired  to  account  for  this  phenomenon;  and 

4.  provide  a  preliminary  demonstration  of  the  viability  of  self-advection  propagation 
models  for  improved  initialization  of  barotropic  TC  track  forecasting  models. 

Because  of  the  mathematical  complexity  of  the  model  to  be  constructed,  Chapter  II 
will  be  devoted  to  preliminary  development  and  analysis.  Objectives  1  and  2  above  will 
be  addressed  in  Chapters  III  and  IV  respectively,  and  the  theoretical  aspect  of  Objective 
3  will  be  addressed  in  Chapter  IV  as  well.  In  Chapter  V,  the  model  will  be  applied  to 
satisfy  the  second  half  on  Objective  3  and  Objective  4.  Chapter  VI  will  conclude  the 
dissertation  with  an  overview  of  the  modeling  approach  and  a  summary  of  results. 
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II.  PRELIMINARY  MODEL  DEVELOPMENT 

A.  REFERENCE  FRAME  TRANSFORMATION 

The  behavior  of  fluid  motion  in  a  NDBT  system  on  a  /J-plane  is  governed  by  the 
conservation  of  absolute  vorticity,  which  may  be  expressed  as 

—  (£+/)  =  0  (2.1) 

£  =  k  •  V  x  V  (2.2) 

/=/*  +  Av,  (2.3) 

where  c  is  the  local  vertical  component  of  relative  vorticity,  /?  is  the  linearized  latitudinal 
derivative  of  the  Coriolis  parameter  /,  and  f0  is  the  average  value  of /  in  the  domain.  In 
terms  of  Eulerian  partial  derivatives,  (2.1)  may  be  expressed  as 

■jp  +  V  .  V(;  +/)  =  0.  (2.4) 

A  formal  partitioning  of  the  total  fluid  flow  that  takes  advantage  of  the  near- 
axisymmetry  of  TC's  as  suggested  by  (1.1)  and  (1.2)  requires  (2.4)  to  be  transformed  to 
a  reference  frame  moving  with  the  TC  at  a  translation  velocity  C (/).  In  terms  of 
cartesian  coordinates,  the  relationship  between  the  position  vectors  for  an  arbitrary 
point  P  in  each  reference  frame  is  illustrated  in  Fig.  2.1,  and  is  mathematically  defined 
by 


R( 


.r,  y)  =  R'(x',y')  +  [  C(t)</t , 


(2.5) 


where  the  (')  symbol  denotes  variables  in  the  moving  reference  frame.  Let  the  (')  symbol 
also  define  certain  derivative  operations  in  the  moving  reference  frame 


V'  EE  i 


ex 


y\?=*const 


+  i 


cy 


x'ti=consr 


(2.6) 


r 

Cl 


r 

cl 


X  ,  y  aeronsi 


(2.7) 
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I:ig.  2. 1  Relationship  between  the  position  vectors  of  an  arbitrary  point  P  in  a  reference 
frame  fixed  to  the  surface  of  the  earth  (unprimed  variables),  and  in  a  reference  frame 
moving  at  velocity  C(t)  (primed  variables). 
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By  using  the  chain  rule  and  (2.5)-(2.7),  it  may  shown  that 

V  =  V'  (2.8) 

—  =  —  -  C.V'  (2  9) 

ct  dt  ^  v  ’ 

which,  when  substituted  into  (2.4),  results  in  the  desired  transformed  expression 

(-fj-  -c +  V.V'(C  +/').  (2.10) 

It  is  important  to  note  that  a  convention  of  leaving  the  dependent  variables 
untransformed  in  (2.10)  has  been  employed.  Although  relative  vorticity  and  the  gradient 
of  the  Coriolis  parameter  are  invariant  with  respect  to  the  coordinate  transformation, 
fluid  velocity  is  not.  Thus,  the  velocity  transformation 

V-V'  +  C(f),  (2.11) 

obtained  from  d.  dt  of  (2.5)  has  not  been  used.  This  results  in  explicit  advection  by  C  in 
(2.10).  but  has  the  desirable  property  of  avoiding  an  implicit  dependence  of  the  fluid  flow- 
on  C  in  the  moving  reference  frame.  Such  a  convention  is  necessary  to  justify  the  use 
of  the  homogeneous  boundary  conditions  that  will  be  employed  in  solving  this  model, 
and  facilitates  later  use  of  the  solutions  in  a  numerical  model  initialization  scheme. 

B.  FLOW  PARTITIONING 

As  outlined  in  Chapter  I,  the  total  fluid  flow-  is  partitioned  into 

v  =  vs  +  V£  +  \A,  (2.12) 

according  to  the  following  definitions:  i)  Vs  is  a  known  symmetric  (S)  TC  flow'  com¬ 
ponent  that  is  steady  in  the  moving  reference  frame;  ii)  V£  is  a  known  environmental 
1 E j  flow;  and  iii)  VA  is  a  predominantly  asymmetric  (A)  flow  component  that  repre¬ 
sents  as  yet  unknown  interactions  between  the  symmetric  TC  and  the  specified  envi¬ 
ronmental  flow.  Subject  to  these  definitions,  substituting  (2.12)  into  (2.10)  gives 

(i  -  C-V'K +  w "  c’v'Cs  +  +  +  /,) 

+  (V£  +  V,).V'(:5  +  iA  +  :E  +  /')  =  0,  (M3) 
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in  which  the  term  Vs  •  V'{s  is  absent  since  Vs  is  by  definition  orthogonal  to  V'(s.  Re¬ 
moving  the  terms  associated  with  the  symmetric  TC  and  interaction  flows  from  (2.13) 
gives 

(-f-  -  C.V')f£  +  V£.V'((£  +/')  -  0,  (2.14) 

which  defines  the  evolution  of  the  environmental  winds  in  the  absence  of  the  TC  with 
respect  to  the  moving  reference  frame.  Subtracting  (2.14)  from  (2.13)  gives 

+  (Vj  +  V£  -  CJ.Vt,  +  (V,  -  Q.VCs 

+  V^.V(sT  +  {£+/)=  -V£.VCi-Vj.V«£+/),  (2.15) 

in  which  the  (')  symbols  have  been  omitted  for  simplicity.  In  all  subsequent  analysis, 
independent  variables  will  be  relative  to  the  moving  reference  frame.  The  various  terms 
of  (2.15)  have  been  grouped  so  that  the  left  side  represents  a  nonlinear  partial  differential 
equation  for  the  evolution  of  the  interaction  flow’  in  response  to  two  asymmetry- 
inducing  forcing  terms  on  the  right  side.  The  first  forcing  term  includes:  i)  distortion  of 
the  symmetric  TC  vorticity  field  due  to  large-scale  horizontal  w’indshear;  and  ii)  steering 
of  the  TC  by  the  environmental  winds.  The  second  forcing  term  represents  the  dispersive 
effect  of  Rossby  wave  radiation  associated  with  the  advection  of  large-scale  absolute 
vorticity  by  the  tangential  wind  of  the  symmetric  TC.  The  TC  translation  velocity  C  also 
appears  on  the  left  side  of  (2. 15),  and  represents  an  additional  unknown. 

The  effect  of  requiring  (2.14)  to  hold  with  the  TC  present  is  that  the  interaction  flow- 
must  represent  all  changes  to  both  the  symmetric  TC  and  the  “basic  state"  environment 
due  to  the  interaction  process.  However,  the  earlier  assumption  that  \s  may  be  regarded 
as  steady  in  the  moving  reference  frame  simplifies  the  interpretation  of  (2.15)  to  a  purely 
asymmetric  interaction  flow’  that  can  cause  TC  propagation  relative  to  the  large-scale 
environment  described  by  (2.14).  The  assumption  of  a  steady  TC  is  regarded  as  a  good 
first-order  approximation  to  lessen  the  considerable  complexity  of  the  TC  motion  prob¬ 
lem  evident  even  in  NDBT  dynamics.  However,  such  an  assumption  excludes  the  po¬ 
tential  impact  of  TC  strength  and  intensity  changes  from  the  scope  of  this  research. 

As  noted  above,  the  first  term  on  the  right  side  of  (2.15)  represents  both  steering  and 
shearing  of  the  TC  by  the  environmental  wind.  As  a  result,  the  TC  translation  velocity 
C  incorporates  both  a  known  advection  by  Vf  as  determined  by  (2.14).  and  an  unknown 


propagation  velocity  associated  with  the  interaction  flow  V„.  To  separate  the  known 
and  unknown  effects,  a  pure  shearing  component  of  the  environmental  wind  relative  to 
the  center  of  the  symmetric  TC  may  be  defined  as 

V£(r,  d,t)  =  V£(r,  6,1 )  -  V£(0,0,/) ,  (2.16) 

where  V^O.O,/)  represents  the  instantaneous  value  of  V£  at  the  TC  center.  Note  that 
V^O.O.r)  provides  a  theoretically  useful  definition  of  environmental  steering  for  this 

A 

model.  A  propagation  velocity  C  then  may  be  defined  by 

C(r)  =  C(t)  -  V^O.0,0 .  (2.17) 

Substituting  (2.16)  and  (2.17)  into  (2.15)  gives 

*s  y 

Sid.  +  (vs  +  v£  _  +  (V,  -  C).vcs 

+  Va'Wa  +  (£+/}=  -  V£.?Cs  -  V5.V(C£  +  /),  (2.18) 

in  which  the  unknowns  are  the  self-advection  flow  field  and  the  associated  propagation 

f 

velocity  C. 

It  should  be  noted  that  (2.14)  and  (2.18)  represent  the  basis  for  a  "propagation  + 
advcction"  track  prediction  model  based  on  self-advection  that  is  analogous  to  the  tra¬ 
ditional  linear  versions  discussed  earlier.  Equation  (2.14)  could  be  numerically  inte¬ 
grated  to  provide  V£(0,0,/)  at  appropriate  time  intervals,  while  a  solution  to  (2.18)  would 
provide  information  on  TC  propagation.  As  discussed  earlier,  this  research  will  focus  on 
solving  and  analyzing  several  approximate  forms  of  (2.18).  The  remainder  of  this  chap¬ 
ter  describes  the  additional  assumptions  and  concepts  on  which  the  solutions  in  Chap¬ 
ters  III  and  IV  are  predicated. 

C.  SIMPLIFYING  ASSUMPTIONS 
I.  Adjustment  time  concept 

Except  for  special  circumstances,  (2.14)  indicates  that  the  environmental 
windficld  will  evolve  with  time.  The  resultant  temporal  variability  of  V£  and  C£  makes 
(2.18)  analytically  intractable  for  use  as  a  theoretical  propagation  model  without  some 
simplifying  assumption  regarding  the  time  dimension  in  (2.18). 
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As  discussed  in  Chapter  I,  an  initially  symmetric  TC  in  a  NDBT  numerical 
model  will  propagate  in  response  to  steady  large-scale  forcing  at  a  nearly  steady  velocity 
and  with  quasi-steadv  asymmetric  structure  after  a  transient  adjustment  period  of  about 
24-48  h  (DM;  CW;  FE).  The  transient  phase  results  from  choosing  an  initial  TC  struc¬ 
ture  that  differs  from  the  steady-state  structure.  Since  TC's  are  continually  subjected 
to  environmental  forcing,  it  may  be  expected  that  mature  TC's  in  a  steady  or  slowly 
varying  large-scale  windfield  possess  an  asymmetric  structure  that  reflects  nearly  com¬ 
plete  adjustment  to  the  asymmetric  forcing  of  the  environment.  Thus,  the  temporal 

A 

varibility  of  V£  and  C£  in  (2.18)  will  be  eliminated  by  regarding  the  t  variable  in  (2.18) 
as  a  hypothetical  "adjustment  time"  during  which  a  "first-guess"  symmetric  TC  will  de¬ 
velop  a  quasi-steady  asymmetry  appropriate  to  the  forcing  defined  by  (2.14)  at  any  in¬ 
stant. 

In  a  rapidly  vary  ing  large-scale  environment,  the  TC  may  not  be  fully  adjusted 
to  the  asymmetric  forcing  of  the  environment,  such  as  during  TC  interaction  with  a 
rapidly  moving  midlatitude  trough.  Without  high  resolution  observational  data  to  ana¬ 
lyze  TC  propagation  and  horizontal  asymmetries  under  such  circumstances,  a  quantita¬ 
tive  estimate  for  the  amount  of  error  associated  with  a  steady-state  propagation  model 
is  not  possible.  Thus,  it  is  simply  noted  that  the  present  model  may  be  less  accurate  in 
such  circumstances. 

2.  Matched  solution  approach 

a.  The  transition  radius 

The  magnitudes  of  the  wind  speed  and  vorticity  of  the  symmetric  TC  com¬ 
ponent  are  much  larger  than  the  corresponding  asymmetric  and  environmental  compo¬ 
nents  over  a  significant  horizontal  area  in  the  lower  and  middle  troposphere.  Since  TC 
tangential  winds  decay  to  essentially  zero  at  large  radius,  a  region  must  exist  in  which 
the  magnitudes  of  the  asymmetric  interaction  and  environmental  flows  are  as  large  as 
the  symmetric  flow.  This  radius  is  called  the  "transition  radius,"  and  will  be  denoted  by 
Rr.  The  value  of  Rr  will  depend  on  the  structure  of  the  symmetric  TC  and  the  envi¬ 
ronmental  windfield  under  consideration,  and  would  vary  with  height  in  baroclinic  situ¬ 
ations.  Subject  to  some  additional  assumptions  delineated  in  Chapter  IV,  a  single 
transition  radius  will  be  assumed  to  exist  for  the  barotropic  model  used  here. 

b.  Self-advcction  Region 

For  r  <  R,.  it  is  assumed  that  advective  terms  in  (2. IS)  not  involving  the 
symmetric  flow  will  be  negligible.  Based  on  such  an  assumption.  (2.18)  simplifies  to 
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-jf-  +  vs.v^  +  <yA-  c).v:5  =  -  v£.vc5  -  vs.v(c£  +/).  (2.19) 

The  circular  area  within  which  (2.19)  applies  will  be  called  the  Self-advection  Region 
since  mutual  advection  by  the  asymmetric  interaction  flow  and  the  symmetric  TC  flow 
dominate  the  left  side  of  (2.19).  It  should  be  noted  that  the  Self-advection  Region  as¬ 
sumption  is  equivalent  to  regarding  the  interaction  flow  as  a  linear  perturbation  to  the 
symmetric  TC  flow. 

c.  Dispersion  Region 

For  r  >  Rr,  it  is  assumed  that  self-advection  terms  on  the  left  side  of  (2.18) 
will  be  negligible  relative  to  advection  terms  involving  the  environment  and  either  the 
symmetric  TC  or  interaction  flow.  Based  on  this  assumption,  (2.18)  simplifies  to 

+  Ve-VC,  +  V*.V(C£  +/)  =  -  V£.V£5  -  VS.V(C£  +/).  (2.20) 

The  outer  area  in  which  (2.20)  applies  will  be  called  the  Dispersion  Region  since  mutual 
advection  by  the  interaction  and  symmetric  TC  flows,  which  enables  barotropic  vortices 
to  resist  dispersion,  is  excluded. 

d.  Asymptotic  assumption 

In  an  annular  region  where  r  %  Rr,  the  assumptions  that  distinguish  the 
Self-advection  Region  and  Dispersion  Region  are  not  valid.  However,  (2.19)  and  (2.20) 
are  suitable  for  a  developing  a  matched  solution  for  a  propagation-inducing  interaction 
flow  field  that  becomes  asymptotically  close  to  an  exact  solution  to  (2.18)  as  r  becomes 
larger  or  smaller  than  Rr.  Such  a  matched-asymptotic-solution  approach  will  be  used 
in  Chapter  IV. 

3.  Isolating  wavenumber  one  processes 

The  diagnostic  analysis  by  FE  showed  that  the  interaction  flow  field  (their 
"ventilation  flow")  responsible  for  TC  propagation  due  to  just  /?  has  a  predominantly 
azimuthal  wavenumber  one  structure.  Chan  and  Williams  (1989)  have  shown  that  the 
interaction  flow  maintains  an  predominantly  wavenumber  one  structure  in  the  presence 
of  horizontal  wind  shear  in  the  environment,  even  though  forcing  at  higher 
wavenumbers  is  present.  The  absence  of  significant  wavenumber  two  and  higher  com¬ 
ponents  in  the  interaction  flow  will  be  addressed  in  the  stability  analysis  of  Chapter  III. 

The  observation  that  only  a  wavenumber  one  flow  structure  should  cause  TC 
propagation  is  intuitively  plausible  since  only  a  wavenumber  one  gyre  will  produce 
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nonzero  flow  across  the  center  of  the  symmetric  TC  (Fig.  2.2).  Formally  establishing 
this  concept  and  simplifying  governing  equations  for  the  Self-advection  and  Dispersion 
Regions  is  facilitated  by  writing  (2.19)  and  (2.20)  in  component  form.  In  a  polar  coor¬ 
dinate  system  centered  on  the  symmetric  TC, 

Vs=vs(r)e,  (2.21a) 

V*  =  +  vA(r,6,t)e,  (2.21*) 

C(/)  =  C(/)  [cos(0  -  a)  r  -  sin(0  -  a)  ©  ]  ,  (2.21c) 

A 

where  r  and  ©  are  unit  vectors  in  the  radial  and  azimuthal  directions  respectively,  and  C 
and  a  are  the  propagation  speed  and  direction  respectively.  The  mathematical  conven¬ 
tion  of  measuring  angles  counter-clockwise  from  east  has  been  used.  In  addition,  let 

V£  =  uE{r,d)i  +  £f(r,0)j,  (2.22) 


in  which  conventional  cartesian  symbology  has  been  used  for  the  environmental  wind 
components  and  unit  vectors.  Note  that  in  view  of  the  adjustment  time  concept,  (2.22) 
contains  no  dependence  on  time. 

a.  Wavenumber  one  Self-advection  Region  equation 
Substituting  (2.21a,b.c)  and  (2.22)  into  (2.19)  gives 


^  ^  y 

r>.1  v5  C~*A 

ct  '  CO 


+  \uA  —  C  cos  (6  —  a)] 


—  uE  - i s  cos  0  —  v£  ^  sin  0  —  v5[  — +  fi)  cos,  6  +  v5  sin  8  .  (2.23) 
ct  or  \  cj  j  cx 


c  *5 


'■E 


e¥ 


*E 


(a)  (b)  (c)  (d) 

Since  both  cCs/cr  and  C  are  independent  of  azimuth,  the  propagation  term 
on  the  left  side  of  (2.23)  is  a  wavenumber  one  process.  Only  projections  onto 
wavenumber  one  by  the  forcing  terms  on  the  right  side  of  (2.23)  can  cause  TC  propa¬ 
gation,  because  projections  onto  other  wavenumbers  are  orthogonal  to  the  propagation 
term.  The  validity  of  this  result  depends  on  the  assumption  that  the  term  C  •  may 
be  neglected  in  both  the  Self-advection  and  Dispersion  Regions.  Since  is  by  definition 
a  function  of  0,  projections  of  forcing  onto  other  wavenumbers  can  potentially  influence 
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Fig.  2.2  Hypothetical  streamfunction  patterns  (solid,  positive;  dashed,  negative)  for 
(lows  that  are  purely  azimuthal  wavenumber:  (a)  one;  (b)  two;  and  (c)  three.  Arrows  in¬ 
dicate  general  direction  of  the  flow  at  that  location. 


TC  propagation  via  this  term.  However,  FE's  finding  that  TC  propagation  was  almost 
completely  accounted  for  by  advection  by  the  wavenumber  one  gyre  provides  numerical 
support  for  regarding  the  interaction  flow  as  a  linear  perturbation  to  the  symmetric  TC. 

Simplifying  (2.23)  to  include  only  the  dominant  wavenumber  one  influence 
of  each  term  is  facilitated  by  writing  the  environmental  flow  factors  as  truncated  Taylor 
series  of  the  form 


A  A 

UE  .  v£  » 


=  A0  +  A,rcos(0  +  <£,). 


(2.24) 


Since  (2.16)  requires  A„  in  (2.24)  to  be  zero  for  uE  and  v£,  the  forcing  terms  (a)  and  (b) 
of  (2.23)  project  onto  wavenumbers  zero  and  two  to  first  order.  In  contrast,  A„  is  gen¬ 
erally  nonzero  for  the  first  derivatives  of  environmental  vorticity.  Thus,  to  first  order  the 
Self-ad vection  Region  equation  governing  TC  propagation  is 


Cl 
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Cl 


v5  e:A  a  dcs 

+  ~~dd~  +  l“A  ~  Ccos(0-a)]  — 


r-0 


cos  0  + 


sin  0  , 

r=0 


(2.25) 


which  may  be  expressed  more  concisely  as 


£k 

*> 

Cl 


+ 


2  Y 

r5  C  V/f 
r  c6 


dr 


+  iuA  -  C  cos(0  -  a)] 
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The  parameter  fi,ff  is  the  magnitude  of  the  environmental  absolute  vorticity  gradient  (i.e., 
an  "effective"  /?)  defined  by 


Peff  = 


(2.27) 


and  (j)  is  the  direction  associated  with  pfff  in  degrees  from  east  as  defined  by 
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b.  Wavenumber  one  Dispersion  Region  equation 

A  wavenumber  one  propagation  term  does  not  appear  explicitly  in  (2.20). 
Nevertheless,  it  may  be  anticipated  that  only  wavenumber  one  processes  in  the 
Dispersion  Region  are  important  to  TC  propagation  since  orthogonality  between 
wavenumbers  would  preclude  the  matching  of  higher  order  wavenumber  solutions  in  the 
Dispersion  Region  to  the  solutions  of  the  wavenumber  one  Self-advection  Region 
equation  (2.26).  In  addition,  FE  and  Chan  and  Williams  (1989)  have  shown  that  the 
self-advection  flow  is  predominantly  wavenumber  one  throughout  a  large  region  where 
the  Dispersion  Region  assumption  is  applicable. 

The  analysis  in  Section  C.3.a  above  has  already  shown  that  the  term 
V£.V£5  is  to  first  order  a  wavenumber  two  process,  and  by  implication  that  the  term 
V,,  •  V(Cf  +  /)  is  also  to  first  order  a  wavenumber  two  process.  Thus,  substituting 
(2.21a,b,c)  and  (2.22)  into  (2.20)  gives 
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=  v5jSe//sin(0  -  0).  (2.29) 


To  further  simplify  (2.29),  let  uE  and  v£  be  approximated  by  the  truncated  Taylor  series 
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in  which  sin0  and  cos 9  terms  have  been  used  in  lieu  of  the  phase  angle  $  in  (2.24).  In 
substituting  (2.30a,b)  into  (2.29),  only  terms  involving  sin:0  and  cosJ0  can  contribute  to 
wavenumber  one.  Retaining  only  the  wavenumber  one  components  of  such  terms  gives 
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where 
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(2.32a,  6) 


Noting  that  D£  and  Z£  represent  first-order  estimates  of  environmental  divergence  and 
vorticity  respectively  in  the  vicinity  of  the  TC,  and  assuming  that  D£xO  simplifies  (2.31) 
to 


,  ZLik 

dt  2  de 


vsPeff  sin(0  -  <f>) , 


(2.33) 


which  is  to  first  order  the  Dispersion  Region  equation  governing  TC  propagation.  Since 
Z£  2  represents  a  constant  angular  velocity  for  the  environmental  wind,  the  second  term 
on  the  left  side  of  (2.33)  represents  a  gyre-rotating  influence  consistent  with  the  numer¬ 
ical  results  of  Chan  and  Williams  (1989). 

D.  SUMMARY  AND  NOND1MENSIONALIZATION 

In  summary,  two  equations  have  been  developed  that  describe  to  first  order  the 
wavenumber  one  interaction  flow  associated  with  barotropic  TC  propagation.  A  trans¬ 
formation  to  a  polar  coordinate  system  moving  with  the  TC  and  a  partitioning  of  the 
total  flow  into  symmetric  TC,  large-scale  environment  and  a  TC-environment  interaction 
flow  have  been  employed.  A  transition  radius  Rr  has  been  defined  inside  of  which  is  a 
Self-advection  Region  and  outside  of  which  is  a  Dispersion  Region  where  mutual 
advections  by  the  symmetric  and  interaction  flows  are  important  and  unimportant  re¬ 
spectively.  The  wavenumber  one  interaction  flow  within  the  Self-advection  Region  is 
governed  to  first  order  by 
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and  within  the  Dispersion  Region  is  governed  to  first  order  by 
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(2.35) 


where  the  asterisk  subscripts  denote  dimensional  variables.  In  all  subsequent  analysis, 
the  absence  of  an  asterisk  will  denote  nondimensional  variables. 

This  preliminary  development  is  concluded  by  deriving  a  nondimensional  form  of 
(2.34)  and  (2.35)  by  scaling  r  by  the  radius  of  maximum  winds  (/?.„.),  all  velocities  by  the 
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maximum  symmetric  wind  ( VM.),  time  by  R^iV".  and  all  vorticities  by  V^:RM..  The  re¬ 
sulting  Self-advection  Region  equation  is 
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and  the  Dispersion  Region  equation  is 
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(2.38a, 6) 


III.  BAROTROPIC  VORTEX  STABILITY 


In  this  chapter,  the  specific  mechanism  responsible  for  barotropic  vortex  stability  to 
asymmetric  forcing  will  be  identified  and  analyzed.  Since  self-advection  has  been  clearly 
identified  to  be  necessary  for  this  stability  to  exist  in  barotropic  numerical  models  (C\V), 
the  scope  of  this  chapter  will  be  limited  to  the  Self-advection  Region  as  governed  by 
(2.32).  Despite  the  already  simplified  form  of  (2.32),  a  complete  solution  for  an  arbitrary 
initial  condition  that  includes  both  the  transient  and  quasi-steadystate  responses  in  TC 
propagation  models  is  in  all  likelihood  analytically  insolvable.  As  a  result,  the  focus  in 
this  chapter  will  primarily  be  on  the  transient  adjustment  process  by  transforming  (2.32) 
into  the  related  unforced  initial  value  problem.  Equation  (2.32)  includes  only  the  dom¬ 
inant  terms  describing  the  wavenumber  one  interaction  flow  associated  with  barotropic 
propagation  of  TC's.  To  explain  the  absence  of  higher  wavenumber  asymmetries  in  TC 
propagation,  the  following  development  will  address  vortex  stability  to  initial  asymmet¬ 
ric  perturbations  of  arbitrary  wavenumber. 

A.  MODEL  DEVELOPMENT 
1.  Background 

The  motivation  for  the  following  analytical  development  is  the  study  of  NDBT 
/-plane  Couette  flow  by  Case  (1960).  He  obtained  an  integral  solution  for  the  time  evo¬ 
lution  of  linear  perturbations  imposed  as  initial  conditions  on  the  NDBT  Couette  flow, 
which  is  a  steady,  zonally  uniform  flow  with  constant  latitudinal  shear  (Fig.  3.1a). 
Case's  result  may  be  interpreted  as  an  infinite  summation  of  a  continuum  of  singular 
solutions,  hereafter  referred  to  as  continuous  spectrum  modes  (Pedlosky  1964).  For  the 
NDBT  /-plane  Couette  flow  problem,  this  continuous  spectrum  forms  a  complete  basis, 
since  discrete  normal  modes  are  eliminated  by  the  lack  of  an  environmental  vorticity 
gradient.  The  superposition  of  these  continuous  spectrum  modes  results  in  a  perturba¬ 
tion  streamfunction  structure  that  has  an  algebraic  time  dependence  (i.e.,  depends  on 
factors  involving  t  to  an  integer  power)  as  the  perturbation  is  tilted  down-shear  by  the 
Couette  (low.  Although  the  baroclinic  Couette  flow  study  of  Farrell  (1982)  showed  that 
initial  growth  of  the  perturbation  is  possible  for  particular  initial  conditions,  both  Farrell 
and  Case  showed  that  the  response  is  asymptotically  proportional  to  r1  for  /-*  oo.  as 
schematically  portrayed  in  Fig.  3.1a.  In  terms  of  instability  theory.  NDBT  /-plane 
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Fig.  3.1  Schematic  portrayal  of  perturbation  damping  with  time  due  to  (a)  a 
meridionally  sheared  Couette  How,  and  (b)  a  radially  sheared  axisymmetric  vortical  flow. 
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Couette  flow  may  be  viewed  as  a  barotropically  stable  "basic  state"  with  respect  to  linear 
perturbations. 

The  NDBT /-plane  Couette  flow  model,  as  in  the  Eadv  or  Charney  models,  may 
be  viewed  as  an  idealization  of  the  response  of  synoptic  disturbances  to  a  particular 
planetary-scale  flow.  Along  with  other  limitations,  the  accuracy  of  this  approach  de¬ 
pends  on  the  extent  to  which  synoptic-scale  disturbances  may  be  regarded  as  linear 
perturbations.  In  the  case  of  an  intense  vortex  such  as  a  tropical  cyclone,  such  an  as¬ 
sumption  would  be  clearly  unjustified  since  the  TC  winds  can  be  substantially  stronger 
than  the  environmental  flow  at  quite  large  distances  from  the  center.  However,  the  large 
magnitude  and  nearly  circular  structure  of  the  TC  windfleld  suggests  that  a  model  anal¬ 
ogous  to  NDBT  /-plane  Couette  flow  could  be  developed  in  a  polar  coordinate  system 
moving  with  the  center  of  the  TC.  In  such  a  model,  the  axisymmetric  component  of  the 
TC  outside  the  radius  of  maximum  winds  ( RM .)  may  be  regarded  as  a  radially  sheared 
NDBT  "basic  state"  that  is  time-invariant  in  the  moving  reference  frame,  and  the  asym¬ 
metric  component  of  the  vortex  may  be  regarded  as  a  perturbation  to  the  symmetric 
basic  state.  If  the  axisymmetric  basic  state  has  constant  vorticity,  thch  an  initial  per¬ 
turbation  can  be  expected  to  damp  completely  as  it  is  tilted  in  the  direction  of  the  sym¬ 
metric  radial  shear  (Fig.  3.1b).  A  vortex  model  based  on  a  constant  vorticity  basic  state 
and  initial  perturbations  is  clearly  a  special  case,  and  the  implications  of  such  approxi¬ 
mations  will  be  addressed  explicitly  below. 

2.  Model  formulation 

The  analysis  of  Chapter  1  Section  C.3  showed  that  shearing  of  symmetric  TC 
vorticity  by  Vf  is  to  first  order  a  wavenumber  two  process.  It  is  desirable  to  include  this 
physically-based  process  in  the  present  analysis  by  rewriting  (2.32)  as 
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(3.1) 


where  the  first  two  terms  on  the  right  side  have  been  obtained  from  (2.23).  The  gov¬ 
erning  equation  for  the  unforced  initial  value  problem  related  to  (3.1)  is 
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C!  r  58 


A  dl  C 

+  [m^  —  C  cos(8  -  a)]  -rr~  =  0. 


(a)  (b)  (c) 


(3.2) 


Equation  (3.2)  describes  the  evolution  of  an  asymmetric  perturbation  initially  imposed 
on  a  symmetric  "basic  state"  vortex  in  a  quiescent  environment  and  on  an  /-plane.  With 
the  addition  of  two  homogeneous  boundary  conditions  at  specified  radii,  (3.2)  becomes 
the  homogenous  counterpart  to  (3.1)  over  the  enclosed  domain.  Since  these  equations 
are  linear  with  respect  to  the  asymmetric  perturbation  flow,  linear  differential  equation 
theory  can  be  employed  to  show  that  certain  relationships  exist  between  the  free- 
perturbation  response  and  the  forced-perturbation  response.  This  issue  will  be  specif¬ 
ically  addressed  in  Section  D.l  within  the  context  of  the  free-perturbation  behavior 
illustrated  in  Section  B  below. 

The  advection  of  perturbation  vorticity  by  the  radially  variable  symmetric  flow 
in  term  (b)  of  (3.2)  is  analogous  to  the  shearing  process  that  causes  perturbation  damp¬ 
ing  in  the  Couette  flow  model.  In  addition  to  manifesting  vortex  motion  through  the 
advection  of  symmetric  vorticity,  term  (c)  permits  the  propagation  of  neutral  and  possi¬ 
bly  exponentially  growing  discrete  normal  modes  on  the  radial  gradient  of  £s.  This  term 
represents  a  serious  obstacle  to  an  analytical  solution  approach  since  Cs  has  a  strong 
radial  dependence  that  varies  with  the  particular  vortex  structure.  To  facilitate  analysis 
of  the  damping  process  associated  with  term  (b),  term  (c)  will  be  removed  by:  i)  re¬ 
quiring  that  vs  has  a  Rankine  radial  dependence  over  a  horizontal  area  bounded  on  the 
inside  by  the  radius  of  maximum  winds  (RM.)\  and  ii)  limiting  the  domain  of  the  imposed 
perturbation  to  that  region.  To  facilitate  specification  of  initial  conditions  in  Subsection 
3  below,  the  modified  Rankine  profile  (see  Anthes  1982,  p.  22)  will  be  used 

fs  ~  r~X  1.  (3.3) 

The  Rankine  profile  is  the  limiting  value  of  (3.3)  as  X-*  I.  Substituting  (3.3)  into  (3.2) 
gives 
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where  os  represents  the  symmetric  angular  wind. 

Choosing  vs  to  be  Rankine  not  only  eliminates  any  discrete  normal  mode  com¬ 
ponents  from  solutions  to  (3.2),  but  also  may  have  a  significant  impact  on  the  contin¬ 
uous  spectrum  response.  Both  these  issues  must  be  addressed  before  the  results  of  this 
model  can  be  used  to  interpret  the  numerical  studies  cited  earlier  that  used  non-Rankine 
vortices.  The  impact  on  the  continuous  spectrum  can  be  assessed  by  drawing  an  analogy 
between  (3.2)  and  NDBT  Couette  flow  on  a  rotating  sphere  in  the  sense  that  the  radial 
gradient  of  has  an  influence  analogous  to  the  latitudinal  variation  of  the  Coriolis  pa¬ 
rameter.  For  a  /?-p lane  approximation,  Kao  (1955)  and  Boyd  (1983)  have  shown  that  P 
has  no  influence  on  the  rate  of  continuous  spectrum  damping,  but  rather  causes  the 
continuous  spectrum  wave  to  retrogress  proportional  to  P  in  a  manner  related  to  the 
familiar  discrete  mode  propagation.  Because  the  retrogression  is  independent  of  latitude 
for  a  constant  p,  the  damping  rate  depends  only  on  the  magnitude  of  the  Couette  flow 
shear.  Since  the  radial  gradient  of  (s  is  inward  in  the  inner  part  of  a  typical  vortical  flow 
and  decreases  rapidly  with  increasing  radius,  it  is  anticipated  that  a  continuous  spectrum 
wave  will  retrogress  (i'.e.,  clockwise  for  a  cyclonic  vortex)  faster  at  smaller  radii  than  it 
would  at  larger  radii.  The  result  would  be  slower  damping  since  the  retrogression  would 
tend  to  counter  the  tilting  induced  bv  vs.  This  effect  should  be  negligible  for  any  modified 
Rankine  vortical  flow  with  A'%1.  In  Section  C.2,  the  numerical  results  of  McCalpin 
(1987)  will  presented  to  show  that  significant  perturbation  retrogression  can  occur  for  a 
vortex  that  has  a  large  and  highly  variable  symmetric  vorticity  gradient.  In  addition, 
comparison  of  the  results  of  this  model  with  those  of  FE  and  DM  will  suggest  that  the 
symmetric  vorticity  gradient  of  a  tangential  wind  profile  that  approximates  a  TC  causes 
only  minor  slowing  of  perturbation  damping  due  to  variable  retrogression.  Finally, 
analysis  of  McCalpin's  results  will  also  suggest  that  virtually  no  energy  from  an  initial 
perturbation  projects  onto  any  discrete  modes.  Thus,  the  use  of  a  Rankine  vortex  will 
be  shown  to  be  reasonable  for  analysis  of  unforced  perturbation  responses. 

The  model  may  now  be  solved  by  defining  a  perturbation  streamfunction 

V^kx  CA  -  V2^  ,  (3.5a, b) 

and  assigning  homogeneous  boundary  conditions 

=  \jiA(b.0,i )  =  0,  (3.6a,/?) 
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where  a  =  1  and  b  must  be  chosen  to  confine  the  perturbation  to  the  Self-advection  Re¬ 
gion.  Although  the  analysis  of  Chapter  11  defined  the  outer  boundary  of  the  Self- 
advection  Region  to  be  a  specific  transition  radius,  here  a  convenient  choice  of  b=  10 
will  be  used.  The  justification  for  this  approach  is  that  the  Self-advection  Region  de¬ 
pends  on  perturbation  magnitude,  and  the  Self-advection  Region  in  an  unforced  model 
can  be  made  arbitrarily  large  by  choosing  sufficiently  small  perturbations.  The 
azimuthal  dependence  of  \I/A  will  also  be  expressed  in  terms  of  the  complex  Fourier  series 
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Substituting  (3.5b)  and  the  k,h  term  of  (3.7)  into  (3.4a)  and  integrating  with  respect  to 
time  gives 
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Substituting  (3.7)  into  (3.6a, b)  gives 


-  r:(M  =  0 , 


(3.8c,  d) 


where  C*  is  the  radial  structure  of  the  azimuthal  wavenumber  k  component  of  the  initial 
asymmetric  vorticity. 

liquations  (3.8a-d)  constitute  a  fully  specified  boundary  value  problem.  The 
solution  may  be  formally  written  in  terms  of  a  Green  s  function 
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(3-9) 


Using  standard  techniques  (Case  1960),  the  Green's  function  is  calculated  to  be 
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Finally,  substituting  (3.9)  into  (3.7)  gives 


iM'.M  =  Re  |  V, p)  C°(p)  e-,k{*~,COs(P))  dp  | ,  '  (3.11) 

which  represents  a  formal  solution  for  the  time  evolution  of  the  perturbation 
streamfunction.  Although  exact  solutions  to  (3.1 1)  do  exist  for  certain  initial  conditions, 
the  integral  generally  must  be  evaluated  numerically.  A  simple  trapezoidal  scheme  will 
be  used. 

3.  Initial  condition  specification 

To  evaluate  (3.11),  an  appropriate  functional  form  and  scale  for  must  be 
specified.  It  is  not  immediately  clear  what  combination  of  radial  and  azimuthal  depend¬ 
encies  to  associate  with  asymmetries  generated  by  the  interaction  of  a  TC  with  its 
environmment.  The  right  side  of  (3.1)  is  used  for  this  purpose  by  writing 

c*  cC  ^ 

-jr-  oc  vsfiejrjf$in(6  -  <f>)  -  uE-jf- cos  d  -  v£-^-sin0 

ci  cr  Cf 


+  cx-~-cosG  +  cy-~-  sin0 
A  cr  cr 


(3.12) 


A  A 

where  Cjf=C  cos  a  and  cr  s  C  sin  a  are  the  zonal  and  meridional  components  of  the  TC 
propagation  velocity  respectively.  Advection  by  C  has  been  included  from  the  left  side 
of  (3.1)  since  the  propagation  induced  by  uA  is  also  an  asymmetry-inducing  influence  (cf. 
Willoughby  19S8).  Since  asymmetries  associated  with  the  two  environmental  shearing 
terms  will  differ  only  with  respect  to  phase,  which  is  unimportant  in  the  unforced  re¬ 
sponse  analysis,  only  the  first  term  will  be  utilized  here.  Applying  similar  logic  to  the 
two  propagation  terms  in  (3.12)  results  in  the  simplified  form 
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(a)  (b)  (c) 


To  determine  the  spatial  structure  of  the  forcing  terms  in  (3.13),  the  radial  gra¬ 
dient  of  symmetric  vorticity  first  is  written  in  terms  of  (3.3).  Now  a  Rankinc  profile 
(.V  =  1 )  has  zero  vorticity  as  desired  to  delete  term  (c)  from  (3.2).  As  noted  in  the  Section 
A. 2  above,  a  small  perturbation  away  from  (A’=  1)  may  also  be  considered  as  a  valid 
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model  extension  since  the  radial  shear  of  vs  would  be  essentially  unchanged,  and  the  as¬ 
sociated  damping  should  be  little  affected  by  the  small  radial  gradient  of  that  would 
be  introduced.  Thus,  using  (3.3)  with  A'xl  gives 


cr 


A2-! 

3 


(3.14) 


Second,  let  uE  be  approximated  by  a  truncated  Taylor  series  expansion  about  the  present 
position  of  the  symmetric  cyclone  center 

uE{y)  =  SEy  =  SEr  sin  6  ,  (3.15a) 


^  A 


(3.156) 


where  S£  represents  the  nondimensionalized  linear  shear  component  of  the  environ¬ 
mental  wind,  and  any  zonal  dependence  of  uE  has  been  ignored  due  to  the  phase  argu¬ 
ment  mentioned  above. 

Considering  only  term  (a)  of  (3.13),  substituting  (3.3)  gives 


sin(0  -  4>) , 


(3.16) 


which  represents  the  generation  of  a  wavenumber  one  asymmetry  from  the  advection 
of  environmental  absolute  vorticity  by  the  symmetric  vortex.  Such  asymmetries  repre¬ 
sent  the  large  gyres  that  have  been  identified  in  numerical  models  (FE)  and  in  observa¬ 
tions  (Chan  19S6).  An  initial  condition  with  equivalent  spatial  structure  (hereafter 
referred  to  as  a  "vorticity-induced  asymmetry")  for  use  in  this  unforced  model  would  be 


Co(')  =  ~y~  k  =  1 


(3.17) 


V»rlRsr  ’ 


(3.18) 


where  a  poleward  gradient  of  absolute  vorticity  has  been  assumed  (  <p  =  90°),  and  ZA. 
represents  a  scale  for  the  perturbation  vorticity  that  remains  to  be  specified. 
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Considering  only  term  (b)  of  (3.13),  substituting  (3.14)  and  (3.15a)  gives 

c*  S 

oc  (1 -*2)— Tsin20,  _  (3.19) 

ct  y  ’  2 r2 

in  which  the  identity  2  sin  6  cos  0  =  sin  26  has  been  used.  The  right  side  of  (3.19)  re¬ 
presents  the  generation  of  a  wavenumber  two  asymmetry  due  to  linear  shearing  of  sym¬ 
metric  vorticity  by  the  TC-relative  environmental  wind.  An  initial  condition  with 
equivalent  spatial  structure  (hereafter  referred  to  as  a  "shear-induced  asymmetry")  for 
use  here  is 


C  *('-)=— T"  *-2,  (3.20) 

r 

where  it  is  has  been  assumed  that  SE  is  positive  corresponding  to  anticyclonic  shear 
equatorward  of  the  subtropical  ridge. 

Considering  only  term  (c)  of  (3. 13),  substituting  (3.14)  gives 

•>  y  ^ 

oc  (A"2  —  l)-y  cos  6 .  (3.21) 

ct  ri 

An  initial  condition  with  equivalent  spatial  structure  (hereafter  referred  to  as  a 
"motion-induced  asymmetry")  for  use  here  is 

£{r)  =  ~T  k  =  1  ,  (3.22) 

r 

where  it  is  assumed  that  cx  is  negative,  which  corresponds  to  the  generally  westward 
movement  of  TC  s  relative  to  environmental  steering  documented  by  composite  studies 
(Chan  and  Gray  1982;  Holland  1984).  In  both  (3.20)  and  (3.22),  it  has  been  assumed 
that  X  is  slightly  less  than  one  so  that  the  symmetric  TC  flow  has  cyclonic  vorticity. 

It  is  inappropriate  in  this  simple  model  to  assign  an  individual  perturbation 
vorticity  scale  ZA .  to  each  of  the  above  initial  conditions  based  on  the  scales  of  the  as¬ 
sociated  forcing  terms.  Instead,  all  but  one  of  the  following  model  solutions  will  use  a 
value  for  ZA .  such  that  c  =  0. 1  to  facilitate  analyzing  the  dependence  of  vortex  stability 
on  perturbation  spatial  structure.  The  radial  structures  of  the  four  initial  perturbations 
defined  above  are  shown  in  Fig.  3.2. 
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Fig.  3.2  Radial  dependence  of  initial  perturbation  vorticity  ( ZA .  C5 )  for  convection- 
induced  (k=2,  dot;  k  =  1,  dash),  motion-induced  (chaindot),  shear-induced  (chaindash) 
and  ^-induced  (solid)  asymmetries. 


36 


B.  MODEL  RESULTS 
1.  Perturbation  damping 

Before  numerically  evaluating  (3.11)  using  the  initial  conditions  developed 
above,  it  is  useful  to  obtain  exact  solutions  by  choosing  initial  conditions  of  the  form 


cSoo 


(3.23) 


An  asymmetry  of  this  type  may  be  viewed  as  an  approximation  to  the  wavenumber  k 
component  of  a  perturbation  produced  by  localized  asymmetric  convection  in  a  TC 
eyewall  cloud  (hereafter  refered  to  as  a  "convection-induced  asymmetry").  Initial  con¬ 
dition  (3.23)  for  k=  1  and  k  —  2  is  also  shown  in  Fig.  3.2,  and  results  in  solutions  to  (3.1 1) 
of 


0-0 


a2(r 2  —  b 2) 

V  -  b 2) 


cos(0  -  tja 2) 


4- 


r\b2  -  a2) 
r{a2  —  b2) 


cos(0  —  t/r2)  + 


b2{a 2  -  r2) 
r{a  -  b2) 


cos(0  —  ifb2) 


(3-24) 


and 


0-0 


+ 


+ 


4  l4\ 

a  (r  -  b  ) 

_L 

2,  4  .4. 

r  {a  -  b  ) 

2 

L  a 

4,. 4  is 

r  (b  —  a  ) 

’  j_ 

2/  4  i  4^ 

r  (a  -  b  ) 

2 

L  r 

I4-  4  “h 

b  (a  -r) 

1 

r  (a  -  b  ) 

.  b2 

cos  2(0  —  r/a2)  + 
cos  2(0  -  t/r2)  + 
cos  2(0  -  1/b1)  + 


sin  2(0  -  i/a 2) 


-Jf  sin  2(0  -  t/r2) 


(3.25) 


-JT  sin  2 (0-tlb2)  | 


respectively. 

Both  (3.24)  and  (3.25)  damp  algebraically  with  time  as  is  characteristic  of  con¬ 
tinuous  spectrum  mode  solutions.  It  should  be  noted  that  the  wavenumber  two  pertur¬ 
bation  damps  four  times  as  rapidly  as  the  wavenumber  one  perturbation,  which  suggests 
that  the  damping  rate  is  proportional  to  the  square  of  perturbation  wavenumber.  This 
result  is  analogous  to  the  dependence  of  damping  on  the  perturbation  zonal 
wavenumber  for  plane  Couette  flow  (Case  1960).  The  evolution  of  (3. 24)  and  (3.25)  over 
a  2  hour  period  are  shown  in  Fig.  3.3a-c  and  Fig.  3.3d-f  respectively.  Unless  otherwise 
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Fig.  3.3  (a)-(c)  Perturbation  strearnfunction  (solid,  positive;  dashed,  negative)  for  a  k  =  1 
convcction-induced  asymmetn’  at  t  =  0,  1  and  2  hours  respectively.  Contour  interval  is 
9.6 x  10 2m2ls.  (d)-(f)  Same  as  (a)-(c)  except  for  a  k=2  convection-induced  asymmetry, 
and  a  contour  interval  of  4.3  x  10J  »i2/s. 
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noted,  all  illustrated  solutions  use  the  following  parameter  specifications:  a=  1,  b=  10, 
Fv.=  50  m  s,  Rm.=  50  km  and  £  =  0.1.  The  contour  interval  is  the  same  for  each  row  of 
panels  to  show  clearly  the  damping  process.  This  comparison  emphasizes  the  strong 
dependence  of  the  damping  process  on  perturbation  wavenumber,  and  also  clearly  illus¬ 
trates  how  the  radial  variability  of  vs  tilts  the  perturbations  downshear. 

The  perturbation  streamfunctions  for  the  motion-induced  and  shear-induced 
asymmetries  (not  shown)  undergo  a  similar,  but  slower  shearing  and  damping  process. 
A  quantitative  comparison  of  damping  rates  is  facilitated  by  Fig.  3.4a,  which  shows  the 
decay  of  perturbation  kinetic  energy  with  time  for  each  of  the  four  cases.  Although 
perturbations  of  the  same  wavenumber  damp  at  nearly  the  same  rate  initially,  the 
damping  rates  quickly  diverge.  As  will  be  shown  in  Section  C.l,  this  property  is  due 
significant  differences  in  the  radial  dependences  of  the  initial  conditions,  e.g.,  r~6  for  a 
wavenumber  two  convection-induced  asymmetry  as  compared  to  r~ 2  for  the  shear- 
induced  asymmetry.  Thus,  NDBT  vortex  stability  to  asymmetric  perturbations  is  quite 
sensitive  to  the  radial  as  well  as  azimuthal  structure  of  the  perturbation. 

The  kinetic  energy  decay  for  a  vorticity-induced  asymmetry  is  shown  in  Fig. 
3.4b.  Perturbation  kinetic  energy  decays  by  approximately  80%  in  24  h,  which  agrees 
quite  well  with  the  adjustment  period  observed  in  the  NDBT  numerical  models  of  TC 
motion  cited  earlier.  The  much  slower  damping  rate  in  Fig.  3.4b  is  consistent  with  the 
discussion  in  the  preceding  paragraph  since  the  perturbation  vorticity  associated  with  a 
vorticity-induccd  asymmetry  decreases  significantly  slower  with  increasing  radius  than 
any  of  the  previous  initial  conditions  (Fig.  3.2).  The  evolution  of  both  the  total  and 
perturbation  streamfunction  components  over  a  36-hour  period  is  shown  in  Fig.  3.5. 
The  value  of  c  has  been  increased  to  0.5  in  Fig.  3.5  to  emphasize  how  an  initially  per¬ 
turbed  NDBT  vortex  is  restored  to  axisymmetry.  By  8  h,  the  inner  part  of  the  vortex 
has  already  regained  an  essentially  axisymmetric  structure  (Fig.  3.5b),  while  the  outer 
vortex  remains  appreciably  distorted.  No  analogy  to  this  behavior  exists  in  horizontal 
plane  Couctte  flow,  since  the  linear  shear  in  that  model  renders  the  damping  process 
independent  of  latitude. 

2.  Influence  of  boundary  conditions 

Boundary  conditions  (3.6a,b)  were  chosen  to  facilitate  development  of  the  model 
and  are  clearly  nonphysical.  Thus,  it  is  important  to  determine  the  extent  to  which  the 
boundary  conditions  may  have  influenced  the  results  presented  above.  Two  aspects  of 
this  problem  must  be  addressed. 
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Fig.  3.4  (a)  Decay  of  perturbation  kinetic  energy  (normalized  by  initial  value)  with  time 
for  convection-induced  (k  =  2,  dot;  k=  1,  dash),  motion-induced  (chaindot)  and  shear- 
induced  (chaindash)  asymmetries,  (b)  Same  as  (a),  except  for  a  /J-induced  (solid)  and  a 
modified  /i-induced  (dash)  asymmetry. 
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Fig.  3.5  (a)-(c)  Perturbation  plus  symmetric  streamlunction  for  a  ^-induced  asymmetry 
at  t  =  0,  8  and  36  hours  respectively.  Contour  interval  is  7.4  x  10s  (d)-(f)  Same  as 

(a)-(c)  except  showing  just  the  perturbation  streamlunction  (solid,  positive;  dashed, 
negative),  and  using  a  contour  interval  of  1,7  x  10s  mJ/s. 
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First,  boundary-  influences  will  become  increasingly  significant  with  time  if  there 
is  a  tendency  for  perturbation  energy  to  propagate  radially  and  cause  a  concentration 
at  either  boundary.  This  problem  may  be  formally  addressed  by  considering  - 

Ur,  0.0  =  Re{C*(r)c^-^>}  ,  (3.26) 

which  is  obtained  by  substituting  (3.5b)  and  (3.8a,b)  into  the  Laplacian  of  (3.7). 
Equation  (3.26)  indicates  that  perturbation  vorticity  is  conserved  following  the  symmet¬ 
ric  angular  wind,  which  is  a  purely  azimuthal  motion.  This  result  depends  on  the  model 
being  linear  and  the  symmetric  vortex  being  Rankine,  but  does  not  depend  on  boundary 
conditions.  Thus,  no  mechanism  exists  in  this  model  to  propagate  perturbation  energy 
radially,  nor  does  an  examination  of  Fig.  3.3  or  3.5  indicate  that  such  a  process  is  taking 
place. 

Second,  the  presence  of  the  inner  domain  boundary  at  the  location  of  maximum 
perturbation  vorticity  may  produce  nonphysical  results.  A  nonzero  inner  boundary  was 
used  because  a  Rankine  vortex  is  singular  at  the  origin.  The  inner  boundary  may  be 
moved  to  the  origin  if  the  singularity  is  removed  by  modifying  the  wind  profile  to  be 
solid  body  rotation  at  all  radii  less  than  Recall  that  the  choice  of  radial  structure 
for  the  various  initial  perturbations  depends  on  the  symmetric  vortex.  Thus,  the  radial 
structure  for  a  vorticity-induced  asymmetry  associated  with  a  Rankine  vortex  altered  as 
just  suggested  would  be 


{fa) 


0  <  r  <  1 
1  <  r  <  b 


(3.27) 


Fig.  3.4b  shows  that  the  kinetic  energy  for  such  an  initial  perturbation  decays  at  virtually 
the  same  rate  as  the  original  vorticity-induced  asymmetry. 

Modifying  the  symmetric  wind  as  described  above  causes  the  radial  gradient  of 
ls  to  be  singular  at  r=  1.  This  singularity  is  equivalent  to  a  large  radial  gradient  of 
vorticity,  which  would  tend  to  cause  clockwise  retrogression  and  somewhat  slower 
damping  of  the  continuous  spectrum  solution  as  discussed  earlier.  Although  the  present 
model  cannot  give  a  quantitatively  precise  solution  for  the  response  of  a  perturbation 
to  a  Rankine  vortex  with  solid  body  inner  rotation,  the  response  of  the  model  to  (3.27) 
is  useful  for  demonstrating  that  solutions  are  not  unduly  sensitive  to  the  location  of  the 
inner  boundary. 
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C.  MODEL  INTERPRETATION 
1.  The  stabilization  mechanism 

Within  the  context  of  fluid  dynamical  instability  theory',  the  present  stability  of 
a  linear  perturbation  with  respect  to  some  "basic  state"  may  be  assessed  by  computing 
the  domain-averaged  local  time  tendency  of  perturbation  energy.  For  the  particular 
boundary  conditions  used  in  this  model,  the  familiar  NDBT  result  is  that  domain- 
averaged  perturbation  kinetic  energy  must  be  presently  increasing  (decreasing)  with  time 
if  the  perturbation  tilts  against  (with)  the  horizontal  shear  of  the  basic  state.  Such  a 
result  is  readily  obtainable  for  this  model  by  substituting  (3.5b)  into  (3.4a),  multiplying 
by  \j/A,  integrating  over  the  domain,  and  performing  a  number  of  integrations  by  parts. 
Using  (3.5a),  the  result  may  be  expressed  as 


cE 


ct 


E  = 


(3.28a) 


(3.286) 


(3.28c) 


Equation  (3.28a)  is  analogous  to  the  familiar  result  for  stationary  cartesian  coordinates 
(e.g.,  Farrell  1987),  except  that  in  this  case  it  is  the  shear  of  the  axisymmetric  angular 
wind  that  controls  the  energy  transfer  process.  In  (3.28a),  a  positive  correlation  develops 
between  the  perturbation  momentum  llux  and  the  symmetric  angular  wind  shear  when 
initially  "upright"  perturbations  are  tilted  downshear.  Since  the  radial  shear  of  cj5  is 
negative  for  a  cyclone,  the  average  perturbation  momentum  flux  must  also  be  negative. 
This  flux  represents  an  inward  transport  of  cyclonic  momentum  that  tends  to  accelerate 
the  symmetric  "basic  state"  at  the  expense  of  asymmetric  perturbation  kinetic  energy. 

Because  the  analysis  thus  far  has  been  strictly  linear,  the  symmetric  vortex  has 
been  regarded  as  steady  in  the  moving  reference  frame.  However,  using  the  azimuthal 
momentum  equation  and  (3.28),  the  influence  of  perturbation  momentum  flux  on  the 
symmetric  basic  state  is  described  by 


C 


(3.29) 


ci 


=  -  (UaIa)- 


This  expression  may  be  evaluated  to  lowest  order  in  e  using  (3.11)  and  (3.26),  and  then 
integrated  with  respect  to  time  to  give 

Av5  =  77 &r)  C C0S  kl((°s(p)  ~  b)s{r))  -  (3JOa) 

Avs(r,i)  e  v5(r,r)  -  vs(r,0) .  (3.306) 

The  time  evolution  of  Ar$  is  illustrated  in  Fig.  3.6  using  the  same  initial  condition  as  in 
Fig.  3.5,  except  that  c  has  been  reduced  to  0.1,  which  corresponds  to  a  more  realistic 
asymmetric  windspeed  of  approximately  2.8  m,'s  near  the  domain  center.  Although  a 
nearly  complete  transfer  of  perturbation  kinetic  energy  occurs  by  48  h,  the  majority  of 
the  transfer  has  taken  place  by  24  h  as  anticipated  from  the  damping  rate  in  Fig.  3.4b. 
The  domain-averaged  sign  of  Ars  is  positive  as  required,  but  a  slight  decrease  is  evident 
for  r> 6.  Interestingly,  the  radius  at  which  Avs  changes  sign  corresponds  precisely  with 
the  maximum  in  perturbation  streamfunction  at  any  time  (see  Fig.  3.50.  It  is  also  evi¬ 
dent  that  the  energy  transfer  begins  initially  at  small  radii  and  spreads  outward  with 
time.  This  aspect  will  be  explained  later  in  this  subsection. 

The  momentum  flux  asssociated  with  the  convection-induced,  motion-induced 
and  shear-induced  asymmetries  (not  shown)  had  a  similar  impact  on  the  symmetric  basic 
state,  with  the  changes  in  Ais  becoming  smaller,  more  concentrated  at  small  radii  and 
occurring  over  a  shorter  period  of  time  for  an  initial  condition  that  decreases  more  rap¬ 
idly  with  radius.  The  extent  to  which  perturbation  flux  tends  to  alter  the  basic  state  is 
one  measure  of  model  linearity.  Thus,  the  results  here  show  that  the  assumptions  made 
in  Section  A. 3  (e  <,  0.1  ,  r  <.  10)  were  reasonable.  This  demonstration  of  how  perturba¬ 
tion  flux  transfers  energy  to  the  symmetric  vortex  will  be  particularly  useful  in  Section 
D.2  for  understanding  how  a  barotropic  vortex  achieves  a  quasi-steady  asymmetric 
structure  in  the  presence  of  steady  asymmetric  forcing. 

The  availability  of  exact  solutions  to  this  model  permits  valuable  insights  into 
the  local  dynamics  of  the  stability  process  that  are  not  readily  evident  from  the 
"domain-averaged"  analysis  conducted  above.  Note  that  by  virtue  of  the  identity 

V->±yw+^)]  =  0,  (3.31) 
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Fig.  3.6  Change  in  symmetric  vortex  windspeed  with  time  (see  legend)  as  a  function  of 
radius  due  to  a  convergence  of  momentum  flux  associated  with  a  jff-induced  asymmetry. 
The  windspeed  of  the  mean  symmetric  vortex  is  50  m/s  at  the  radius  of  maximum  winds 
and  decreases  with  inverse  dependence  on  radius  to  5  m/s  at  ten  times  the  radius  of 
maximum  winds. 
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the  first  and  third  terms  on  the  right  side  of  (3.24)  and  (3.25)  have  no  vorticity,  and  thus 
are  necessary  solely  to  satisfy  the  boundary  conditions.  Therefore,  substituting  either  of 
the  exact  solutions  into  (3.5b)  results  in  a  form 


Ur,  6,t )  =  V2  Re  {A(r,t)  elk(B~  “*'>}  ,  (3.32) 

where  A(r,t )  is  an  appropriate  amplitude  function  based  on  the  second  term  on  the  right 
side  of  either  (3.24)  or  (3.25).  Since  a>s  is  a  function  of  radius,  (3.32)  will  include  an  ex¬ 
plicit  time  dependence  in  some  of  the  terms  associated  with  radial  derivatives  of  the 
Laplacian.  For  sufficiently  large  t,  only  the  term  possessing  the  highest  order  explicit 
dependence  on  t  needs  to  be  retained,  which  gives 

Ur.M*  |  +  (t)2  (3.33) 

in  which  the  term  generated  by  the  azimuthal  derivatives  in  the  Laplacian  has  been  re¬ 
tained  for  comparison.  Note  that  the  expression  k/r  inside  the  brackets  represents  the 
inverse  of  the  perturbation  azimuthal  length  scale  at  radius  r.  If  the  first  term  in 
brackets  is  treated  similarly,  and  the  length  scales  are  defined  as 

L*  =  i  L*  -  —hr  •  (3'3Ja''’1 

kt  —r— 
cr 


The  local  dynamics  of  the  damping  process  may  be  explained  in  terms  of  (3.35) 
as  follows.  Under  the  influence  of  the  radial  shear  of  cos,  the  radial  length  scale  of  the 
perturbation  at  any  point  in  the  domain  exhibits  an  inverse  dependence  on  t,  while  the 
azimuthal  length  scale  remains  unaltered.  Since  this  model  conserves  perturbation 
vorticity,  the  left  side  of  (3.35)  maintains  a  constant  magnitude.  A  continuing  decrease 
in  Lr  thus  requires  the  amplitude  of  the  perturbation  streamfunction  to  decrease  pro¬ 
portional  to  r2  so  that  the  right  side  of  (3.35)  maintains  a  constant  magnitude.  This 
constraint  may  also  be  argued  conceptually  from  the  viewpoint  that  vorticity  is  a  vcloc- 
itv  chance  over  some  length  scale.  If  the  lencth  scale  is  decreasing,  then  the  velocitv 


change  must  also  decrease  to  conserve  vorticity.  This  process  is  shown  in  Fig.  3.7  which 
is  a  time  sequence  of  perturbation  vorticity  using  (3.26)  with  a  vorticity-induced  initial 
asymmetry.  The  reduction  in  the  radial  length  scale  due  to  the  shearing  process  is  clearly 
illustrated  by  the  decreasing  radial  spacing  of  the  isolines  of  vorticity  with  time. 

The  results  in  Section  B.l  show  that  vortex  stability  to  asymmetric  perturba¬ 
tions  is  strongly  dependent  on  the  spatial  structure  of  the  initial  perturbation.  From 
(3.34b),  this  behavior  can  be  associated  with  two  factors  that  influence  LR's  inverse  de¬ 
pendence  on  t.  The  first  is  the  perturbation  wavenumber  k,  which  results  in  perturba¬ 
tions  with  a  higher  wavenumber  damping  faster  than  those  with  a  low  wavenumber,  if 
the  radial  structures  are  similar.  This  is  the  case  for  the  wavenumber  one  and  two 
convection-induced  asymmetries  (Fig.  3.3),  which  differ  little  in  radial  dependence  (Fig. 
3.2).  The  second  factor  in  (3.34b)  is  the  radial  shear  of  the  symmetric  angular  wind, 
which  is  strongly  dependent  on  radius.  Thus,  the  shearing  process  reduces  LR  much 
more  rapidly  in  the  inner  part  of  the  vortex  relative  to  the  outer  regions,  which  explains 
why  the  perturbed  vortex  in  Fig.  3.5b  has  achieved  an  essentially  axisymmetric  state  in 
the  inner  region  while  the  outer  region  remains  distorted.  This  shearing  distribution  also 
explains  why  the  speed  of  the  damping  process  depends  on  the  radial  structure  of  the 
perturbation  vorticity.  An  initial  perturbation  that  decreases  in  magnitude  more  slowly 
with  increasing  radius  has  a  greater  fraction  of  its  kinetic  energy  at  larger  radii  than  docs 
a  perturbation  that  deci  cases  rapidly  with  radius.  Since  the  radial  shear  ofws  decreases 
rapidly  with  radius,  it  will  take  longer  to  transfer  the  same  amount  of  kinetic  energy  from 
the  first  perturbation.  Nevertheless,  energy  transfer  still  takes  place  rapidly  near  the  in¬ 
ner  boundary.  This  is  illustrated  in  Fig.  3.6  since  the  symmetric  vortex  windspeed  in¬ 
creases  almost  immediately  near  the  inner  boundary  and  is  followed  by  speed  increases 
at  larger  radii  with  increasing  time. 

2.  Verification  with  independent  numerical  results 

Insight  into  the  impact  of  using  a  Rankine  basic  state  vortex  in  this  model  can 
be  gained  by  comparing  the  present  results  with  a  similar  model  that  uses  a  non-Rankine 
vortex.  McCalpin's  ( 1 9 S 7 )  quasi-gcostrophic  reduced-gravity  numerical  model  of  an  in¬ 
itially  perturbed  ocean  eddy  will  be  used  for  this  purpose.  The  difference  between  the 
dynamical  frameworks  of  the  two  models  is  not  a  significant  issue  since  this  NDBT 
model  can  be  readily  reformulated  as  a  quasi-geostrophic  reduced-gravity  model.  The 
result  is  that  ZA  is  replaced  by  perturbation  potential  vorticity  in  (3.4a)  and  that  vs  be¬ 
comes  a  K,  exponential  Bessel  function  having  mo  potential  vorticity. 
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lTig.  3.7  (a)-(c)  Perturbation  vorticity  (solid,  cyclonic,  dashed,  anticyclonic)  for  a 

(l- induced  asymmetry  at  t  =  0,  1  and  2  hours  respectively.  Contour  interval  is 
9.1xl0-‘s-'. 
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The  main  difficulty  in  comparing  the  two  models  is  McCalpin's  choice  of 
Gaussian  radial  dependence  for  both  the  symmetric  and  asymmetric  components  of  the 
eddy.  Additionally,  McCalpin  characterized  the  damping  of  perturbation  energy  in 
terms  of  an  exponential  decay  time  scale,  whereas  the  continuous  spectrum  response  in 
this  model  is  algebraic  in  time.  Nevertheless,  if  the  convection-induced  asymmetries  used 
here  are  considered  to  be  roughly  equivalent  to  Gaussian  perturbations,  then 
McCalpin's  wavenumber  two  decay  time  scale  of  1.5  times  the  symmetric  flow  circu¬ 
lation  time  (at  the  radius  of  maximum  velocity)  may  be  compared  to  the  0.8  h  "expo¬ 
nential"  time  scale  given  by  Fig.  3.4a  for  a  wavenumber  two  convection-induced 
asymmetry-.  Using  the  1.75  h  circulation  time  for  at  r.=  in  this  model,  the 
damping  rate  here  is  approximately  3  times  faster  than  in  McCalpin's  model.  This  dif¬ 
ference  is  not  surprising  since  a  Gaussian  vortex  has  an  extremely  strong  and  rapidly 
decreasing  radial  gradient  of  symmetric  vorticity  just  outside  /?*.,  that  should  cause  a 
significant,  radially  variable  retrogression  that  can  substantially  reduce  the  rate  of  per¬ 
turbation  tilting  by  the  symmetric  flow.  This  assertion  is  supported  by  McCalpin's  ob¬ 
servation  that  perturbations  were  advected  around  the  eddy  at  only  20%  of  the 
maximum  tangential  velocity.  Tangential  winds  in  a  TC  decay  much  more  slowly  with 
increasing  radius  than  a  Gaussian  vortex,  and  thus  will  have  a  symmetric  vorticity  gra¬ 
dient  that  is  comparatively  much  smaller  and  decreases  more  slowly  with  radius.  Thus, 
it  may  be  reasonably  argued  that  the  damping  rates  of  this  model  are  more  represen¬ 
tative  of  those  associated  with  typical  TC  wind  profiles.  This  argument  is  also  suported 
by  the  good  agreement  between  the  24  h  adjustment  period  cited  in  NDBT  numerical 
studies  of  TC  motion  on  a  /?-plane  and  the  kinetic  energy  transfer  rates  in  Figs.  3.4b  and 
3.6.  With  regard  to  the  influence  of  perturbation  wavenumber,  McCalpin's  wavenumber 
three  decay  time  scale  is  2.6  times  shorter  than  the  wavenumber  two  time  scale,  which 
indicates  a  wavenumber  dependence  similar  to  that  above. 

Earlier  it  was  noted  that  choosing  Rankine  symmetric  vortex  excludes  poten¬ 
tially  important  discrete  normal  modes.  When  McCalpin's  model  was  run  on  an  /-plane, 
approximately  99%  of  the  energy  initially  present  in  the  imposed  perturbation  was 
transferred  to  the  symmetric  component.  Such  a  response  indicates  that  virtually  no 
perturbation  energy  was  projected  onto  any  discrete  normal  modes  that  might  exist  due 
to  the  presence  of  a  symmetric  vorticity  gradient.  Thus,  the  exclusion  of  discrete  mode 
processes  by  the  choice  of  a  Rankine  vortex  for  this  model  appears  justifiable  for  tran¬ 
sient  perturbation  analysis.  However,  this  does  not  rule  out  the  potential  importance 
of  discrete  modes  in  the  steady-state  component  of  a  forced  perturbation  response. 
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D.  MODEL  APPLICATION 

1.  Free  versus  forced  system  relationships 

The  motivation  for  the  present  modeling  approach  has  been  to  exploit  the  ana¬ 
lytical  tractabilitv  of  an  unforced  model  for  analyzing  transient  responses.  The  closed- 
form  solutions  obtained  via  this  technique  have  permitted  a  rigorous  and  illuminating 
analysis  of  the  perturbation  damping  process.  To  apply  these  results  to  the  forced 
problem,  it  is  necessary  to  establish  what  aspects  of  a  forced-perturbation  response  can 
be  reasonably  inferred  from  the  free-perturbation  response. 

First,  it  can  be  shown  by  superposition  that  the  full  solution  of  a  linear  partial 
differential  equation  with  steady  forcing  and  nonhomogeneous  boundary  conditions  can 
be  represented  as  a  sum  of  a  steady-state  part  that  satisfies  the  original  equation  with 
forcing  and  boundary  conditions,  and  a  transient  part  that  satisfies  the  homogeous 
equation  and  that  results  from  an  initial  condition  that  is  different  from  the  steady-state 
solution.  Implicit  in  such  a  partitioning  of  the  full  solution  is  the  assumption  that  the 
dynamical  system  actually  supports  a  non-trivial  steady-state  condition.  In  the  case  of 
vorticity-induced  asymmetries,  the  NDBT  numerical  study  of  FE  indicates  that  the 
asymmetric  component  of  the  vortex  does  tend  toward  a  steady-state  condition  over  the 
region  where  the  symmetric  vortex  is  significantly  stronger  than  the  asymmetric  com¬ 
ponent  (i.e.,  over  the  region  where  the  present  linear  model  is  valid).  Although  DM  did 
not  comment  on  vortex  asymmetries,  his  illustration  of  a  slowly  varying  TC  motion 
track  after  24  h  indicates  that  vortex  asymmetries  associated  with  steady,  spatially  vari¬ 
able  environmental  winds  also  tend  toward  a  quasi-steady-state.  Thus,  the  properties 
(e.g.,  damping  time  scales,  azimuthal  wavenumber  dependence,  etc.)  of  the  unforced 
transient  responses  shown  here  may  be  expected  to  be  relevant  to  the  initial  adjustment 
phase  of  a  symmetric  NDBT  vortex  on  a  /?-plane  with  steady  environmental  winds. 

Second,  in  a  linear  dynamical  system  that  is  first  order  in  time  and  exhibits  a 
damped  transient  response  to  steady  forcing,  the  magnitude  of  the  steady-state  condition 
can  be  expected  to  be  proportional  to  the  magnitude  of  the  forcing,  but  inversely  pro¬ 
portional  to  the  magnitude  of  any  parameter  that  acts  to  increase  the  rate  of  transient 
damping.  Such  an  assertion  is  formally  verifiable  in  the  case  of  a  constant  coefficient 
ordinary  differential  equation,  and  may  be  reasonably  extended  to  NDBT  vortex  dy¬ 
namics  for  situations  where  numerical  evidence  (i.e..  FE:  DM)  confirms  a  damped  tran¬ 
sient  response  to  steady  asymmetric  forcing.  Thus,  the  impact  of  perturbation  azimuthal 
and  radial  structure  on  the  damping  process  in  this  model  can  be  expected  to  influence 
the  steady-state  response  to  asymmetric  forcing. 
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2.  Barotropic  vortex  adjustment  to  steady  forcing 

The  process  by  which  vortex  adjustment  to  steady  asymmetric  forcing  occurs  in 
the  various  NDBT  simulations  of  TC  motion  (Anthes  and  Hoke  1975;  Kitade  1981;  DM; 
CW;  FE)  may  be  explained  as  follows.  As  the  initially  symmetric  vortex  advects  plane¬ 
tar}'  and  environmental  vorticity  or  is  distorted  by  environmental  wind  shear  and  vortex 
motion,  an  asymmetric  component  of  vortex  vorticity  is  generated  from  the  symmetric 
component  as  indicated  in  (3.12).  This  process  is  a  transfer  of  kinetic  energy  from  the 
symmetric  vortex  to  the  growing  asymmetry.  If  the  self-advection  process  is  omitted  as 
in  CW,  then  the  energy  transfer  continues  unchecked  and  causes  rapid  dispersion  of  the 
vortex.  However,  the  analysis  in  Section  C.l  showed  that  advection  of  the  asymmetric 
vortex  vorticity  by  the  radially  sheared  symmetric  angular  wind  acts  to  transfer  pertur¬ 
bation  energy  to  the  symmetric  vortex  (Fig.  3.6)  in  this  model.  In  the  \DBT  models  just 
cited,  the  rate  of  energy  transfer  to  the  symmetric  vortex  apparently  grows  as  the  forced 
asymmetry  grows  until  a  quasi-steady  balance  is  achieved  over  the  region  in  which  the 
symmetric  vorticity  is  significantly  greater  than  the  asymmetric  vorticity.  This  balance 
is  in  large  measure  achieved  by  about  24  h  for  asymmetries  associated  with  vortex 
advection  of  environmental  absolute  vorticity  (Fig.  3.4b),  and  occurs  much  faster  for 
asymmetries  associated  with  distortion  of  the  vortex  by  the  environment,  vortex  motion, 
or  asymmetric  convection  (Fig.  3.4a).  This  variability  in  adjustment  time  scale  for  dif¬ 
ferent  asymmetries  is  due  primarily  to  the  difference  in  radial  dependence  between  vs 
(/-1)  and  the  radial  gradient  of  Cs  (r3)  from  which  the  asymmetries  are  generated.  Al¬ 
though  these  particular  dependences  are  specific  to  a  near-Rankine  vortex,  the  principal 
applies  generally  since  the  symmetric  vorticity  gradient  must  decrease  faster  than  the 
symmetric  wind  for  all  vortical  flows  that  tend  toward  zero  with  increasing  radius. 

The  steady-state  phase  and  amplitude  toward  which  the  asymmetric  structure 
tends  clearly  cannot  be  addressed  with  a  homogeneous  model  that  also  excludes  the 
potentially  important  influence  of  discrete  normal  modes.  However,  two  aspects  of  the 
steady  state  may  reasonably  be  inferred.  First,  the  steady-state  asymmetry  will  likely  re¬ 
tain  a  down-shear  tilt  to  maintain  the  vortex  against  the  continuous  dispersive  effect  of 
the  asymmetric  forcing.  Such  a  feature  was  noted  by  FE  in  the  structure  of  a  vorticity- 
induccd  asymmetry.  The  second  inference  concerns  the  combination  of  radial  and 
azimuthal  dependence  that  is  likely  to  be  present  in  the  steady-state  asymmetry.  In  the 
absence  of  environmental  winds,  the  /J-effect  results  in  a  quasi-steady  vortex  asymmetry 
that  is  essentially  wavenumber  1  in  structure  (FE).  However,  as  shown  in  Chapter  II. 
a  horizontally  variable  environmental  windlield  will  act  to  induce  higher  wavenumbers 
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through  distortion  of  the  symmetric  vortex.  The  spatial  structure  of  the  resultant 
asymmetry  will  depend  on  both  the  spatial  structure  of  the  environmental  forcing  and 
the  dependence  of  the  damping  mechanism  on  perturbation  structure.  The  author  is  not 
aware  of  any  detailed  analyses  of  the  wavenumber  distribution  and  radial  structure  of 
TC  asymmetries  in  either  composite  observations  or  model  runs  using  realistic 
windfields.  Thus,  it  is  merely  noted  that  the  wavenumber  and  radial  dependences  of  the 
stabilization  mechanism  shown  here  should  contribute  to  the  predominance  of  low 
wavenumbers  and  increasing  axisymmetry  toward  the  vortex  center  respectively. 

E.  SUMMARY  AND  DISCUSSION 

The  principal  objective  of  this  chapter  has  been  to  identify  the  asymmetry-damping 
influence  of  symmetric  angular  windshear  as  the  mechanism  by  which  a  NDBT  vortex 
counters  dispersive  and  distorting  influences  over  the  region  dominated  by  nonlinear 
self-advection.  The  present  "linear"  model  captures  the  essence  of  the  self-advection 
process  by  linearizing  with  respect  to  a  nonzero  symmetric  basic  state.  In  the  NDBT 
models  cited  in  Section  D.2,  the  asymmetry-damping  mechanism  acts  as  a  negative 
feedback  process  in  which  a  kinetic  energy  transfer  from  the  asymmetric  to  the  sym¬ 
metric  component  of  the  vortex  occurs  as  a  result  of,  but  in  opposition  to,  the  kinetic 
energy  transfer  from  symmetric  to  asymmetric  component  induced  by  external  forcing. 

In  the  previous  section,  the  stabilization  mechanism  was  applied  to  the  initial  ad¬ 
justment  of  a  symmetric  NDBT  model  vortex  subjected  to  steady  asymmetric  forcing. 
However,  the  principle  can  also  explain  NDBT  vortex  adjustment  to  changes  in  the  ex¬ 
ternal  forcing  with  time.  For  example,  assume  that  a  quasi-steady  vortex  asymmetry 
exists  due  to  previously  steady  asymmetric  forcing,  and  that  a  change  now  occurs  in  the 
environmental  windfield.  e.g.,  in  fit/,  (3.16)  or  S£  (3.19).  If  the  change  is  such  that  the 
magnitude  of  asymmetric  forcing  at  a  particular  wavenumber  is  reduced  (increased),  then 
the  shear-induced  feedback  of  energy  from  the  existing  asymmetry  to  the  symmetric 
vortex  will  be  greater  (less)  than  the  environmentally-forced  transfer  of  energy  from  the 
symmetric  vortex  to  the  asymmetry.  As  a  result,  the  vortex  will  adjust  toward  a  less 
(more)  asymmetric  state  at  that  particular  wavenumber  until  a  quasi-steady  balance  is 
reestablished.  A  similar  adjustment  process  would  take  place  if  the  radial  distribution 
of  the  forcing  at  any  wavenumber  is  altered  by  changes  in  symmetric  vortex  structure. 
If  the  duration  of  the  external  forcing  change  is  brief  compared  to  the  time  scale  of  the 
stabilization  mechanism  (i.c.,  approximating  a  step-function),  then  the  adjustment  time 
should  be  on  the  order  of  the  stabilization  time  scale.  Conversely,  if  the  forcing  is  slowly 
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varying  in  time  compared  to  the  stabilization  time  scale  ( e.g .,  a  TC  propagating 
poleward  through  a  steady,  but  latitudinally  variable  windfield  as  in  DM),  then  the  ad¬ 
justment  process  should  have  a  time  scale  appropriate  to  the  "apparent  variability"  of 
the  environment  from  a  reference  frame  moving  with  the  vortex.  Such  a  scenario  is  not 
intended  to  be  all  inclusive,  since  dynamical  situations  may  exit  in  which  the  vortex 
might  be  barotropically  unstable  to  asymmetric  forcing,  or  in  which  temporary  contin¬ 
uous  spectrum  growth  might  occur  analogous  to  the  temporary  baroclinic  growth 
mechanism  studied  by  Farrell  (1982). 

Symmetric  angular  windshear  outside  the  radius  of  maximum  winds  can  be  expected 
to  exert  a  dominant  influence  on  the  stability  of  any  barotropic  vortex.  Shapiro  and 
Ooyama  (1989)  have  shown  that  divergence  had  negligible  influence  on  TC  propagation, 
and  thus  stability  to  asymmetric  forcing  within  a  barotropic  context.  In  baroclinic 
model  vortices  or  in  a  TC,  the  role  of  a  barotropic  stability  mechanism  in  influencing 
vortex  asymmetries  will  depend  on  the  competing  influence  of  the  inertial 
stability,  instability  made  possible  by  the  introduction  of  a  secondary  circulation  into  the 
dynamics.  Modeling  studies  of  ocean  eddies  indicate  that  coupling  between  vertical 
modes  can  also  be  expected  to  alter  vortex  stability  and  associated  motion  (McWilliams 
and  Flierl  1979).  Significant  radial  shear  in  the  tangential  winds  exists  to  large  heights 
in  a  mature  TC  (e.g.,  Hawkins  and  Imbembo  1976,  their  Fig.  13;  Frank  1977,  his  Fig. 
9).  Thus,  the  essential  element  that  enables  the  barotropic  vortex  stability  mechanism 
to  operate  is  certainly  present  in  TC's.  As  mentioned  earlier,  the  qualitative  similarity 
of  TC  motion  tracks  in  baroclinic  models  (Madala  and  Piacsek  1975;  Kitadc  1980)  to 
barotropic  results  provides  at  least  circumstantial  evidence  to  suggest  that  barochnic 
vortex  stability  is  a  modification  to,  rather  than  being  fundamentally  dilfcrent  from, 
barotropic  vortex  stability.  Since,  the  magnitude  of  TC  tangential  windshear  decreases 
with  height  above  the  boundary  layer  and  with  increasing  radius  outside  the  radius  of 
maximum  winds,  a  barotropic  stability  mechanism  should  be  most  influential  in  the  re¬ 
gion  where  the  TC's  convective  forcing  is  initiated. 


IV.  BAROTROPIC  VORTEX  SELF-ADVECTION 


In  this  chapter,  an  analytical  XDBT  model  of  TC  propagation  due  to  planetary  and 
environmental  influences  will  be  developed  from  the  Self-advection  Region  and 
Dispersion  Region  equations  of  Chapter  II.  The  three  principal  obstacles  that  must  be 
overcome  are:  i)  the  temporal  dependence  of  the  equations;  ii)  the  strong  radial  vari¬ 
ability  of  the  symmetric  flow  variables;  and  iii)  the  additional  unknowns  represented  by 
the  speed  and  direction  of  TC  propagation. 

The  problem  of  temporal  variability  will  be  eliminated  by  seeking  only  a  steady-state 
solution  for  TC  propagation.  Thus,  it  will  be  assumed  that  the  model  TC  is  fully  ad¬ 
justed  to  the  asymmetric  forcing  of  the  environment  at  any  time.  The  complex  radial 
variability  of  the  symmetric  flow  will  be  simplified  by  approximating  the  tangential  wind 
by  "piecewise-defined"  function.  The  presence  of  the  propagation  velocity  in  the  Self- 
advection  Region  equation  will  be  addressed  in  two  ways:  i)  using  externally  generated 
values  for  C  and  a;  and  ii)  devising  an  internal  closure  scheme  that  will  enable  the  pres¬ 
ent  model  to  predict  C  and  a. 

A.  MODEL  DEVELOPMENT 

1.  Solution  for  the  Self-advection  Region 

The  analysis  of  Chapter  III  has  shown  that  the  damping  influence  of  symmetric 
angular  windshear  will  permit  (2.36)  to  evolve  toward  a  steady-state  from  an  unbalanced 
initial  condition.  Thus,  the  analysis  here  will  focus  on  solving 
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The  propagation  term  appears  as  a  forcing  process  in  (4.1),  and  thus  introduces  the  un¬ 
known  parameters  C  and  a  as  noted  above.  The  streamfunction  and  complex  Fourier 
scries  definitions  of  Chapter  III  will  also  be  used  here.  Thus,  substituting  (3.5a,b)  and 
(3.7)  for  k  =  1  into  (4.1)  gives 
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in  which  the  identities 


cos  0  =  Re(c^}  sin  6  =  Re{  — /  e19] , 


(4.3a, b) 


have  been  used,  and  the  superscript  on  ¥  has  been  omitted  since  all  subsequent  analysis 

will  deal  exclusively  with  wavenumber  one  asymmetries.  Equation  (4.2)  is  put  in  stan- 
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dard  form  by  dividing  through  by  — ,  which  gives 
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To  obtain  closed-form  analytical  solutions  to  (4.4),  a  convenient  but  realistic 
functional  form  must  be  found  for  the  symmetric  flow.  In  particular,  requiring 


j__£k 

v5  cr 
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(4.5) 


has  the  desirable  property  of  making  the  left  side  of  (4.4)  equidimensional  similar  to 
(3.8a).  It  may  be  readily  verified  that 

v5('')  =  +  IV- X\  (4.6) 


has  an  associated  vorticity  gradient 
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and  that  the  product  of  (4.6)  and  (4.7)  satisfies  (4.5).  Thus,  substituting  (4.6)  and  (4.7) 
into  (4.4)  gives 
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T(r)  =  T,(r), 


(4.8a) 


iT(xr;-n 


(4.86) 


It  is  important  to  note  that  the  radial  gradient  of  the  symmetric  component  of 
the  vorticity  for  any  finite  vortex  such  as  a  TC  must  change  sign  at  least  once.  Subse¬ 
quent  calculations  of  Rr  will  show  that  a  vorticity  gradient  sign  change  occurs  within  the 
Self-advection  Region  for  TC  wind  profiles  typically  used  in  barotropic  models.  In 
contrast,  the  vorticity  gradient  given  by  (4.7)  is  monodirectional  for  a  single  value  of  X„ 
in  a  region  where  vs  is  purely  cyclonic.  Thus,  at  least  tw'o  segments  of  vs  as  defined  by 
(4.6)  for  different  parameters  (i.e.,  n=  1,2,  .  .)  must  be  combined  to  approximate  a  TC 
symmetric  windfield.  The  two  linear  (A„  and  B„)  and  one  nonlinear  (X„)  degrees  of 
freedom  in  (4.6)  permit  such  a  "piecewise-defined"  vortex  to  have  continuous  values  of 
velocity  and  vorticity,  but  a  only  piecewise-continous  vorticity  gradient.  As  a  result, 
(4.8a, b)  must  be  solved  in  at  least  two  annular  regions  for  particular  values  of  X„,  and 
then  a  total  solution  must  be  constructed  by  applying  matching  conditions  at  the 
interface(s)  as  shown  in  Subsections  A. 3  and  A.4  below.  The  approach  here  is  to  cir¬ 
cumvent  the  difficulty  in  obtaining  closed-form  analytic  solutions  (i.e.,  all  derivatives 
continuous)  to  (4.4),  by  seeking  piecewise-analytic  solutions  based  on  a  piecewise- 
analvtic  approximation  for  the  symmetric  TC. 

Consider  one  of  the  annular  regions  with  r=a  and  r=b  representing  the  inner 
and  outer  boundary  respectively.  Within  this  region,  a  solution  of  (4.8a)  may  be  ob¬ 
tained  by  noting  that  the  functions  r=x»~'  are  integrating  factors  for  the  left  side. 
Multiplying  (4.8a)  by  r^*1  and  integrating  from  a  to  r  gives 
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and  multiplying  (4.8)  by  r~x"*'  and  integrating  from  r  to  b  gives 
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in  which  a  n-subscript  has  been  added  to  T  to  denote  the  n’h  annulus.  Performing  the 
indicated  integrations  and  eliminating  dVJcr  from  the  equations  gives 
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in  which  G,  results  from  the  right  sides  of  (4.9)  and  (4.10),  and  is  defined  by 
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Solutions  for  as  many  annulae  as  necessary  to  approximate  the  symmetric  wind  profile 
with  segments  defined  by  (4.6)  may  be  constructed  from  (4.11). 

2.  Solution  for  the  Dispersion  Region 

Since  (2.37)  is  a  first-order  hyperbolic  equation  in  variable  CA,  the  method  of 
characteristics  may  be  used  to  give  an  initial  value  solution  of 


kM,/)  = 


^'sPeff 


Z  ft  1 

cos(0  —  — 2~  —  0)  -  cos(0  -  0)  , 


(4.12) 


in  which  CA{r,  6,0)  =  0  has  been  assumed.  The  response  is  an  undamped  oscillation  with 
period  47r/Z£  about  the  steady-state  solution  to  (2.37),  which  is 


2v5/Lr 
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(4.13) 


It  should  be  noted  that  no  steady  solution  exists  in  the  limit  as  Z£-*0,  and  that  (4.12) 
then  reduces  to 


tA(r,0.i)  =  tv5%sin(0  -  0) , 


(4.14) 


with  pe//=  /?  and  0  =  90°. 

In  view  of  abundant  numerical  evidence  for  quasi-steadv  TC  propagation,  the 
physical  relevance  of  (4.12)  and  (4.14)  for  large  t  is  doubtful.  Such  continually  evolving 
solutions  are  consistent  with  the  Dispersion  Region  assumption  that  completely  excludes 
the  stabilizing  influence  of  self-advcctivc  processes.  Within  the  present  modeling  con- 
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text,  asymmetry  damping  due  to  self-advection  can  impact  the  Dispersion  Region  only 
via  dynamical  matching  conditions  at  the  transition  radius.  Consequently,  the  evolution 
of  (4.12)  will  limited  by  determining  an  "adjustment  time"  (Te<(,)  such  that  the  Dispersion 
Region  solution  obtained  below'  and  the  solution  for  the  outermost  annulus  in  the  Self- 
advection  Region  satisfies  an  appropriate  matching  condition  at  the  transition  radius 
(Subsection  3). 

Making  the  substitution  and  applying  the  streamfunction  (3.5a, b)  and 

complex  Fourier  series  definitions  (3.7)  to  (4.12)  gives 

tt  +  tw  -  -r  w  =  F,w  (4.i5a) 
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Given  an  appropriate  form  for  vs  in  the  Dispersion  Region,  (4.15)  may  be  solved  using 
the  integrating  factors  r 1  and  t*.  A  convenient  choice  is 


vs('  )  =  +  B„r  5 , 


(4.16) 


since  A„  and  B„  provide  the  necessary  degrees  of  freedom  to  match  both  the  value  of  the 
function  and  the  first  derivative  of  (4.6)  and  (4.16)  at  r=  Rr.  Also,  the  choice  of  expo¬ 
nents  in  (4.16)  ensures  closed-form  solutions  to  (4.15)  that  remain  bounded  as  r  ->  oo. 
By  the  same  procedure  used  for  (4.11),  the  solution  to  (4.15)  is 
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Only  one  annulus  will  be  used  to  represent  the  Dispersion  Region.  Thus,  r  =  a  is  the 
transition  radius  and  r=  b  is  the  outer  boundary  of  the  model  in  (4.l7a,b).  In  the  case 
of  a  quiescent  environment,  evaluating  the  limit  of  (4.17b)  as  Z£  -*  0  gives 
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3.  Matching  annular  solutions 

The  solutions  presented  in  Section  B  below  will  use  at  most  three  annular  re¬ 
gions  to  define  the  model  domain.  Fig.  4.1  gives  a  schematic  view  of  the  model  domain 
and  annular  subregions,  and  also  shows  the  notation  to  be  used  to  denote  radial 
boundaries.  Limiting  values  for  R„  and  R3  are  the  origin  and  infinity  respectively.  With 
R2  =  Rr,  two  of  the  annulae  are  within  the  Self-advection  Region,  and  the  transition 
radius  is  a  matching  interface.  The  alternate  use  of  the  symbol  R2  for  the  transition  ra¬ 
dius  will  permit  a  concise  statement  of  the  matching  conditions  in  the  following  sub¬ 
section. 

a.  Transition  radius  specification 

A  scale  analysis  of  (2.36)  will  be  used  to  provide  an  empirical  formula  for 
Rr.  Implicit  in  the  transition  radius  concept  is  the  existence  of  an  annular  region  about 
Rr  where  interaction  and  symmetric  flow  variables  are  roughly  equal.  For  r2rRr .  let 

a  vv(Rr) 
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Substituting  these  relationships  into  (2.36)  and  assuming  a  quasi-steady  balance  gives 


vs 

r  cO 


v\s  ( R  7~) 
Rr 


C*S 


cr 


+  [uA  -  C  cos(0-a)]-^-  =  vspejy  sin{d  -  <f>) 

l's(Rr)  Peff 


vsi  R  r) 


R 


where  the  scales  are  shown  below  the  respective  terms.  Assuming  that  the  two  self- 
advection  terms  do  not  cancel,  it  is  necessary  that 
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Fig.  4. 1  An  illustration  of  the  model  domain  and  three  annular  subdomains.  The  inner 
and  outer  boundaries  of  the  model  domain  are  denoted  by  R*  and  R,  respectively.  The 
interface  between  the  inner  (n  =  1)  and  outer  (n=  2)  annulus  of  the  Self-advection  Region 
is  denoted  by  R,.  The  interface  between  the  Dispersion  Region  annulus  (n  =  3)  and  the 
inner  annulae  is  denoted  by  R2  and  corresponds  to  the  transition  radius. 
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from  which  the  empirical  formula 
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Pe/f 


(4.19) 


is  obtained.  Given  an  observed  or  specified  radial  profile  of  TC  tangential  winds  and  the 
magnitude  of  the  environmental  absolute  vorticity  gradient  in  the  vicinity  of  the  TC,  a 
transition  radius  may  be  calculated  from  (4.19). 

The  parameter  8  has  been  included  in  (4.19)  to  reflect  the  approximation  in 
(4.18),  as  well  as  the  uncertainty  that  will  accompany  any  estimate  of  Pt/r  Thus,  6  is  a 
"tunable"  factor  that  could  be  adjusted  in  an  "propagation  +  advection"  forecasting 
model  to  minimize  forecast  error  in  some  statistical  sense.  The  extent  to  which  a  sta¬ 
tistically  determined  d  approximates  unity  would  tend  to  provide  operational  confirma¬ 
tion  of  the  dynamical  validity  of  the  transition  radius  and  Self-advection,  Dispersion 
Region  concepts.  It  will  be  shown  below  that  assuming  8  =  1  provides  realistic  results, 
and  will  be  altered  only  to  test  model  sensitivity  in  Section  C.3  below. 
b.  Boundary  I  Interface  Condition  Specifications 

Eight  boundary  and  interface  conditions  must  be  specified  to  construct  a 
solution  for  the  entire  model  domain  from  solutions  to  (4.1  la)  or  (4.17a)  in  three  annular 
subregions.  The  asymmetric  streamfunction  will  be  set  to  be  zero  at  the  boundaries  of 
the  model  domain.  R0  and  R3,  and  the  asymmetric  streamfunction  and  velocity  will  be 
made  continuous  at  the  interface  radii  R,  and  R2.  The  inner  boundary  condition  is  dy¬ 
namically  required  if  R0  is  the  origin.  The  rationale  for  displacing  Ro  from  the  origin  is 
related  to  functional  properties  of  (4.6),  as  will  be  explained  later. 

These  conditions  may  be  concisely  expressed  in  terms  of  the  complex  vari¬ 
able  T  as 
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Practical  implementation  of  conditions  (4.20a-d)  involves  evaluating  (4.11a,b)  and 
(4.17a,b)  at  radii  R„  R2,  and  R3  to  form  a  6x6  linear  system  of  equations  in  the  variables 
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in  which  double  subscripts  denote  a  matched  quantity.  The  linear  equation  system  will 
be  solved  using  a  standard  computer  algorithm. 

Because  (4.8a)  and  (4.15a)  are  second  order  in  ¥,  the  associated 
boundary,  interface  conditions  provide  only  enough  degrees  of  freedom  to  achieve  con¬ 
tinuous  asymmetric  streamfunction  and  velocity  at  R,  and  R2.  Thus,  is  free  to  be 
piecewise-continuous  at  the  interface  radii,  which  is  generally  the  case  to  dynamically 
balance  (4.1)  given  the  piecewise-continuous  nature  of  the  vorticity  gradient  associated 
with  (4.6)  and  (4.16).  As  discussed  above,  the  parameter  T^  provides  an  additional  de¬ 
gree  of  freedom  to  permit  matching  of  asymmetric  vorticity  at  Rr.  In  general,  it  will  be 
possible  to  match  only  the  magnitude  of  CA  at  Rr,  since  varying  TtdJ  only  influences  G2, 
and  because  the  interface  condition  term  in  (4.17a)  is  irrotational.  Nevertheless,  the 
additional  degree  of  smoothness  given  to  model  solution  at  Rr  should  be  beneficial.  To 
be  consistent  with  this  quasi-matching  of  £*  at  Rr,  parameters  in  (4.6)  and  (4.16)  should 
be  chosen  to  eliminate  any  discontinuity  in  o£s/ct  there.  Not  only  is  it  possible  to  do 
this,  but  it  is  also  desirable  to  facilitate  closure  of  the  model. 

4.  Outline  of  Solution  Procedure 

A  NDBT  solution  for  the  propagation-inducing  wavenumber  one  gyre 
streamfunction  associated  with  known  symmetric  TC  and  environmental  windfields  has 
been  obtained  for  an  arbitrarily  large  domain  R<,  <  r  <>  R3  by  matching  three  partial  sol¬ 
utions  that  arc  separately  valid  in  annular  subdomains.  In  terms  of  the  real 
streamfunction  variable  *}/A,  the  solution  is 

e)  -  Rej^oy*}  Rn_,<r^Rn  ft-  1,2,3  (4.21) 

where  T*,  and  T2  are  defined  by  (4.1  la,b)  and  Tj  is  defined  by  (4.17a,b).  The  procedure 
to  compute  the  model  results  shown  in  the  remainder  of  this  dissertation  is: 

1.  Given  an  analytic  function  that  approximates  a  TC  symmetric  wind  profile,  deter¬ 
mine  R,  and  R2  (  =  R2)  based  on  where  the  symmetric  vorticity  gradient  changes 
sign  and  where  the  TC  winds  satisfy  (4.19)  respectively.  Also,  determine  values  for 
Rm.  and  1',,.. 

2.  Determine  parameters  A„,  B„  and  X,  by  constructing  from  (4.6)  and  (4.16)  a 
piccewise-analytic  approximation  to  the  wind  profile  used  in  Step  1.  Because  of  the 
irrational  function  form  of  (4.6),  this  task  is  accomplished  by  successive  trials  using 
interactive  computer  graphics.  The  procedure  is  outlined  in  Appendix  B. 

3.  Assign  values  to  the  environmental  parameters  /?,,,,  0  and  Zr. 
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A 

4.  Determine  propagation  parameters  C  and  a  using  either  external  information,  or  a 
closure  scheme  internal  to  the  dynamics  of  this  model. 

5.  Solve  the  8x$  linear  system  of  equations  for  the  value  and  first  derivative  of  the 
complex  streamfunction  as  determined  by  the  boundary/ interface  conditions 
(4.20a-d). 

6.  Iteratively  compute  the  complete  model  solution  (4.21)  by  varying  until  the 
magnitude  of  asymmetric  vorticity  is  continuous  at  r=  Rr. 

B.  MODEL  RESULTS  PART  I:  EXTERNAL  CLOSURE 

As  indicated  in  Chapter  I,  many  observations  of  TC  propagation  are  clearly  sug¬ 
gestive  of  the  barotropic  theories,  although  little  or  no  information  is  available  con¬ 
cerning  the  TC  characteristics  or  the  environmental  windfields  associated  with  those 
vectors.  Additional  problems  are  the  unknown  influence  of  baroclinic  processes  on  TC 
propagation  and  possible  random  systematic  errors  associated  with  the  compositing 
process.  Thus,  only  numerical  predictions  of  TC  propagation  are  suitable  for  externally 
closing  this  model  to  the  degree  of  accuracy  that  accurate  interpretation  and  analysis  of 
the  results  can  be  made.  The  XDBT  results  of  Fiorino  and  Elsberry  (FE)  (1989)  will  be 
used  since  that  study  provides  both  an  accurate  calculation  of  TC  propagation  with 
which  to  close  this  model,  and  detailed  illustrations  of  the  interaction  flow  from  which 
the  accuracy  of  closure  can  be  evaluated.  Since  FE  used  only  a  quiescent  environment, 
the  parameter  assignments  /?,,,=  /?,  </>  =  90°  and  Zf=0  will  apply  throughout  this  section. 
The  response  of  the  model  to  environmental  forcing  will  be  addressed  within  the  context 
of  the  internal  closure  scheme  developed  in  Section  C  below. 

1.  Symmetric  TC  Specification 

The  analytic  function  used  by  FE  to  approximate  the  symmetric  windfield  of  a 

TC  is 


v5(r)  =  re 


(1  -rbY,b 


(4.22) 


in  which  the  parameter  b  varies  the  strength  of  the  TC,  but  not  the  maximum  intensity. 
The  windspeed  and  vorticity  gradient  given  by  (4.22)  for  FE's  "basic  vortex"  (b  =  0.96) 
arc  shown  in  Fig.  4.2a  and  4.2b  respectively.  According  to  (4.18),  the  function  p,f/r2 
(Fig.  4.2a)  intersects  v s(r)  at  the  transition  radius  Rr=  4.75.  The  sign  change  of  vorticity 
slope  in  Fig.  4.2b  provides  a  good  first  estimate  for  R,.  The  location  of  the  inner 
boundary  R(J  depends  of  properties  of  the  piecewisc-analytic  wind  profile  as  described 
below. 
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VORT  GRAD  (VM./RM.2)  VELOCITY  (VM.) 


Fig.  4.2  (a)  Radial  profiles  of  the  analytic  TC  windHeld  for  the  Fiorino  and  Elsberry 
(1989)  basic  vortex  (dashed)  and  the  piccewise-analytic  function  used  (solid)  to  approx¬ 
imate  it.  Parameters  that  define  the  piecewise-analytic  function  according  to  (4.6)  and 
(4.16)  are  shown  in  the  inset.  Vertical  dotted  lines  (left  to  right)  correspond  to  the  radial 
boundaries/interfaces  R0,  R,  and  Rj  (=Rr).  The  chain-dashed  curve  represents  the 
function  ft,/rr}.  (b)  The  analytic  and  piecewise-analytic  vorticity  gradient  profiles  asso¬ 
ciated  with  the  wind  profiles  in  (a). 
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A  piecewise-analytic  wind  profile  is  defined  by  (4.6)  and  (4.16)  for  the  parame¬ 
ters  shown  in  the  inset  in  Fig.  4.2a.  Because  the  functional  properties  of  (4.6)  do  not 
permit  the  piecewise-analytic  wind  to  have  negative  curvature,  the  two  wind  curves  de¬ 
part  markedly  near  the  radius  of  maximum  winds.  As  a  result,  the  values  for  R*  and 
v5(R0)  have  been  chosen  to  improve  the  fit  of  the  two  curves  at  larger  radii  at  the  expense 
of  poor  fit  at  smali  radii,  which  is  justified  based  on  the  insensitivity  of  TC  propagation 
on  TC  intensity  demonstrated  by  FE.  It  should  be  noted  that  the  negative  curvature 
region  given  by  (4.22)  near  r  =  1  is  not  generally  representative  of  TC  winds,  which  are 
much  better  fit  by  modified  Rankine  profiles  just  outside  the  TC  eye  (Anthes  1982;  p. 
24).  Adjusting  the  piecewise-analytic  profile  to  accomodate  such  an  uncharacteristic 
wind  curvature  is  at  present  justified  by  the  need  to  compare  this  theoretical  model  with 
the  results  of  FE  that  were  based  on  (4.22). 

The  fit  of  the  corresponding  analytic  and  piecewise-analytic  vorticity  gradient 
profiles  is  illlustrated  in  Fig.  4.2b.  As  discussed  earlier,  the  piecewise-analytic  curve  has 
a  discontinuity  at  R,  due  to  the  change  in  X„  in  (4.7)  at  that  location.  In  contrast,  the 
piecewise-analytic  vorticity  gradient  is  nearly  continuous  at  the  transition  radius 
(Rr=R,)  in  anticipation  of  matching  the  magnitude  of  the  wavenumber  one  gyre 
vorticity  there.  The  choice  of  R,  =  3.5,  instead  of  where  the  analytic  vorticity  gradient 
changes  sign  (r=  3),  was  made  to  facilitate  the  matching  of  gyre  vorticity  at  Rr.  As  in 
Fig.  4.2a.  the  piecewise-analytic  curve  in  Fig.  4.2b  accurately  approximates  the  analytic 
profile  at  larger  radii  at  the  expense  of  a  good  fit  at  smaller  radii.  Because  the  vorticity 
gradient  depends  on  the  second  derivative  of  the  corresponding  windficld,  two  very 
similar  wind  profiles  can  have  quite  dissimilar  vorticity  gradients.  This  fact  has  impor¬ 
tant  implications  for  the  dependence  of  TC  propagation  on  outer  wind  strength,  and  will 
be  addressed  in  Section  C.3.b  of  this  chapter  and  in  Chapter  V. 

2.  /(-induced  gyre  structure 

Using  an  initial  symmetric  TC  wind  defined  by  (4.22)  with  b  =  0.96,  and  choos¬ 
ing  1’M.  =  35  m  s,  Rm.=  WO  km  and  /?.  =  2  x  10  "wr’-r',  FE  noted  a  /(-induced  propa¬ 
gation  velocity  of  2.65  m  s  at  330°  measured  clockwise  from  North.  Using  their 
parameter  specifications  and  the  piecewise-analytic  symmetric  TC  parameters  in  Fig. 
4.2a,  the  present  theoretical  model  for  gyre  structure  may  be  "closed"  by  requiring  C.  = 
2.65  m  s  and  a=  120°  (measured  counterclockwise  from  East).  The  resulting  asymmetric 
streamfunction  for  an  outer  model  boundary  at  10,000  km  (R3=  100)  is  shown  in  Fig. 
4.3  for  the  central  2400x2400  km  of  the  domain.  For  comparison,  the  gyre  structure 
extracted  from  FE's  numerical  mode!  after  48  h  is  shown  in  Fig.  4.4. 


X-DISTANCE  (Rm.) 

Fig.  4.3  Model-predicted  wavenumber  one  gyre  strcamfunction  (solid,  positive;  dashed, 
negative)  using  the  piecewise-defined  wind  profile  of  Fig.  4.2a,  and  the  parameter 
specifications:  FM.  =  40  m/s,  RKI. -  100  km,  /?.=  2x  10_"  nr's~',  4>  =90,  R0=  1.5,  R3  =  100, 

C.  =  2.65  m/s  and  a  =120.  The  contour  interval  is  2xl05/;j 2js,  and  only  the  inner 
2400x2400  of  the  domain  is  shown  to  correspond  to  the  illustration  from  FE  in  Fig.  4.4. 
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Fig.  4.4  Numerically  predicted  asymmetric  streamfunction  at  48  h  due  to  ^-induced  dis¬ 
tortion  of  a  TC  wind  profile  initially  defined  as  in  Fig.  4.2a.  Contour  interval  is 
2  x  lO'm-’/s,  and  the  distance  between  axis  tick  marks  is  40  km  (Fiorino  and  Elsberry 
1989). 
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The  analytically  and  numerically  generated  gyres  have  central  uniform  flow  re¬ 
gions  that  agree  quite  closely  with  regard  to  flow  strength  and  orientation,  although 
some  distortion  due  to  R0^O  is  evident  in  Fig.  4.3.  Such  agreement  is  of  fundamental 
importance  for  the  analytical  gyres  to  be  consistent  with  the  propagation  velocity  input 
from  the  FE  numerical  model.  The  radial  positions  of  the  streamfunction  extrema  are 
also  accurate.  However,  the  magnitudes  of  the  extrema  are  about  20%  less  in  Fig.  4.3. 
This  is  still  a  rather  remarkable  degree  of  agreement  considering  that  the  two  potentially 
adjustable  parameters  Rr  and  lc<tj  have  not  been  adjusted  to  improve  the  fit.  Thus,  the 
accuracy  of  the  analytical  gyre  of  Fig.  4.3  provides  preliminary  evidence  that  the  con¬ 
cepts  on  which  this  model  is  based  are  dynamically  sound.  Further  evidence  will  be 
provided  in  Chapter  V  by  showing  that  wavenumber  one  gyres  generated  from  this  an¬ 
alytical  model  can  accurately  initialize  a  numerical  model  so  that  steady  TC  propagation 
begins  almost  immediately. 

One  of  the  differences  between  Figs.  4.3  and  4.4  is  that  the  analytical  gyres  are 
symmetric  relative  to  a  line  drawn  between  the  extrema  and  antisymmetric  relative  to  the 
zero  contour  line,  whereas  the  numerical  solution  deviates  somewhat  from  these  condi¬ 
tions.  The  first-order  wavenumber  one  approximations  in  the  theoretical  model  and  the 
cyclic  east-west  boundary  conditions  used  in  the  numerical  model  are  likely  responsible 
for  this  dissimilarity.  By  retaining  only  the  wavenumber  one  processes,  the  theoretical 
analysis  has  excluded  the  interaction  process 
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which  results  in  a  westward  propagation  of  the  interaction  flow  on  the  gradient  of 
planetary  vorticity.  Such  a  propagation  could  explain  the  stretching  weakening  of  the 
cyclonic  gyre  and  compressing  strengthening  of  the  anticyclonic  gyre  in  Fig.  4.4.  This 
process  is  in  effect  generating  a  wavenumber  two  component  from  the  wavenumber  one 
gyres.  The  cyclic  boundary  condition  in  the  numerical  model  artificially  enhances  this 
process  by  enabling  the  opposing  gyres  to  interact  in  the  eastern  portion  ol  the 
anticyclonic  gyre,  which  is  displaced  southward  relative  to  the  analytical  counterpart  in 
Fig.  4.3.  However,  FE  noted  essentially  negligible  impact  of  these  wavenumber  two 
processes  on  TC  propagation  out  to  144  hrs.  which  is  further  confirmation  that  a 
barotropic  theory  of  TC  propagation  needs  to  retain  only  wavenumber  one  processes. 


3.  Influence  of  boundary  conditions 

As  noted  above,  the  choice  of  a  nonzero  inner  boundary  for  the  model  has 
perturbed  the  uniform-flow  region  of  the  analytical  gyre  (Fig.  4.3).  It  is  desirable  to 
eliminate  such  an  unphvsical  disturbance  by  moving  R0  to  the  origin.  For  a  wind  profile 
defined  by  (4.6),  it  is  not  possible  to  require  the  associated  vorticity  to  be  continuous  for 
all  r>0  and  also  be  nonsingular  at  r=0.  As  in  a  geostrophic  point  vortex,  the  velocity 
in  (4.6)  and  the  higher  derivatives  are  singular  at  the  origin.  Moving  Rq  to  the  origin  can 
be  justified  by  noting  that  it  is  the  ratio  of  the  symmetric  vorticity  gradient  to  the  sym¬ 
metric  velocity  that  influences  the  interaction  flow  within  the  Self-Advection  Region 
(4.4).  Since  the  ratio  of  (4.6)  and  (4.7)  contributes  to  a  regular  singular  point  at  the  or¬ 
igin.  the  interaction  flow  streamfunction  remains  defined  as  r-*  0,  even  though  %  as 
defined  by  (4.6)  is  singular  at  the  origin.  The  predicted  gyre  structure  for  R0^O  (Fig.  4.5) 
eliminates  the  disturbed  portion  of  the  uniform  flow  region,  but  leaves  the  rest  of  the 
gyre  structure  essentially  unchanged  from  that  in  Fig.  4.3.  Since  the  TC  associated  with 
Fig.  4.5  has  essentially  infinite  intensity,  this  result  is  >n  effect  the  ultimate  extension  of 
FE's  numerical  demonstration  that  barotropic  TC  propagation  is  virtually  independent 
of  TC  intensity.  The  inner  boundary  of  this  model  will  be  placed  at  the  origin  for  all 
subsequent  solutions. 

The  outer  boundary  of  the  theoretical  model  (R3=  100)  is  much  larger  than  that 
typically  used  by  numerical  models.  The  impact  of  choosing  R3  =  2S.  which  is  about  the 
distance  from  the  center  to  the  corners  of  the  4000x4000  km  domain  used  by  FE,  is 
shown  in  Fig.  4.6.  Although  the  peak  ulues  of  the  gyre  streamfunction  have  been  re¬ 
duced  about  8°o  compared  to  Fig.  4.5,  the  uniform  flow  region  of  Fig.  4.6  is  virtually 
unchanged,  which  confirms  the  domain  size  tests  of  Fiorino  (19S7).  Increasing  the  outer 
boundary  of  the  theoretical  model  had  no  significant  effect  on  the  gyre  structure,  which 
indicates  the  Ru=  100  is  effectively  infinity  as  far  as  the  theoretical  model  response  is 
concerned. 

C.  INTERNAL  CLOSURE  FORMULATION 

In  the  absence  of  external  information  on  TC  propagation  under  given  environ¬ 
mental  conditions,  determining  the  correct  wavenumber  gyre  structure  depends  on  a 
closure  scheme  to  predict  the  associated  propagation.  Such  a  scheme  will  be  termed 
"internal  closure"  since  it  is  based  solely  on  inherent  dynamical  characteristics  of  the 
theory.  The  analysis  that  follows  is  founded  on  the  basic  hypothesis  that  the  "correct" 
propagation  velocity  will  be  associated  with  a  particular  gyre  structure  that  is  unique  and 
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Fig.  4.5  As  in  Fig.  4.3,  except  for  R0s:0. 
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Fig.  4.6  As  in  Fig.  4.5,  except  for  R3=  28. 
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readily  discernable.  Since  the  work  of  FE  will  be  used  to  provide  both  motivation  for 
and  validation  of  this  approach,  the  closure  scheme  will  be  developed  within  the  context 
of  TC  propagation  in  a  quiescent  environment,  and  then  applied  in  Section  D  below  to 
illustrate  the  response  of  the  model  to  variable  environmental  winds. 

1.  Preliminary  Analysis 

As  noted  in  Chapter  I,  FE  found  that  the  advection  of  the  wavenumber  one  gyre 
by  the  symmetric  winds  of  the  TC  produces  a  streamfunction  tendency  that  nearly  can¬ 
cels  the  tendency  associated  with  the  advection  of  planetary  vorticity  by  the  symmetric 
winds.  Shapiro  and  Ooyama  (1989)  have  recently  refined  this  concept  further  by  show¬ 
ing  that  tangential  windshear  produces  nearly  uniform  asymmetric  absolute  vorticity 
within  2:350  km  of  the  TC  center.  As  shown  in  Chapter  II,  this  process  is  fundamentally 
a  kinetic  energy  feedback  loop  in  which  a  transfer  of  energy  from  the  asymmetry  to  the 
symmetric  TC  due  to  symmetric  angular  windshear  occurs  as  a  result  of,  and  in  oppo¬ 
sition  to,  the  transfer  of  energy  from  symmetric  TC  to  the  asymmetry  induced  by  a  gra¬ 
dient  of  planetary  vorticity  (/?). 

The  closure  scheme  developed  here  will  implement  the  above  concepts  by  seek¬ 
ing  a  particular  balance  among  the  various  processes  that  contribute  to  interaction  flow' 
vorticity  tendency.  Using  (4.1),  the  balance  in  the  Self-Advection  region  must  be 
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in  which  a  quiescent  environment  has  been  assumed,  and  the  identity  on  the  right  em¬ 
phasizes  the  steady-state  nature  of  the  present  model.  In  light  of  Shapiro  and  Ooyama's 
observation,  <4.24)  will  be  scaled  by  -  c^Jdr  to  give 
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which  is  valid  at  all  radii  except  where  the  symmetric  vorticity  gradient  goes  to  zero. 
Terms  (a)-(c)  of  (4.25)  represent  contributions  to  the  vorticity  tendency  of  the  inter¬ 
action  flow'  associated  with  symmetric  angular  windshear.  /?- induced  dispersion,  and 
advection  of  symmetric  TC  vorticity  by  the  uniform  portion  of  this  interaction  (low  re¬ 
spectively  (hereafter  the  "Shear",  "Beta"  and  "Advection"  terms).  'Jcrm  (d)  represents  a 
vorticity  tendency  that  is  an  artifact  of  the  movement  of  the  model  reference  frame  with 


the  propagating  TC  (hereafter  the  "Motion"  term).  The  quantity  within  the  brackets  in 
(4.25)  is  in  fact  the  asymmetric  absolute  vorticity  gradient,  which  according  to  Shapiro 
and  Ooyania  must  become  nearly  zero  inside  about  350  km  for  a  TC  undergoing  quasi¬ 
steady  propagation.  Thus,  the  cancellation  of  the  Shear  and  Motion  terms  results  in  the 
Motion  term  being  balanced  solely  by  the  Advection  term,  as  observed  by  FE. 

Since  all  the  terms  on  the  left  side  of  (4.25)  are  wavenumber  one  processes,  a 
particular  direction  may  be  associated  with  the  maximum  vorticity  tendency  of  each 
term  at  any  radius.  The  vectors  in  Fig.  4.7a, b  show  the  magnitude  and  direction  of  such 
vorticity  tendency  maxima  at  r=2  and  4  respectively  associated  with  the  wavenumber 
one  gyre  streamfunction  (Fig.  4.5)  predicted  by  this  model  using  a  piecewise-analytic 
approximation  to  FE  "basic"  vortex  (Fig.  4.2a)  and  their  reported  propagation  velocity 
of  2.65  m  s  at  330°  (a=  120).  Since  (4.25)  has  been  scaled  by  the  symmetric  vorticity 
gradient,  the  Shear,  Beta  and  Advection  terms  appear  in  Fig.  4.7a, b  as  velocity  compo¬ 
nents  that  must  together  balance  the  imposed  Motion  vector.  Since  the  symmetric 
vorticity  gradient  changes  sign  in  the  Self-advection  Region,  two  radii  have  been  chosen 
to  permit  comparison  of  the  balance  between  the  inner  (r=2)  and  outer  (r=4)  annulus 
of  the  Self-advection  Region. 

A  radial  profile  (Fig.  4.7c)  of  the  amplitude  and  phase  of  the  interaction  flow 
vorticity  (cj  corresponding  to  the  streamfunction  in  Fig.  4.5  highlights  distinct  ampli¬ 
tude  and  phase  changes  between  the  inner  and  outer  annulus  of  the  Self-advection  Re¬ 
gion.  Identifying  specific  patterns  in  both  this  profile  of  interaction  flow  vorticity  and 
the  orientation  of  vorticity  tendency  vectors  above  is  the  basis  for  the  following  closure 
scheme. 

2.  The  Closure  Scheme 

Because  the  Shear  and  Beta  vector  tend  to  oppose  each  other  in  Fig.  4.7a  and 
b,  the  Advection  vector  balances  most  of  the  imposed  Motion  vector.  Such  a  vector 
pattern  may  be  regarded  as  roughly  consistent  with  FE's  observation  that  the 
streamfunction  tendencies  associated  with  the  Shear  and  Beta  terms  of  (4.25)  nearly 
cancel  in  their  numerical  solutions.  If  barotropic  TC  propagation  is  indeed  precisely 
steady  propagation,  then  the  Beta  and  Shear  vector  would  be  expected  to  cancel  exactly. 
The  inexact  cancellation  in  Fig.  4.7  may  be  attributed  to  inaccuracies  in  the  theoretical 
model  and  or  the  possibility  that  barotropic  TC  propagation  is  more  a  quasi-steady 
process  than  a  truly  steady  process.  Assuming  that  the  second  factor  is  primarily  re¬ 
sponsible  for  the  vector  pattern  in  Fig.  4.7.  a  closure  scheme  based  on  precisely  steady- 
propagation  will  be  feasible  if:  i)  a  propagation  velocity  can  always  be  found  such  that 
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i  ,v  4.7  (a)-(b)  Vector  diagrams  showing  direction  and  magnitude  of  C<»  tendency 
maxima  that  have  been  scaled  by  -  6{sldr  at  (a)  r*2  and  (b)  r*4  respectively.  The 
vectors  represent  the  terms  of  (4.24)  (see  key  above),  and  were  computed  using  the 
piccewisc-defined  symmetric  1C  wind  profile  of  Fig.  4.2a  and  the  wavenumber  one  gyre 
solution  of  Fig.  4.5.  (c)  Radial  profiles  of  the  amplitude  (solid)  and  phase  (dashed)  of  (4 
corresponding  to  the  strcamfunction  field  of  Fig.  4.5.  The  amplitude  curve  has  been 
scaled  by  the  amplitude  of  £4  at  r=U.2.  Vertical  dotted  lines  (left  to  right)  correspond 
to  R,  and  Rj. 
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theoretical  model  predicts  precise  cancellation  of  the  Shear  and  Beta  vectors;  and  ii)  such 
a  propagation  velocity  is  a  reasonable  approximation  to  the  quasi-steady  result  obtained 
from  a  dynamically  equivalent  numerical  model.  As  shown  in  Fig.  4.8a,  the  propagation 
velocity  C=  1.9,  a  =  130  satisfies  the  first  of  the  above  conditions.  Such  a  closure  scheme 
is  clearly  radius  dependent  since  the  Shear  and  Beta  terms  in  Fig.  4.8b  are  opposite  in 
phase,  but  are  somewhat  unequal  in  magnitude.  Although  the  theoretically  predicted 
propagation  direction  differs  by  only  10°  from  the  propagation  direction  reported  by  FE, 
the  speed  difference  of  0.75  m  s  represents  a  significant  error. 

Both  the  speed  error  and  the  radius  dependent  nature  of  the  above  closure 
scheme  may  be  addressed  by  recalling  FE's  demonstration  that  TC  propagation  and 
advection  by  \A  in  the  uniform  flow  region  are  nearly,  but  not  precisely  equal.  They 
observed  vector  differences  between  the  propagation  velocity  of  the  TC  and  the  uniform 
portion  of  the  interaction  flow  (their  "ventilation  flow")  that  are  persistently  westward 
at  about  0.3  ms  (Fig.  4.9).  This  difference  is  readily  evident  quite  early  (12  h)  in  the 
integration,  which  suggests  a  linear  Rossby  wave-type  propagation  associated  with  the 
subsynoptic  scale  Fourier  components  of  the  symmetric  TC.  Such  a  hypothesis  is  sup¬ 
ported  by  Shapiro  and  Ooyama's  (1989)  Fourier-Bessel  spectral  analysis,  which  showed 
that  the  squared  angular  momentum  spectrum  for  a  typical  TC  wind  radial  profile  has 
a  distinct  peak  at  a  radial  length  scale  of  200  km.  If  this  radial  length  scale  is  regarded 
as  approximately  equivalent  to  the  inverse  of  a  sinusoidal  wavenumber,  then  it  may  be 
shown  that  a  symmetric  nondivergent  Rossby  wave  of  such  a  length  scale  will  propagate 
westward  at  0.4  m  s. 

The  requirement  that  the  Shear  vector  exactly  cancels  the  Beta  vector  in  the 
closure  scheme  is  inconsistent  with  these  numerical  results  that  only  the  component  of 
the  Beta  vector  over  about  0.4  m  s  actually  contributes  to  distortion  of  the  TC  circu¬ 
lation.  That  is,  only  the  relative  linear  propagation  speeds  of  the  large  and  small  scale 
Fourier  components  contributes  to  the  distortion  of  the  TC  circulation.  Thus,  a  better 
closure  formulation  is  to  choose  C.  and  a  so  that  at  a  radial  location  200  km  from  the 
TC  center  the  Beta  and  Shear  vectors  are  1S0°  out  of  phase  and  differ  in  magnitude  by 
0.4  m  s  (Fig.  4.10a).  The  resulting  values  for  C.  and  a  for  this  closure  scheme  are  2.4 
m  s  and  132r  respectively,  which  reduces  the  speed  error  to  0.25  m  s.  Such  an  error  is 
considered  quite  acceptable  in  view  of  the  typical  magnitudes  of  total  TC  motion 
(advection  +  propagation).  The  Beta  and  Shear  vectors  in  Fig.  4.10b  also  differ  by 
about  0.4  nt  s,  which  results  in  Advection  and  Motion  vectors  at  r=4  that  closely  ap¬ 
proximate  those  at  r=  2.  In  contrast  to  Fig.  4.7c.  the  phase  of  interaction  llow  vorticity 
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I;ig.  4.8  (a)-(c)  As  in  Fig.  4.7,  except  using  the  analytical  model  solution  based  on 
C.*  1.9  m  s  and  a  -  130. 
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Fig.  4.9  Vector  differences  between  the  TC  propagation  and  an  averaged  interaction  flow 
velocity  (0-300  km)  at  various  times  (see  legend)  during  a  120  h  integration  of  a 
barotrcpic  numerical  model  (Fiorino  and  Elsberry  1989). 
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Fig.  4.10  As  in  Fig.  4.7,  except  using  the  analytical  model  solution  based  on  C.~  2.4  m/s 
and  a  -  132. 
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is  constant  for  1.5  <;  r  <,  Rr  in  Fig.  4.10c,  which  indicates  that  the  vector  patterns  in  Fig. 
4.10  are  characteristic  of  most  of  the  Self-advection  Region.  Such  consistency  of  the 
vorticity  tendency  balance  is  to  be  expected  for  a  TC  propagating  as  a  quasi-steady  en¬ 
tity,  and  supports  the  dynamical  validity  of  both  the  theoretical  model  and  the  modified 
closure  scheme. 

It  should  be  noted  that  the  phase  of  changes  abruptly  for  r<1.5.  In  addition, 
the  amplitude  of  CA  increases  quite  rapidly  and  actually  becomes  singular  at  the  origin. 
Both  of  these  properties  are  characteristic  of  an  intense  "inner  gyre"  in  the  vorticity  field 
of  the  interaction  flow,  which  must  exist  to  balance  the  singular  nature  of  the 
piecewise-analytic  symmetric  flow  at  the  origin.  Thus,  the  inner  vorticity  gyre  may  be 
purely  an  artifact  of  this  theoretical  model.  By  contrast,  numerical  modeling  studies 
suggest  that  such  inner  vorticity  gyres  result  from  a  dynamical  instability  which  is  sup¬ 
ported  by  the  barotropically  unstable  structure  of  typical  TC  tangential  wind  profiles 
( e.g .,  Willoughby  1988;  Peng  and  Williams  1989).  Similar  gyre  patterns  have  been  sug¬ 
gested  by  aircraft  observations  of  the  inner  wind  fields  of  a  TC,  although  these  may  be 
to  some  extent  due  to  mislocating  the  TC  center.  In  this  model,  the  inner  gyres  have 
no  significant  impact  on  the  propagation  prediction  capability  of  the  internal  closure 
scheme  just  demonstrated. 

The  usefulness  of  the  internal  closure  scheme  depends  on  the  "closure"  vector 
pattern  (i.c.,  Fig.  4.10a)  occurring  at  a  well-defined  and  unique  location  in  the  C.  and  a 
parameter  space.  By  showing  that  the  vorticity  tendency  vector  pattern  changes  dis¬ 
tinctly  in  response  to  small  adjustments  of  C.  and  a  away  from  the  "closure”  values, 
Figs.  4.11-4.14  verify  that  the  closure  point  may  be  quite  precisely  located.  Changes  in 
a  to  137°  (Fig.  4.11)  and  127°  (Fig.  4.12)  induce  a  meridional  component  into  the  Shear 

A 

vector.  In  contrast,  changes  in  C.  to  2.1  m's  (Fig.  4.13)  and  to  2.7  m  s  (Fig.  4.14)  tend 
to  alter  the  relative  magnitudes  of  the  Beta  and  Shear  vectors,  although  the  processes 
are  coupled  to  a  moderate  extent.  The  phase  uniformity  of  (,A  in  the  Self-advection  Re¬ 
gion  is  rapidly  altered  (Fig.  4.11c  and  4.12c)  when  the  specified  propagation  direction 
differs  by  only  +5°  from  the  correct  value.  Although  all  closures  in  this  research  are 
obtained  via  interactive  computer  graphics,  the  well-defined  nature  of  the  closure  vector 
pattern  should  be  amenable  to  computer  automation.  A  formal  proof  that  only  one 
closure  pattern  exists  for  a  particular  set  of  TC  and  environment  parameters  has  not 
been  attempted.  However,  a  careful  search  over  a  wide  range  of  imposed  propagation 
velocities  has  never  revealed  more  than  one  closure  point  in  C.  and  a  parameter  space. 
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Fig.  4.12  As  in  Fig.  4.10,  except  using  C.=  2.4  m/s  and  a=  127 
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Fig.  4.14  As  in  Fig.  4.10,  except  using  C.  =  2.7  m/s  and  a=  132 


The  closure  scheme  quite  precisely  locates  a  "predicted"  direction  of  propagation 
that  is  12°  to  the  left  of  the  actual  value  found  by  FE.  The  cause  for  most  of  this  bias 
is  that  the  phase  of  \j/A  converges  to  a  value  of  9°  at  large  radii  (Fig.  4.15),  rather  than 
to  zero,  as  would  occur  in  a  quiescent  environment  if  the  Dispersion  Region  equation 
was  solved  over  the  entire  model  domain.  The  mathematical  basis  for  this  property  of 
the  solution  is  that  in  (4.17a)  the  terms  associated  with  the  inner  boundary  condition  and 
to  first  order  the  forcing  term  G2/2r  both  decay  proportional  to  r[  for  r>Rr.  As  a  result, 
the  symmetric  angular  windshear  that  produces  the  barotropically  stable  downshear 
phase  shift  in  and  near  the  Self-advection  Region  (Fig.  4.15;  r<8)  still  has  a  small  but 
measureable  influence  in  the  Dispersion  Region  at  large  radii.  Since  the  phase  shift  of 
\pA  at  large  radii  may  always  be  computed  from  the  model  solution,  the  closure-predicted 
propagation  direction  could  be  corrected  by  this  angle  to  increase  the  accuracy  of  the 
model  for  applications  in  which  a  precise  direction  prediction  would  be  of  particular 
importance  (e.g.,  a  propagation  +  advection  track  forecasting  aid).  Such  directional 
corrections  will  not  be  used  here.  As  will  be  shown  in  Chapter  V,  the  corrections  are 
not  required  when  using  this  model  to  initialize  a  numerical  model. 

The  accuracy  of  the  internal  closure  scheme  developed  above  may  be  put  into 
perspective  by  computing  the  24  h  forecast  errors  that  would  result  solely  from  speed 
and  direction  errors  of  0.25  m  s  and  12°.  which  represent  the  difference  between  the 
propagation  predicted  here  via  the  internal  closure  scheme  and  that  observed  by  FE. 
The  speed  error  translates  into  3:25  km  of  forecast  error.  Assuming  a  typical  TC  prop¬ 
agation  velocity  of  2.5  m  s,  the  direction  error  can  cause  2:45  km  of  forecast  error. 
These  errors  are  about  one  quarter  of  the  typical  24  h  official  forecast  error  (e.g., 
Thompson  ei  al.  1981),  and  thus  would  be  of  secondary  importance  compared  to  errors 
induced  by  poor,  inadequate  observations. 

3.  Sensitivity  testing 

a.  Transition  radius  adjustments 

This  model  contains  two  potentially  adjustable  features:  the  location  of  the 
transition  radius  as  determined  by  <5  in  (4.19),  and  the  requirement  for  a  relatively  small 
change  in  ICJ  at  r  =  Rr,  which  may  be  defined  in  terms  of  the  parameter  y  by 


y  s  lim 


!  C,i<'+)  I 


(4.26) 
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Fig.  4.15  Radial  profile  of  streamfunction  amplitude  (solid)  and  phase  (dashed)  corre 
sponding  to  vorticity  profile  shown  in  Fig.  4.10c.  The  amplitude  curve  has  been  nor 
malized  by  the  maximum  amplitude  of  the  streamfunction. 


Thus,  it  is  important  to  test  the  sensitivity  of  the  closure  scheme  to  moderate  adjust¬ 
ments  of  these  conditions,  and  also  to  determine  whether  any  physically  unreasonable 
changes  occur  in  the  predicted  wavenumber  one  gyre  structure. 

Table  4.1  is  a  summary  of  the  sensitivity  tests  performed,  including  the  re¬ 
sultant  changes  in  l,dr.  Test  1  and  (2)  shows  that  decreasing  (increasing)  Rr  results  in  a 
less  (more)  westward  direction  of  propagation,  and  a  slower  (faster)  speed  of  propa¬ 
gation.  In  contrast,  decreasing  (increasing)  y  in  Test  3  (4)  results  in  a  less  (more) 
northward  direction  of  propagation,  and  a  slower  (faster)  speed  of  propagation.  Since 
changes  in  y  are  accomplished  by  adjusting  laij.  in  (4.17b),  Table  4.1  includes  the  values 
of  T*,,.  to  show  that  changes  in  Ta<<,„  are  approximately  proportional  to  the  changes  in 
C.,  which  is  consistent  with  the  idea  that  a  stronger  outer  gyre  structure  will  contribute 
to  stronger  interaction  flow  across  the  center  of  the  TC.  In  each  test,  the  phase  of  \J/A 
at  large  radii  was  found  to  be  similar  in  magnitude  to  that  noted  above  (e.g.,  Test  3, 
11.1°;  Test  4,  7.5°).  In  Test  4,  the  predicted  propagation  speed  and  direction  are  quite 
close  to  FE's  2.65  m/s  at  330°  (a=  120°).  The  wavenumber  one  gyre  structure  associated 
with  the  C.  and  a  of  Test  4  is  shown  in  Fig.  4.16.  If  the  wavenumber  two  component 


Table  4.1  Response  of  the  model-predicted  propagation  velocity  (columns  5,6)  to  four 
combinations  of  the  parameters  S  (column  2)  and  y  (column  3)  as  defined  by  (4.19)  and 
(4.26)  respectively. 
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Fig.  4.16  As  in  Fig.  4.5,  except  using  the  parameter  specifications  of  Test  4  of  Table  4.1. 
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is  ignored,  then  the  numerical  result  of  FE  (Fig.  4.4)  is  quite  closely  approximated  by  the 
theoretical  result  in  Fig.  4.16.  Thus,  the  sensitivity  tests  demonstrate  that  the  theoretical 
model  depends  of  the  adjustable  parameters  in  a  well-defined  and  physically  reasonable 
manner. 

b.  Piecewise-analytic  symmetric  TC  specification 

The  piecewise-analytic  TC  winds  profile  (Fig.  4.2a)  used  above  was  con¬ 
structed  so  that  the  associated  vorticity  gradient  (Fig.  4.2b)  is  continuous  and  equal  to 
the  analytic  vorticity  gradient  at  the  transition  radius.  In  contrast,  the  piecewise-analytic 
vorticity  gradient  is  discontinuous  at  the  interface  and  only  approximates  the  analytic 
vorticity  gradient  in  an  average  sense.  The  motivation  for  approximating  the  analytic 
vorticity  profile  in  this  manner  may  be  illustrated  by  determining  the  model-predicted 
propagation  for  the  piecewise-analytic  profiles  in  Fig.  4.17.  Note  that  the  piecewise- 
analytic  wind  profile  approximates  the  analytic  wind  profile  more  closely  than  in  Fig. 
4.2a.  In  contrast,  the  piecewise-analytic  vorticity  gradient  is  now  discontinuous  and  only 
poorly  approximates  the  analytic  profile  at  both  Rj  and  Rr.  The  model-predicted  prop- 

A 

agation  using  the  parameters  of  the  figure  inset  is  C.=  1.75  and  a=  144°,  which  deviates 
significantly  from  FE's  results  of  2.65  m  s  and  330°  (a  =  120). 

The  cause  of  this  important  error  in  model-predicted  propagation  may  be 
determined  by  eliminating  the  processes  that  cannot  be  responsible.  Since  the 
piecewise-analytic  wind  profiles  in  Figs.  4.2a  and  4.17a  are  nearly  identical,  the  change 
in  predicted  propagation  cannot  be  attributed  to  the  Beta  or  Shear  terms  of  (4.25).  That 
is,  the  asymmetric  forcing  associated  with  a  planetary  vorticity  gradient  and  the  sym¬ 
metric  angular  windshear-induced  stability  of  the  TC  remain  essentially  unchanged.  The 
change  in  the  piecewise-analytic  vorticity  gradient  in  the  Dispersion  Region  cannot  be 
responsible  since  self-advection  processes  are  excluded  there.  The  rough  agreement  be¬ 
tween  the  piecewise-analytic  and  analytic  vorticity  profiles  inside  R,  are  also  quite  similar 
in  Figs.  4.2b  and  4.17b,  and  FE  have  shown  conclusively  that  changes  in  the  symmetric 
TC  inside  300  km  have  a  negligible  effect  on  propagation.  Thus,  it  must  be  concluded 
that  the  significant  reduction  in  the  average  magnitude  of  the  piecewise-analytic  vorticity 
gradient  in  the  outer  annulus  of  the  Self-advection  Region  must  be  responsible  for  the 
change  in  predicted  propagation.  Since  the  piecewise-analytic  vorticity  gradient  roughly 
brackets  the  analytic  curve  in  Fig.  4.17b  near  R„  the  underestimation  of  the  analytic 
vorticity  gradient  just  inside  the  transition  radius  is  inferred  to  be  the  source  of  the  error 
in  predicted  propagation. 
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Fig.  4.17  (a)-(b)  As  in  Fig.  4.2,  except  that  the  parameters  (inset)  of  the  piecewise- 
analytic  TC  wind  structure  have  been  changed  so  that  the  associated  vorticity  gradient 
is  no  longer  continuous  and  equal  to  the  analytic  value  at  the  transition  radius  (r  =  4.75). 


This  analysis  suggests  that  the  sensitivity  of  TC  propagation  on  wind 
strength  in  a  300-800  km  "critical  annulus"  demonstrated  by  FE  is  actually  a  manifesta¬ 
tion  of  the  sensitivity  of  the  self-advcction  process  to  the  radial  gradient  of  symmetric 
TC  vorticity  via  the  Advection  term  of  (4.25).  Thus,  barotropic  TC  propagation  actually 
depends  on  the  second  derivative  of  the  symmetric  windfield  in  the  critical  annulus, 
which  explains  why  seemingly  modest  changes  in  outer  wind  strength  as  used  by  FE  can 
cause  significant  changes  in  propagation.  The  piecewise-defmed  nature  of  the  present 
theoretical  model  tends  to  concentrate  this  dependence  near  the  transition  radius,  and 
also  exaggerate  the  dependence  by  enabling  very  similar  piecewise-analytic  wind  profiles 
to  have  associated  vorticity  profiles  that  differ  more  severely  than  similar  analytic  pro¬ 
files. 


D.  MODEL  RESULTS  PART  II:  INTERNAL  CLOSURE 

The  availability  of  an  internal  closure  scheme  permits  this  model  to  predict  TC 
propagation  for  various  TC  wind  profiles  and  environmental  influences.  Efforts  are 
presently  underway  to  compute  TC  propagation  relative  to  horizontally  variable  envi¬ 
ronmental  winds  from  numerical  models  (personal  communication,  R.  T.  Williams 
1989).  However,  no  published  results  are  available  against  which  the  above  internal 
closure  may  be  compared.  Thus,  it  will  be  assumed  that  the  closure  scheme  with  6  =  1 
and  y  =  1  will  give  results  of  comparable  accuracy  when  applied  with  spatially  variable 
environmental  winds.  An  internal  check  of  the  likely  accuracy  of  the  prediction  can  be 
made  by  comparing  the  closure  vector  plot  for  dynamical  consistency  with  the  associated 
gyre  structure  predicted  by  the  model  for  various  environmental  windficlds. 

1.  Influence  of  TC  structure  change 

Chan  and  Williams  (1987;  CW)  also  used  (4.22)  to  specify  the  initial  TC 
windfield,  and  varied  TC  intensity  and  strength  simultaineously  by  adjusting  the  pa¬ 
rameter  VM.  with  b  =  1.  Three  of  their  parameter  specifications  (Table  4.2;  columns  2-4) 
will  be  used  here.  In  light  of  the  closure  testing  above,  the  piecewise-analytic  profiles 
(Figs.  4.18-4.20)  have  been  constructed  so  that  the  analytical  vorticity  gradients  are  ac¬ 
curately  approximated  near  the  transition  radius.  Note  also  the  increase  in  the  transition 
radius  (Table  4.2;  column  5)  in  response  to  larger  values  of  Fw..  The  comparison  of  the 
theoretically-predicted  propagation  speeds  and  directions  with  those  reported  by  CW 
(Table  2;  columns  6.7)  reveals  slow  speed  biases  and  left  direction  biases  quite  similar  to 
those  in  Section  C  above,  which  indicates  that  the  model  results  are  quite  consistent  with 
regard  to  changes  in  1C  structure.  As  shown  in  subsection  3. a  above,  the  persistent 
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Fig.  4.18  As  in  Fig.  4.2,  except  using  Chan  and  Williams  (1987)  analytic  TC  wind  profile 
parameters  of  VM.-20  m/s,  RM*=  100  km  and  b=  1.0. 
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Fig.  4.19  As  in  Fig.  4.17,  except  for  VM.  =  40  m/s. 
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Fig.  4.20  As  in  Fig.  4.17,  except  for  FM.  =  60  m/s. 


speed  and  direction  biases  could  be  reduced  by  moderately  adjusting  the  parameters  3 
and  y. 

It  should  be  noted  that  the  theoretically-predicted  propagation  directions  in 
Table  4.2  are  essentially  independent  of  changes  in  TC  structure,  whereas  FE  noted  a 
more  westward  propagation  direction  in  response  to  increases  in  TC  outer  wind  strength. 
The  track  of  FE's  "large-weak  vortex"  matched  closely  that  of  the  "basic  vortex"  for  the 
first  24  h,  but  became  increasingly  westward  as  the  integration  proceeded.  Thus,  it  is 
possible  that  the  numerical  model  includes  some  non-steady  propagation  process  that  is 
excluded  from  this  steady-state  model.  An  alternate  explanation  is  that  cyclic  boundary 
conditions  in  the  numerical  model  may  have  a  significant  impact  on  the  propagation 
vector  associated  with  a  large  and  purely  cyclonic  vortex  that  radiates  substantial 
amounts  of  Rossbv  wave  energy  (cf.  Shapiro  and  Ooyama  1989).  Althougth  domain  size 
tests  were  conducted  by  FE  with  the  basic  vortex,  it  does  not  appear  that  such  tests  were 
conducted  for  the  large-weak  vortex.  The  combined  evidence  of  the  theoretical  results 
here  and  the  numerical  results  of  CW  suggest  that  changes  in  TC  strength  result  prima¬ 
rily  in  changes  in  propagation  speed,  but  not  propagation  direction. 


Table  4.2  TC  propagation  velocity  (columns  6,7)  predicted  by  the  theoretical  model  and 
the  numerical  model  (in  parentheses)  of  Chan  and  Williams  (1987)  for  three  values  of 
maximum  symmetric  wind  (column  2).  The  analytic  and  piecewise-analytic  curves  used 
by  the  numerical  and  theoretical  models  respectively  are  shown  in  Figs.  4.18-4.20. 
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2.  Influence  of  uniform  environmental  vorticity 

The  next  step  in  complexity  from  ^-induced  TC  propagation  in  a  quiescent  en¬ 
vironment  is  to  include  a  uniform  environmental  vorticity  field,  i.e.,  ZE^0,  =  0.  Ex¬ 

amples  of  this  situation  would  be:  i)  approximating  the  trade  winds  as  a  nondivergent 
zonal  flow  with  linear  meridional  shear;  and  ii)  approximating  the  subtropical  ridge 
(monsoon  trough)  region  as  a  large-scale  anticyclone  (cyclone)  that  is  in  solid  body  ro¬ 
tation  (cf.  Willoughby  1988). 

Table  4.3  is  a  summary  of  the  closure  results  for  four  different  values  of  envi¬ 
ronmental  relative  vorticity  using  the  same  piecewise-analytic  TC  as  in  Case  2  of  Table 
4.2.  Recalling  that  closure  occurs  at  C.=  2.65  m  s,  a=  132°  for  Case  2  of  Table  4.2,  a 
cyclonic  (anticyclonic)  environmental  vorticity  induces  a  counterclockwise  (clockwise) 
adjustment  in  the  direction  of  TC  propagation  that  is  proportional  to  the  magnitude  of 
Z£..  Note  that  the  theoretical  model  also  predicts  a  greater  TC  propagation  speed  as  ZP 
progresses  from  anticyclonic  to  cyclonic,  which  suggests  that  TC  propagation  may  be  a 
more  dominant  component  of  total  TC  motion  in  regions  of  large-scale  cyclonic  vorticity 
such  as  the  monsoon  trough.  Although  no  numerical  model  confirmation  exists  for  this 


Table  4.3  Theoretically-predicted  propagation  velocity  (columns  3,4)  for  four  values  of 
uniform  environmental  vorticity  (column  2). 
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aspect  of  the  prediction,  it  is  consistent  with  the  faster  propagation  of  westward  moving 
TC's  in  composite  data  (Fig.  1.1b).  The  wavenumber  one  gyre  patterns  associated  with 
Cases  1  and  4  of  Table  4.3  are  shown  in  Fig.  4.21a  and  b  respectively.  Chan  and 
Williams  (1989)  have  noted  gyre  orientation  and  prppagation  changes  in  an  equivalent 
numerical  simulation  that  are  consistent  with  the  results  given  here. 

It  is  important  to  note  that  this  model  cannot  be  internally  closed  if  Zf.  is  much 
less  or  greater  than  the  range  of  values  in  Table  4.3.  The  mathematical  explanation  for 
this  problem  is  that  the  solution  for  £A  in  the  Dispersion  Region  depends  inversely  on 
Z£..  For  sufficiently  large  values  of  environmental  vorticity,  it  becomes  impossible  to 
iterate  to  a  value  of  lcdr  such  that  the  magnitude  of  is  continuous  at  the  transition 
radius.  From  a  dynamical  perspective,  the  assumptions  V£<VS  and  become  un¬ 
acceptably  inaccurate  in  the  outer  part  of  the  Self-advection  Region.  As  a  result,  the 
present  form  of  the  theoretical  model  can  accomodate  only  moderate  values  of  large- 
scale  vorticity,  and  should  not  be  applied  to  situations  involving  interaction  of  the  TC 
with  an  intense  cyclonic  system  such  as  another  TC.  Such  situations  might  be  addressed 
by  extending  the  model  to  include  some  of  the  higher-order  wavenumber  one  processes 
that  were  excluded  here  (Chapter  II). 

3.  Influence  of  environmental  vorticity  gradients 

To  determine  the  model  response  to  gradients  of  absolute  vorticity,  combina¬ 
tions  of  Z£.,  /?,«.  and  <f>  are  devised  that  approximate  the  horizontal  variability  of  the 
typical  large-scale  environments  through  which  TC's  move.  Figs.  4.22a  and  b  illustrate 
parameter  combinations  that  correspond  to  the  regions  south  and  north  of  the  NW 
Pacific  subtropical  ridge  respectively.  As  indicated,  Z£.  is  assumed  to  be  moderately 
anticyclonic  in  both  regions.  The  westward  component  of  the  environmental  vorticity 
gradient  reflects  the  influence  of  the  intense  Asian  heat  low  throughout  much  of  the 
troposphere  during  the  NW  Pacific  typhoon  season.  Similarly,  the  southward  and 
northward  components  reflect  the  influence  of  the  equatorial  trough  and  midlatitude 
baroclinic  regions  respectively.  Such  parameter  combinations  produce  an  absolute 
vorticity  gradient  that  is  larger  and  is  directed  more  poleward  north  of  the  subtropical 
ridge  compared  to  south  of  the  subtropical  ridge  (Table  4.4;  columns  3,4).  Although  the 
analytic  TC  wind  profile  of  Table  4,3  Case  2  is  used  for  both  cases  here,  an  additional 
piecewise-analytic  profile  (Fig.  4.23)  has  to  be  constructed  to  account  for  the  increased 
magnitude  of  J3t/r  north  of  the  subtropical  ridge. 

If  the  phase  shift  of  the  environmental  vorticity  gradient  is  taken  into  account, 
then  the  tendency  vector  patterns  (Figs.  4.24a.b  and  4.25a, b)  from  which  the  propa- 
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Fig.  <1.21  Theoretically-predicted  wavenumber  one  streamfunction  fields  for  (a)  Case  1 
and  (b)  Case  4  of  Table  4.3.  Contour  interval  is  2  x  10 
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Fig.  4.22  Idealized  planetary  and  environmental  vorticity  gradients  for  TC  positions  in 
an  anticyclonic  vorticity  region  (a)  south  and  (b)  north  of  the  subtropical  ridge  in  the 
lower  to  middle  troposphere  of  the  western  North  Pacific  TC  region.  The  units  of  the 
vorticity  gradients  are  lCH'/ir's-',  and  the  units  of  Z£.  are  10"5s*'. 
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Fig.  4.23  As  in  Fig.  4.18,  except  that  the  piecewise-defined  wind  parameters  (see  inset) 
have  been  recalculated  to  account  for  p,fr=  3.5x  10-"  /rr's-1  in  Fig.  4.21b. 


gation  velocities  in  Table  4.4  were  obtained  are  consistent  with  those  shown  earlier  (e.g., 
Fig.  4.10).  In  each  case,  the  Shear  vector  phase  is  opposite  that  of  the  Beta  vector, 
which  is  now  oriented  90°  to  the  left  of  the  direction  of  the  environmental  absolute 
vorticity  gradient.  The  characteristic  patterns  in  the  radial  profiles  of  interaction  flow 
vorticity  magnitude  and  phase  are  also  evident  (Figs.  4.24c  and  4.25c).  The  orientation 
and  amplitude  of  the  wavenumber  one  gyre  streamfunction  are  also  dynamically  con¬ 
sistent  (Fig.  4.26),  which  suggests  that  the  accuracy  of  the  model  with  an  environmental 
vorticity  gradient  included  should  be  similar  to  the  results  shown  earlier.  Although 
equivalent  numerical  results  are  not  available,  the  propagation  directions  predicted  here 
are  quite  similar  to  the  observed  propagation  directions  for  westward  (pre-recurvature) 
and  eastward  (post-recurvature)  moving  TC's  (Fig.  1.1b). 

Table  4.4  Theoretically-predicted  propagation  velocities  (columns  5.6)  in  response  to 
environmental  parameter  combinations  (columns  2-4)  representing  TC  locations  south 
(Case  1)  and  north  (Case  2)  of  the  subtropical  ridge  during  the  western  North  Pacific 
typhoon  season. 


Case 

Zf. 

(X  10-*  5-') 

Pr't' 

(x  10“"  mls~{) 

4> 

(deg) 

C. 

(ms) 

a 

(deg) 

1 

-0.5 

2.0 

150 

2.3 

183 

2 

-0.5 

3.5 

120 

2.9 

156 

E.  SUMMARY 

An  analytical  NDBT  model  based  on  the  principle  of  nonlinear  self-ad vecti on  has 
been  developed  to  predict  the  steady-state  TC  propagation  and  associated  wavenumber 
one  interaction  flow  induced  by  planetary  and  environmental  forcing.  The  model  has 
been  made  analytically  tractable  by  dividing  the  highly  complex  and  generally  intractable 
TC  propagation  process  into  a  number  of  manageable  pieces  via  a  number  of  reasonable 
assumptions.  A  pieccwisc-analytic  solution  for  the  wavenumber  one  interaction  flow 
and  the  associated  propagation  velocity  was  then  constructed  from  analytic  solutions  to 
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Fig.  4.24  (a)-(b)  Vorticity  tendency  vector  diagrams,  and  (c)  vorticity  amplitude  and 
phase  profiles  as  in  Fig.  4.7,  except  for  the  model  solution  for  Case  1  of  Table  4.4. 
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the  individual  pieces  of  the  problem.  The  two  key  aspects  of  this  piecewise-analytic 
modeling  approach  are  to: 

•  divide  the  model  domain  relative  to  a  transition  radius  within  which  is  a  Self- 
advection  Region  where  mutual  advections  by  the  symmetric  TC  and  the  asym¬ 
metric  interaction  flow  are  important,  and  outside  of  which  is  a  Dispersion  Region 
where  such  advections  are  considered  unimportant;  and 

•  approximate  the  symmetric  TC  windfield  by  a  piecewise-analytic  modified  Rankine 
profile  that  closely  matches  the  analytic  TC  wind  profile,  except  near  the  radius  of 
maximum  winds,  and  closely  matches  the  analytic  vorticity  gradient  in  the  Self- 
advection  region  near  the  transition  radius. 

This  theoretical  model  is  distinguished  from  previous  efforts  by  the  capabilities  to: 

•  accurately  predict  both  the  zonal  and  meridional  components  of  TC  propagation 
by  integrating  linear  and  nonlinear  mechanisms  into  a  single  self-advection  process; 

•  accurately  predict  the  wavenumber  one  gyre  structure  responsible  for  TC  propa¬ 
gation  based  on  either  a  knowledge  of  the  propagation  velocity,  or  determining  the 
propagation  from  the  model  dynamics  via  a  closure  scheme; 

•  include  the  influence  of  changes  in  the  symmetric  wind  of  the  TC; 

•  include  the  first-order  effects  of  large-scale  relative  vorticity  gradients  of  arbitrary 
magnitude  and  direction;  and 

•  include  the  first-order  effects  of  moderate  values  of  uniform  relative  vorticity  of  the 
large-scale  environment. 

The  close  agreement  between  the  TC  propagation  velocities  predicted  by  this  model 
and  the  equivalent  numerical  solutions  verifies  that  the  piecewise-analytic  construction 
technique  is  dynamically  sound.  The  propagation  vector  errors  from  this  model  have 
small,  well-defined  biases  that  depend  on  the  two  adjustable  parameters  of  the  model  in 
a  predictable  and  dynamically  reasonable  manner.  This  property,  combined  with  the 
capability  of  the  model  to  predict  TC  propagation  and  gyre  structure  for  many  realistic 
combinations  of  TC  structure  and  environmental  windfields,  suggests  that  this  model  has 
significant  potential  for  use  as  either:  i)  an  initialization  scheme  for  barotropic  numerical 
forecast  models  such  as  SANBAR;  ii)  or  as  part  of  a  "propagation  +  advection"  track 
prediction  aid.  A  demonstration  of  the  potential  usefulness  of  this  model  as  an  initial¬ 
ization  tool  is  provided  in  Chapter  V. 

The  success  of  this  theoretical  model  in  reproducing  the  numerical  results  of  CW 
and  FE  provides  strong  evidence  that  nonlinear  self-advection.  rather  than  linear  prop¬ 
agation.  is  the  barotropic  mechanism  that  contributes  to  the  observed  propagation  of 
TC's.  It  is  important  to  emphasize  that  the  individual  effects  of  /?,  environmental  shear 
and  environmental  vorticity  gradients  are  integrated  within  this  self-advection  model, 
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since  those  processes  together  determine  the  phase  and  amplitude  of  a  wavenumber  one 
interaction  flowfield  that  advects  the  TC  relative  to  the  large-scale  environment.  By 
providing  a  unified  theory  for  barotropic  TC  propagation,  this  model  represents  an  im¬ 
portant  step  toward  the  development  of  a  general  theory  of  TC  motion.  For  the  sake 
of  brevity  and  to  emphasize  the  unifying  aspect  of  self-advection  theory,  this  model  will 
hereafter  be  referred  to  as  the  Barotropic  Self-advection  Model  (BSAM). 


105 


V.  BAROTROPIC  SELF-ADVECTION  MODEL  APPLICATION 


In  this  chapter,  the  usefulness  and  versatility  of  the  BSAM  will  be  demonstrated  in 
three  important  areas.  First,  the  issue  of  how  accurately  the  outer  wind  strength  must 
be  measured  to  adequately  account  for  propagation  of  TC's  will  be  addressed  by  using 
tangential  wind  profiles  based  on  composited  data  as  input.  An  alternate  approach  of 
using  the  BSAM  to  predict  an  "effective"  outer  wind  strength  based  on  past  TC  propa¬ 
gation  will  also  be  outlined.  Second,  a  quantitative  assessment  of  the  extent  to  which 
the  asymmetric  interaction  flow  will  be  accounted  for  in  steering  flow  calculations  will 
be  made  based  on  the  wavenumber  gyre  structure  predicted  by  the  BSAM  for  the  com¬ 
posite  data.  Suggestions  on  how  to  isolate  the  environment  and  interaction  flow  com¬ 
ponents  will  also  be  included.  Finally,  the  feasibility  of  using  the  BSAM  to  initialize 
barotropic  numerical  forecast  models  will  be  demonstrated. 

A.  BSAM  PREDICTIONS  USING  COMPOSITE  DATA 
I.  Preliminary  analysis 

In  Chapter  IV,  the  BSAM  propagation  predictions  were  based  on  piecewise- 
analytic  TC  wind  profiles  that  closely  approximated  the  exponential  wind  profile  (4.22) 
used  in  recent  numerical  studies  (CW  and  FE).  A  similar  exponential  profile  was  also 
used  by  DeMaria  (19S5;  DM).  Exponential  TC  wind  profiles  were  used  in  the  BSAM 
to  verify  the  accuracy  of  the  theoretical  predictions  relative  to  equivalent  numerical  sol¬ 
utions.  In  light  of  the  importance  of  TC  structure  on  barotropic  propagation,  it  is  im¬ 
portant  to  examine  how  well  exponential  profiles  approximate  TC  winds,  and  whether 
some  other  functional  form  might  be  better  for  propagation-prediction  purposes. 

As  shown  in  Fig.  5.1a,  FE  selected  a  "basic  vortex"  profile  to  have  a  radius  of 
15  ms  winds  at  300  km  to  agree  with  typical  composite  observations  of  TC  surface 
tangential  winds,  such  as  those  given  by  Merrill  (19S4)  for  large  Atlantic  hurricanes. 
For  comparison,  a  "composite  pressure-averaged  typhoon"  was  calculated  by  taking  a 
950-150  mb  pressure-weighted  average  of  the  western  North  Pacific  composite  TC 
tangential  wind  data  of  Frank  (1977;  his  Fig.  9).  The  basic  vortex  has  greater  intensity 
and  less  outer  wind  strength  than  typical  typhoons  or  large  hurricanes,  as  represented 
by  either  the  surface  or  the  pressure-averaged  tangential  winds  (Fig.  5.1a).  The  com¬ 
posite  pressure-averaged  wind  profile  for  N\V  Pacific  typhoons  is  significantly  weaker 
at  all  radii  than  the  composite  surface  wind  profile  for  large  Atlantic  hurricanes  because 
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Fig.  5.1  (a)  Radial  profiles  of  Fiorino  and  Elsberry's  "basic  vortex"  (dashed),  the  com¬ 
posite  surface  tangential  winds  of  large  Atlantic  hurricanes  (solid  squares;  Merrill  1984) 
and  the  taneential  winds  of  the  composite  pressure-averaged  typhoon  (open  squares; 
Frank  1977)7  Here  Vsr=  35  m/s  and  RM.=  100  km.  (b)  Analytic  approximations  to  the 
radial  profile  of  the  composite  pressure-averaged  typhoon  (open  squares)  as  defined  by 
(4.22)  with  b  =  0.6  (dashed)  and  (5.1)  (solid;  parameters  in  inset).  Here  rJf.=  20  m  s. 


the  tangential  wind  components  decrease  with  elevation  in  TC's.  Using  surface  wind 
profiles  to  initialize  barotropic  models  may  cause  an  overestimation  of  TC  propagation 
associated  with  barotropic  processes.  The  composite  pressure-averaged  typhoon  wind 
profile  will  be  used  here. 

Both  FE  and  DM  adjusted  the  maximum  wind  scale  VM.  and  the  b  parameter 
in  (4.22)  to  better  fit  the  larger  and  weaker  average  typhoon/hurricane.  Although  the 
exponential  profile  can  be  adjusted  to  approximate  the  wind  profile  of  the  composite 
pressure-averaged  typhoon  (Fig.  5.1b),  the  exponential  profile  does  not  represent  the 
curvature  (and  thus  the  vorticity  gradient)  of  the  composite  profile  nearly  as  w’ell  as  the 
analytic  function 

vs(r)  =  D0r~X  +  D ,(r-  1)  +  D2(r  —  l)2  +  D3(r-  if  ,  (5.1) 

which  is  a  modified  Rankine  profile  with  a  three-term  Taylor  series.  Thus,  (5.1)  will  be 
used  as  an  analytic  representation  for  composite  TC  tangential  windfields  below. 

2.  Propagation  speed  versus  composite  TC  strength 

By  adjusting  the  curvature  of  the  piecewise-analytic  wind  profile  via  the  param¬ 
eter  X,  in  (4.6),  the  outer  wind  strength  of  the  composite  pressure-averaged  typhoon 
may  be  underestimated,  closely  approximated  and  overestimated  in  the  300-800  km 
annulus  (Figs.  5. 2-5.4  respectively),  even  though  the  inner  wind  speeds  remain  essentially 
unchanged.  In  each  case,  the  interface  radius  R,  is  located  at  300  km  as  in  FE  so  that 
the  same  strength  change  methodology  is  used.  An  additional  parameter  C3  is  present 
(see  insets)  because  the  functional  form  of  vs  in  the  Dispersion  Region  has  been  modified 
to 

v5(r)  =  A3,--4  +  B3,  -5  +  C3,  _6  ,  (5.2) 

to  require  that  cCs/cr  be  continuous  at  the  transition  radius.  Although  such  a  change 
adds  two  additional  terms  to  (4.17b),  it  has  a  negligible  influence  on  the  BSAM-predicted 
propagation,  and  serves  only  to  avoid  overly  large  jumps  in  the  piecewise-analytic 
vorticity  gradient  at  Rx.  The  magnitude  of  the  analytic  vorticity  gradient  at  Rr  is  about 
four  times  smaller  for  the  composite  profile  (Fig.  5.2b)  than  for  the  exponential  profiles 
used  earlier  (e.g.,  Fig.  4.2b).  Also,  the  change  of  vorticity  gradient  sign  occurs  well  into 
the  Dispersion  Region  for  the  composite  data  due  to  the  significantly  greater  outer  wind 
strength.  The  vorticity  gradient  of  the  piecewise-analytic  profile  results  in  some  under¬ 
estimation  of  the  analytic  vorticity  gradient  at  R7.  However,  it  is  not  certain  whether 
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Fig.  5.2  (a)  Radial  wind  profiles  of  the  composite  pressure-averaged  typhoon  (dashed) 
as  approximated  by  (5.1)  in  Fig.  5.1b  and  a  piecewisc-analytic  profile  (parameters  in  in¬ 
set)  that  overestimates  outer  wind  strength,  (b)  Piecewise-analytic  and  analytic  vorticity 
gradient  profiles  corrcspondine  to  the  winds  in  (a).  In  both  cases,  FM.=  20  m,'s  and 
Rm.=  100  km. 
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Fig.  5.3  As  in  Fig.  5.2,  except  for  a  piecewisc-analytic  wind  profile  that  closely  approx¬ 
imates  the  outer  wind  strength  of  the  composite  pressure-averaged  typhoon  in  the  300 
to  800  km  radius  "critical  annulus." 
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this  will  introduce  a  speed  bias  into  BSAM  propagation  predictions  because  the  coarse 
radial  resolution  and  limited  accuracy  of  the  composite  data  allows  for  some  subjectivity 
in  selecting  the  curvature  and  fit  of  (5.1). 


Table  5.1  TC  propagation  velocities  (columns  5,6)  predicted  by  the  BSAM  for 
piecewise-analytic  wind  profiles  (see  Figs.  5.2a-5.4a)  that  underestimate,  approximate 
and  overestimate  the  outer  wind  strength  of  the  composite  pressure-averaged  typhoon 
(Cases  1-3  respectively).  In  each  case,  a  quiescent  environment  has  been  assumed.  The 
wind  speed  at  the  transition  radius  R7  (column  3)  and  at  550  km  (column  4)  are  also 
shown. 
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The  propagation  speeds  and  directions  predicted  by  the  BSAM  (Table  5.1)  in¬ 
crease  with  increased  TC  outer  wind  strength,  which  is  in  general  agreement  with  the 
results  of  Chapter  IV  and  previous  studies.  Quasi-linear  relationships  exist  between  the 
propagation  speed  (column  5)  and  the  wind  speed  at  the  transition  radius  (column  3), 
and  the  wind  speed  at  550  km  (column  4).  The  first  relationship  has  more  relevance  to 
the  dynamics  of  the  BSAM  since  is.(Rr.)  is  a  rough  measure  of  the  amount  the  linear 
asymmetric  forcing  generated  in  the  Dispersion  Region,  which  acts  to  strengthen  the 
wavenumber  gyres  and  thus  increase  the  propagation  speed.  The  second  relationship  is 
also  dynamically  relevant  because  550  km  would  be  the  center  of  PC's  "critical  annulus" 
between  300  and  800  km.  The  1  m  s  propagation  speed  increase  for  a  2.4  m  s  increase 
in  TC  tangential  winds  at  550  km  represents  a  potential  for  about  85  km  of  forecast  error 
in  24  h.  This  result  suggests  that  TC  outer  wind  strength  must  be  measured  to  within 
±1  ms  in  order  for  numerical  models  to  be  initialized  with  sufficient  accuracy  to  avoid 
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significant  forecast  errors  due  to  misrepresentation  of  TC  strength.  Important  related 
questions  are:  i)  do  numerical  model  heating  and  dissipative  parameterizations  tend  to 
maintain  the  tangential  wind  profile  to  this  accuracy  throughout  the  forecast  integration; 
and  ii)  are  changes  in  TC  strength  that  occur  during  a  forecast  period  accurately  re¬ 
produced  by  numerical  forecast  models? 

This  sensitivity  of  TC  propagation  on  outer  wind  strength  poses  a  significant 
challenge  for  present  observing  system  technologies.  However,  the  method  just  used  to 
demonstrate  the  wind-strength/propagation-speed  relationship  offers  a  potential  alter¬ 
nate  approach  to  determine  an  "effective"  TC  strength.  In  Table  5.1,  the  independent 
adjustment  of  a  single  parameter  X,  (Rr  varies  as  a  result)  causes  related  changes  to  TC 
outer  wind  strength  and  propagation  speed.  It  is  therefore  possible,  in  principle,  to  seek 
a  value  ofX,  such  that  the  propagation  speed  predicted  by  the  BSAM  matches  the  ob¬ 
served  value  at  the  initial  time.  The  associated  wavenumber  one  gyre  structure  predicted 
by  the  BSAM  could  then  be  used  to  initialize  a  barotropic  forecast  model  like  SANBAR, 
or  provide  guidance  in  selecting  an  appropriate  bogus  vortex  for  baroclinic  forecast 
models. 

Implicit  in  such  a  scheme  to  predict  TC  strength  are  the  requirements  that: 
i)  the  large-scale  environment  is  known  sufficiently  well  to  provide  P,rr,  <f>  and  for  the 
BSAM;  ii)  the  track  of  the  TC  is  resolved  well  enough  to  compute  propagation  relative 
to  the  barotropic  steering  of  the  environment;  and  iii)  an  estimate  of  TC  speed  at  about 
300  km  is  available.  Since  it  is  likely  that  many  combinations  of  inner  intensity  and 
outer  wind  strength  can  give  similar  propagation  velocities,  such  a  scheme  might  not 
need  to  reproduce  the  outer  wind  structure  of  the  TC,  but  only  an  inner- 
intensity 'outer- wind-strength  combination  (i.e.,  an  "effective"  TC  strength)  that  will 
produce  the  correct  propagation.  In  addition,  possible  baroclinic  contributions  to  TC 
propagation  might  also  be  parameterized  to  some  extent  by  such  an  approach. 

B.  INTERPRETATION  OF  COMPOSITE  TC  PROPAGATION  VECTORS 

In  Chapter  I,  comparison  of  observations  of  TC  propagation  with  numerically  de¬ 
rived  depictions  of  /?-induced  wavenumber  one  gyre  structure  (FE)  led  to  the  preliminary 
conclusion  that  the  uniform  flow  region  of  the  gyres  would  not  be  accounted  for  in 
steering  flows  calculated  over  an  annulus  of  5-7°  lat.  radius.  However,  the  wavenumber 
one  gyres  for  the  composite  pressure-averaged  typhoon  due  to  /^-forcing  only  (Fig.  5.5) 
are  somewhat  larger  than  the  gyres  associated  with  the  basic  vortex  profile  of  FE  (Fig. 
4.4).  The  well  resolved  gyre  structure  provided  by  the  BSAM  will  be  used  to  calculate 
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Fig.  5.5  Wavenumber  one  gyre  streamfunction  (dashed,  negative)  associated  with  Table 
5.1  Case  2.  Contour  interval  is  2.0  x  10-W/s.  Dotted  circles  have  radii  of  2°,  4°,  6°,  8° 
and  10°  lat.  The  middle  circle  corresponds  to  the  center  of  the  5°-7°  lat.  radius  annulus 
typically  used  to  compute  steering  from  composite  data  (e.g.,  Chan  and  Gray  1982; 
Holland  1984).  The  parameter  RM.~  100  km. 
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the  extent  to  which  TC  self-advection  is  included  in  such  steering  flow  calculations  or  is 
manifested  as  propagation. 

By  averaging  the  velocities  around  various  circles  centered  on  the  BSAM -predicted 
gyres  (Fig.  5.5;  dotted  lines),  a  set  of  velocity  vectors  representing  interaction  flow  con¬ 
tributions  to  environmental  steering  may  be  obtained  (Fig.  5.6a).  Subtracting  such 
vectors  from  the  BSAM-predicted  propagation  velocity  for  the  composite  pressure- 
averaged  typhoon  (Table  5.1;  columns  5,6)  gives  a  set  of  vectors  that  represent  the 
propagation  of  the  storm  relative  to  the  steering  computed  at  a  particular  radius  (Fig. 
5.6b).  For  composite  steering  computed  at  a  2°  lat.  radius,  most  of  the  interaction  flow 
is  accounted  for  in  the  steering,  whereas  the  propagation  vector  is  about  0.4  m/s  to  the 
west.  In  contrast,  at  the  6°  lat.  radius  where  steering  is  typically  computed,  all  of  the 
westward  component  and  about  one  quarter  of  the  northward  component  are  manifested 
as  propagation.  This  result  is  consistent  with  the  structure  of  the  wavenumber  one  gyre 
uniform  flow  region  in  Fig.  5.5,  which  is  oriented  essentially  north-south  at  6°  lat.  radius, 
but  still  has  a  velocity  nearly  as  large  as  the  innermost  portion  of  the  flow.  This  result 
may  then  explain  why  TC  propagation  computed  from  composite  data  over  a  5-7° 
annulus  tends  to  have  a  predominantly  westward  orientation  (Fig.  1.1).  Although  the 
analysis  here  has  used  gyres  associated  with  only  the  influence  of  planetary  vorticitv,  the 
results  would  also  apply  to  the  more  realistic  environmental  situations  of  Chapter  IV 
Section  D.3,  if  appropriate  phase  shifting  of  the  vector  patterns  is  taken  into  account. 

As  discussed  in  Chapter  III,  the  angular  windshear  in  the  symmetric  TC  produces 
a  down-shear  tilting  of  the  inner  portion  of  the  wavenumber  one  gyre  that  barotropically 
stabilizes  the  TC  to  asymmetric  forcing.  It  is  this  down-shear  tilt  that  produces  the 
characteristic  cyclonic  rotation  and  increasing  magnitude  of  TC  propagation  relative  to 
steering  flows  computed  at  increasingly  larger  radii  (Fig.  5.6b).  This  propagation  vector 
pattern  was  evident  in  early  compositing  studies  (George  and  Gray  1976;  their  Figs.  12 
and  13),  and  is  a  persistent  feature  of  the  recent  compositing  results  based  on  much 
larger  data  sets  presented  by  W.  Gray  at  a  recent  Office  of  Naval  Research  Workshop 
on  Tropical  Cyclone  Motion  (Elsberry  1989).  Such  vector  patterns  may  be  interpreted 
as  evidence  that  observed  TC  propagation  is  a  manifestation  of  self-advection  by 
wavenumber  one  gyres  produced  by  asymmetric  forcing.  Although  baroclinic  processes 
may  also  contribute  to  the  total  wavenumber  one  gyre  pattern,  the  general  agreement 
of  observed  propagation  (Figs.  1.1 )  with  the  predictions  of  barotropic  theory  (Table  4.4) 
suggests  that  barotropic  processes  are  the  dominant  influence. 
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Fig.  5.6  (a)  BSAM -predicted  propagation  (C)  for  the  composite  pressure-averaged 
typhoon  and  contributions  to  environmental  steering  by  the  associated  BSAM-gencratcd 
wavenumber  one  gyre  How  evaluated  around  the  2, 4, 6, 8  and  10°  lat.  radius  circles  of  Fig. 
5.5.  (b)  Differences  between  composite  pressure-averaged  typhoon  propagation  (C)  and 
the  wavenumber  one  gyre  contributions  to  steering  of  (a). 
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For  the  BSAM  to  be  used  as  a  numerical  model  initialization  tool  or  as  part  of  a 
propagation  +  advection  forecasting  aid,  the  interaction  flow  must  be  effectively  excluded 
from  the  environmental  windfield.  Although  filtering  is  often  used  in  an  attempt  to  re¬ 
move  the  influence  of  the  TC  from  the  large-scale  environment,  alternate  methods  to 
compute  steering  that  better  exclude  the  uniform  flow  region  of  the  gyres  should  also 
be  investigated.  For  example,  the  steering  might  be  computed  along  circles  as  in  Fig. 
5.5,  but  including  only  the  90°  arcs  centered  on  a  line  normal  to  the  direction  of  the 
large-scale  absolute  vorticity  gradient  (Fig.  5.7).  Using  this  approach  (Fig.  5.8a),  the 
interaction  flow  contribution  to  the  6°  lat.  radius  steering  accounts  for  only  half  the 
northward  component  of  propagation,  and  the  contribution  to  the  8°  lat.  radius  steering 
is  essentially  zero.  As  a  result,  nearly  all  of  the  interaction  flow  influence  is  manifested 
as  propagation  relative  to  steering  computed  at  8°  lat.  radius  or  larger  (Fig.  5.8b).  Thus, 
it  may  be  desirable  to  depart  from  traditional  steering  computation  methods  to  better 
separate  the  large-scale  environment  from  the  influence  of  the  TC. 

C.  NUMERICAL  MODEL  INITIALIZATION 

The  wavenumber  one  gyre  structures  provided  by  the  BSAM  (based  on  either  an 
external  specification  or  an  internal  prediction  of  TC  propagation)  should  be  ideal  for 
initializing  a  nondivergent.  barotropic  (NDBT)  numerical  model  so  that  quasi-steady 
motion  occurs  immediately.  Comparisons  with  the  NDBT  numerical  model  of  CW  will 
be  used  to  provide  a  preliminary  verification  of  this  assertion.  The  model  uses  a 
/?-plane  approximation,  a  4000x4000  km  channel  domain  with  an  east-west  cyclic 
boundary  condition  and  a  resolution  of  20  km. 

1.  Quiescent  environment  predictions 

Using  the  symmetric,  exponential  TC  wind  profile  of  Table  4.2  Case  2.  the  CW 
model  predicts  a  nearly  constant  /^-induced  propagation  direction  of  about  330°  during 
the  first  4S  h  (Fig.  5.9).  In  contrast,  the  propagation  speed  increases  from  zero  initially 
to  about  2.65  m  s  at  48  h.  During  this  period,  the  asymmetric  gyres  develop,  and  then 
the  propagation  speed  becomes  almost  constant  (2.8  m  s)  after  the  gyres  are  established. 
In  a  real  storm,  these  gyres  arc  presumably  present  continually,  and  change  only  in  re¬ 
sponse  to  variations  in  TC  structure  or  environmmental  forcing.  The  "forecast  error" 
associated  with  the  transient  acceleration  period  evident  in  Fig.  5.9  may  be  estimated 
by  comparing  the  U-48  h  track  with  a  48  h  displacement  between  say  36  h  and  84  h 
during  which  the  propagation  speed  is  quasi-steady,  l  hc  comparison  is  facilitated  by 
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Fig.  5.7  As  in  Fig.  5.5,  except  for  including  in  steering  flow  computation  only  the  winds 
in  the  90°  arcs  to  the  right  and  left  of  the  planetary  vorticity  gradient. 
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I  ig.  5.9  /?-induced  propagation  tracks  and  6-hour  positions  predicted  by  the  numerical 
model  of  Chan  and  Williams  (1987)  for  an  initial  TC  wind  profile  corresponding  to  Table 
4.2  Case  2  from  0-4S  h  (solid  circles),  and  for  the  period  36-84  h  (crossed  circles)  except 
translated  so  that  the  36  h  position  corresponds  to  the  initial  position  of  the  TC. 


translating  the  36-84  h  track  so  that  the  36  h  position  of  the  TC  corresponds  to  the  in¬ 
itial  position. 

The  along-track  error  of  about  150  km  suggested  by  this  approach  may  be  an 
underestimation  since  the  TC  has  slowed  to  about  2.6  m/s  by  84  h  due  to  a  weakening 
of  outer  wind  strength  as  the  TC  is  displaced  northward  on  a  /J-plane  (cf.,  FE).  The 
much  smaller  cross-track  error  is  caused  by  the  slightly  more  northward  track  of  the  TC 
from  36-84  h  as  compared  to  0-48  h.  Since  this  error  is  small,  and  may  be  associated 
with  different  TC  strength  during  the  two  48  h  periods  or  boundary  influences,  it  will 
be  ignored. 

The  numerical  model  is  now  modified  to  include  the  wavenumber  one  gyre 
structure  generated  by  the  BSAM  given  the  propagation  %'elocity  of  2.8  m  s  at  330° 
(a  =120°)  reported  by  CW,  or  the  BSAM-predicted  propagation  velocity  (Table  4.2, 
Case  2),  as  shown  in  Figs.  5.10a  and  5.10b  respectively.  In  both  cases,  essentially  steady 
propagation  of  the  model  TC  occurs  immediately  (Fig.  5.11).  The  propagation  tracks 
for  the  two  asymmetric  initializations  agree  remarkably  well,  despite  the  12°  difference 
between  the  initial  direction  of  propagation  predicted  by  the  BSAM  and  the  value  ob¬ 
served  by  CW.  This  indicates  that  the  numerical  model  has  rapidly  adjusted  for  differ¬ 
ences  between  the  analytically-predicted  gyres  and  those  that  would  be  generated  from 
an  initially  symmetric  TC  wind  profile  {e.g.,  Fig.  4.4.).  This  is  consistent  with  the  results 
of  Chapter  III  regarding  short  adjustment  times  near  the  radius  of  maximum  winds, 
since  only  the  central  uniform  flow  region  of  the  BSAM-generated  gyres  is  altered  by 
modest  changes  in  propagation  direction.  Thus,  the  more  important  aspect  of  the  in¬ 
itialization  with  the  BSAM-generated  gyres  included  is  that  the  proper  propagation 
speed  is  established  immediately. 

Using  the  36-84  h  track  in  Fig.  5.11  as  a  benchmark,  the  two  initializations  of 
the  numerical  model  with  ^-induced  gyre  structures  produce  an  along-track  forecast  er¬ 
ror  of  about  60  km.  As  indicated  above,  this  error  is  likely  an  overestimation,  because 
the  36-S4  h  track  includes  an  anomalously  slow  portion  around  84  h.  Thus,  the  along- 
track  error  at  4S  h  has  been  reduced  from  greater  than  150  km  to  a  value  under  60  km. 
Part  of  the  apparent  cross-track  error  of  38  km  relative  to  the  translated  84  h  position 
of  the  TC  is  an  artifact  of  the  more  northward  track  of  the  TC  from  36  h  to  84  h,  as 
mentioned  above.  Comparing  the  0-48  h  track  in  Fig.  5.9  with  the  0-4S  h  tracks  in  Fig. 
5.1 1  indicates  that  the  cross-track  error  is  actually  about  20  km. 
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Fig.  5.10  Wavenumber  one  gyre  streamfunction  patterns  generated  by  the  BSAM  for 
the  symmetric  TC  of  Table  4.2  Case  2  for  (a)  Chan  and  Williams  propagation  velocity 
of  2.8  m  s  at  330°  (a=  120°),  and  (b)  a  BSAM-predictcd  propagation  velocity  of  2.65 
m  s  at  <z  =  132°.  Areal  extent  of  the  figures  corresponds  to  the  domain  of  the  numerical 
model,  and  the  streamlunctions  have  been  linearly  adjusted  to  zero  within  20  gridpoints 
of  the  domain  boundaries. 
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Fig.  5.11  The  /{-induced  propagation  tracks  and  6-hour  positions  predicted  by  the  nu¬ 
merical  model  initialized  with  the  symmetric  TC  as  in  Fig.  5.9,  except  including  the  gyre 
structures  of  Fig.  5.10a  (solid  circles)  and  Fig.  5.10b  (open  circles).  For  comparison,  the 
translated  36-84  h  track  (crossed  circles)  from  Fig.  5.9  is  also  shown. 


2.  Results  in  a  linearly-sheared  environmental  current 

As  shown  in  Table  4.3,  the  BSA\1  predicts  significant  changes  in  TC  propa¬ 
gation  speed  and  direction  for  different  values  of  uniform  environmental  vorticity.  In¬ 
itialization  tests  may  also  be  made  with  the  CW  numerical  model  by  including  a 
linearly- sheared  zonal  current  defined  by 

«£*  =  -ZE.y.,  (5.3) 

which  has  zero  velocity  at  the  initial  position  of  the  TC.  The  numerical  model  is  then 
integrated  for  48  h  using  the  symmetric  TC  vortex  of  Table  4.2  Case  2  with  and  without 
the  BSAM -predicted  gyres  for  Table  4.3  Cases  1  and  2.  Only  BSAM -predicted  propa¬ 
gation  velocities  can  be  used  here  since  numerical  predictions  of  TC  propagation  relative 
to  a  sheared  zonal  current  are  not  presently  available,  although  work  is  proceeding  in 
this  area  (Williams  1989;  personal  communication). 

Comparison  of  the  0-48  h  tracks  of  the  initially  symmetric  and  BSAM -initialized 
TC's  in  the  uniformly  anticyclonic  (Fig.  5.12)  and  cyclonic  (Fig.  5.13)  zonal  currents 
reveals  a  peculiar  aspect  that  warrants  detailed  analysis.  In  the  anticyclonic  shear  case, 
the  distance  between  the  24  h  and  30  h  positions  of  the  BSAM -initialized  TC  is  equiv¬ 
alent  to  a  speed  of  3.3  m/s.  Since  the  TC  is  moving  essentially  due  north  during  this 
period,  such  a  movement  can  only  be  due  to  a  /(-induced  propagation  that  is  substan¬ 
tially  faster  than  the  2.3  m,  s  initially  provided  by  the  BSAM-generated  gyre.  This  result 
may  be  due  to  the  BSAM  having  a  slow  bias  in  the  anticyclonic  shear  case.  However, 
other  evidence  suggests  the  propagation  speed  increase  may  be  associated  with  a  non¬ 
steady  process  present  in  the  numerical  model.  For  example,  the  42  h  to  48  h  movement 
of  the  BSAM-initialized  TC  is  equivalent  to  3.9  m  s,  which  is  significantly  faster  than  the 
24  h  to  30  h  speed  computed  above.  In  additional,  the  42-48  h  speed  of  the  initially 
symmetric  TC  in  Fig.  5.12  is  faster  (3.0  m  s)  than  the  counterpart  in  a  quiescent  envi¬ 
ronment  (Fig.  5.9),  which  was  noted  earlier  to  have  a  42-48  h  speed  of  2.65  m/s.  Thus, 
the  evidence  strongly  suggests  that  the  influence  of  uniform  anticyclonic  environmental 
shear  causes  a  nonstcady  alteration  of  /?- induced  TC  propagation  that  cannot  be  ac¬ 
counted  for  by  a  steady-state  model  such  as  the  BSAM. 

A  similar  result  is  not  found  for  the  cyclonic  shear  case.  The  42-48  h  speed  of 
4.2  m  s  for  the  BSAM-initialized  TC  in  cyclonic  shear  (Fig.  5.13)  also  is  substantially 
faster  than  the  propagation  of  3.0  nt  s  generated  by  the  BSAM  gyres  (Table  4.3;  Case 
2.  column  3).  However,  since  this  track  is  largely  westward,  a  significant  fraction  of  this 
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Fig.  5.12  Tracks  and  6-hour  positions  predicted  by  the  numerical  model  for  /?-induced 
TC  propagation  in  an  zonal  current  with  linear  anticyclonic  shear  of 
ZE.  =  —  5.0  x  10-‘s_1  using  the  symmetric  initial  TC  wind  profile  of  Fig.  5.9  with  (open 
circles)  and  without  (solid  circles)  BSAM-generated  wavenumber  one  gyres. 
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total  motion  is  due  to  the  1.5  m,'s  westward  speed  of  the  environmental  current  at  that 
location.  Subtracting  the  environmental  velocity  vector  results  in  a  //-induced  propa¬ 
gation  velocity  of  2.9  m  s  in  a  direction  as:  144°.  The  difference  between  this  result  and 
initial  propagation  specified  in  Table  4.3  Case  2  is  insignificant  compared  to  the  large 
change  in  propagation  velocity  that  occurs  in  the  anticyclonic  case  above. 

The  solutions  of  Chan  and  Williams  (1989)  for  //-induced  TC  propagation  in  a 
linearly-sheared  zonal  current  also  contain  anomalously  fast  propagation  in  anticyclonic 
shear,  but  no  measurable  changes  in  cyclonic  shear.  Thus,  this  unexplained  behavior  is 
not  limited  only  to  the  BSAM.  The  increased  propagation  speed  in  anticyclonic  shear 
is  certainly  not  intuitively  plausible,  since  the  results  of  FE  (if  applicable)  would  suggest 
that  the  weakening  of  the  outer  wind  strength  of  the  TC  due  to  a  superposition  of  an 
anticyclonically  sheared  current  should  cause  slower  propagation.  Additional  numerical 
integrations  are  being  performed  to  gain  a  better  understanding  of  the  dynamical  or 
perhaps  numerical  basis  for  this  behavior  (personal  communication,  R.  T.  Williams 
1989). 

The  potentially  nonsteady  influence  of  linear  environmental  windshear  on 
//-induced  TC  propagation  make  a  precise  determination  of  how  well  the  BSAM  has  in¬ 
itialized  the  CW  numerical  model  difficult.  However,  several  qualitative  observations 
and  comments  can  be  made.  It  is  evident  in  Figs.  5.12  and  5.13  that  significant  cross¬ 
track  as  well  as  along-track  errors  exist  between  the  positions  of  initially  symmetric  and 
BSAM -initialized  TC's  when  a  spatially  variable  environment  is  present.  As  in  the 
quiescent  environment  case  above,  the  along-track  differences  are  due  to  the  nonzero 
initial  speed  of  the  BSAM -initialized  TC.  The  cross-track  differences  are  in  part  due  to 
a  slightly  different  starting  direction  for  the  BSAM-initialized  TC,  which  in  effect  ac¬ 
counts  for  phase  shifting  of  the  wavenumber  one  gyres  by  the  "past"  influence  of  envi¬ 
ronmental  shear.  Additional  research  will  be  necessary  to  understand  to  what  extent  the 
present  asymmetric  structure  of  TC's  reflects  past  environmental  influences.  The  cross¬ 
track  differences  are  also  due  to  the  nonzero  initial  propagation  speed  of  the 
BSAM-initialized  TC  which,  continually  subjects  the  TC  to  different  environmental 
winds  compared  to  the  initially  symmetric  TC.  Thus,  the  tracks  of  the  BSAM-initialized 
and  initially  symmetric  TC  diverge  with  time.  The  divergence  is  much  greater  in  the 
cyclonic  environment  case  due  to  the  greater  speed  of  propagation  and  the  cyclonic  ro¬ 
tation  of  the  /(-induced  gyres  such  that  a  larger  zonal  motion  component  (propagation 
+  environmental  advcction)  is  produced.  It  should  also  be  emphasized  that  the  4S  h 
along-track  and  cross-track  differences  of  up  to  2S0  km  and  80  km  respectively  in  Figs. 


5.12  and  5.13  indicate  that  initializing  a  numerical  model  with  only  a  symmetric  TC  wind 
profile  has  much  more  serious  consequences  when  spatially  spatially  variable  environ¬ 
mental  winds  are  present  than  for  a  quiescent  environment  (Fig.  5.11). 

D.  SUMMARY 

In  this  chapter,  three  basic  but  important  applications  of  the  BSAM  have  been 
demonstrated.  First,  the  BSAM  has  been  used  with  composite  observations  to  predict 
the  barotropic  propagation  due  to  /?  for  a  composite  pressure-averaged  typhoon,  and  to 
predict  the  change  in  propagation  speed  that  might  be  expected  for  a  typhoon  with 
larger  or  smaller  outer  wind  strength.  The  results  suggest  that  misrepresenting  the  outer 
wind  strength  by  about  2.5  m/s  will  cause  a  1  m/s  error  in  barotropic  propagation,  which 
represents  a  potential  forecast  error  of  85  km  in  24  h.  This  application  also  outlined  a 
method  by  which  the  BSAM  might  be  employed  to  predict  an  ''effective"  outer  wind 
strength  based  on  a  priori  knowledge  of  TC  propagation. 

Second,  analysis  of  the  /^-induced  gyre  structure  predicted  by  the  BSAM  for  the 
composite  pressure-averaged  typhoon  provided  evidence  that  the  TC  propagation  vec¬ 
tors  in  composite  data  (Fig.  1.1)  are  indeed  manifestations  of  wavenumber  one  gyres 
associated  with  self-advection  processes.  An  alternate  method  for  computing  environ¬ 
mental  steering  that  better  separates  the  large-scale  environment  from  the  self-advection 
flow  was  also  demonstrated. 

Finally,  and  most  importantly,  the  potential  usefulness  of  the  BSAM  as  an  initial¬ 
ization  tool  for  barotropic  numerical  forecast  models  was  demonstrated.  In  the  case  of 
a  quiescent  environment,  initializing  the  numerical  model  with  BSAM-predicted 
wavenumber  one  gyres  reduced  the  along-track  forecast  error  at  48  h  from  more  than 
150  km  to  less  than  60  km.  When  a  linearly-sheared  zonal  environmental  current  was 
included,  significantly  larger  along-track  and  cross-track  differences  of  up  to  280  km  and 
80  km  respectively  were  observed  between  the  48  h  positions  of  an  initially  symmetric 
TC  and  a  BSAM-initialized  TC. 
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VI.  CONCLUSION 


The  principal  achievement  of  this  research  has  been  the  development,  testing  and 
preliminary  application  of  a  comprehensive  theoretical  model  for  predicting  tropical 
cyclone  (TC)  propagation  and  associated  asymmetries  based  on  the  barotropic  self- 
advection  process  identified  in  recent  numerical  studies.  The  model  is  comprehensive  in 
the  sense  that  it  includes  the  first-order  effects  of  all  of  the  dynamical  processes  that  are 
presently  understood  to  be  important  to  barotropic  TC  propagation:  gradients  of  plan¬ 
etary  and  environmental  vorticitv,  changes  in  TC  wind  structure,  and  environmental 
windshear.  Since  the  model  is  based  on  a  single  self-advection  principle  through  which 
individual  external  forcing  effects  collectively  determine  the  phase  and  amplitude  of  a 
wavenumber  one  gyre  flow  that  advects  the  TC,  the  nomenclature  Barotropic  Self- 
Advection  Model  (BSAM)  has  been  adopted.  Because  of  the  dynamical  complexity  of 
barotropic  TC  self-advection,  a  rather  lengthy  chain  of  hypotheses,  assumptions  and 
approximations  were  employed  to  make  the  problem  analytically  tractable.  The  princi¬ 
pal  "links"  in  this  development  chain  are  summarized  in  Section  A  and  a  summary'  of 
results  follows  in  Section  B. 

A.  OVERVIEW  OF  PRINCIPAL  ANALYSIS  TECHNIQUES 
1.  Dissection  of  the  problem  and  model  formulation 

Based  on  prior  numerical  findings  that  TC  propagation  relative  to  steering  re¬ 
sults  from  wavenumber  one  asymmetries  induced  by  TC-environment  interaction,  the 
total  flow  field  has  been  partitioned  into  three  components:  i)  a  specified  large-scale  en¬ 
vironment;  ii)  a  symmetric  cyclone  circulation;  iii)  and  an  unknown  asymmetric  inter¬ 
action  flow  component  that  is  presumed  to  be  responsible  for  TC  propagation.  A 
transformation  to  a  reference  frame  moving  with  the  TC  has  been  utilized,  and  the 
symmetric  TC  component  is  assumed  steady  in  such  a  reference  frame  to  further  simplify 
the  problem.  By  subtracting  from  the  total  problem  an  equation  governing  the  large- 
scale  environment  in  the  absence  of  the  TC,  an  equation  for  the  evolution  of  the  inter¬ 
action  flow  in  response  to  external  forcing  has  been  obtained.  Thus,  the  partitioning  and 
transformation  process  significantly  simplified  the  problem  by  focusing  the  analysis  on 
that  portion  of  total  flow  field  that  is  relevant  to  TC  propagation. 

1  lie  mathematical  complexity  associated  with  the  polar  symmetry  of  the  prob¬ 
lem  and  the  extreme  radial  anisotropy  of  the  symmetric  TC  flow  variables  makes  it  very 
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difficult  to  find  a  single  solution  that  is  analytic  in  the  whole  model  domain.  Thus,  it 
has  been  assumed  that  a  piecewise-analytic  solution  obtained  using  a  matched- 
asymptotic  approach  will  be  a  reasonable  approximation.  In  the  preliminary  model  de¬ 
velopment  (Chapter  II),  this  assumption  entailed  dividing  the  model  domain  relative  to 
a  “transition  radius."  Inside  this  radius  is  a  Self-advection  Region  in  which  mutual 
advections  by  the  symmetric  TC  and  interaction  flow  dominate,  and  outside  this  radius 
is  a  Dispersion  Region  in  which  those  advections  are  unimportant  relative  to  the  ex¬ 
ternal  asymmetric  forcing.  After  formally  showing  that  only  wavenumber  one  processes 
can  contribute  significantly  to  TC  propagation,  the  problem  has  been  further  simplified 
by  retaining  only  the  first-order  contributions  to  wavenumber  one  processes  in  the 
interaction  flow  governing  equation. 

To  overcome  the  problem  of  radial  anisotropy  noted  above,  the  model  domain 
has  been  subdivided  further  (Chapter  IV)  by  assuming  that  the  symmetric  TC  winds  in 
the  Self-advection  Region  can  be  reasonably  represented  by  a  set  of  piecewise-analytic, 
modified  Rankine  segments.  A  solution  valid  for  the  whole  model  domain  then  has  been 
assembled  from  solutions  valid  in  individual  annular  regions  by  imposing  matching 
conditions  on  solutions  at  the  annulae  interfaces.  Thus,  the  intractability  of  the  full 
problem  has  been  overcome  by  subdividing  the  problem  sufficiently  to  make  the  indi¬ 
vidual  pieces  analytically  tractable.  The  success  of  this  systematic  dissection  approach 
to  modeling  TC  self-advection  has  been  demonstrated  by  several  comparisons  with 
complete  nu^rical  solutions  and  with  some  observational  results. 

2.  Free  versus  forced  transient  analysis 

An  examination  of  the  response  of  initially  imposed  perturbations  on  a  Rankine 
vortex  provides  important  theoretical  insights  into  why  TC-like  vortices  in  the 
barotropic  numerical  models  propagate  steadily  in  a  slightly  deformed  state  in  response 
to  persistent  environmental  forcing,  rather  than  rapidly  distorting  and  dispersing.  The 
unforced  initial  value  problem  analyzed  in  Chapter  III  represents  the  homogeneous 
counterpart  to  the  forced  steady  propagation  problem  addressed  in  Chapter  IV.  By 
choosing  initial  conditions  with  the  same  spatial  structure  as  in  the  forcing  terms  of  the 
steady-state  problem,  the  forced  transient  adjustments  in  numerical  TC  models  may  in¬ 
terpreted  in  terms  of  the  analytically  tractable  unforced  problem.  This  modeling  tech¬ 
nique  permits  a  formal  mathematical  analysis  of  barotropic  vortex  adjustment  to 
imposed  asymmetries  that  was  not  possible  for  the  related  forced  problem. 
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3.  Closing  the  Barotropic  Self-advection  Model  (BSAM) 

Although  the  transformation  to  a  reference  frame  moving  with  the  TC  was  an 
important  step  in  dividing  the  TC  self-advection  problem  into  manageable  pieces,  the 
presence  of  the  propagation  speed  and  direction  in  the  interaction  flow  equation  re¬ 
presented  two  additional  unknowns.  Thus,  a  closure  scheme  based  on  numerical  results 
has  been  devised  in  which  it  is  hypothesized  that  the  correct  propagation  velocity  will 
be  associated  with  a  particular  wavenumber  one  gyre  structure  that  has  an  approximate 
balance  between  vorticity  tendencies  due  to  the  advection  of  large-scale  absolute 
vorticity  by  the  symmetric  TC  and  the  shearing  of  interaction  flow  vorticity  by  the 
symmetric  TC  winds.  Specifically,  the  advection  of  the  vorticity  of  the  symmetric  TC 
by  the  wavenumber  one  interaction  flow  is  assumed  to  account  for  all  of  the  propagation 
vector,  except  for  a  0.4  m  s  component  associated  with  linear  Rossby  wave  propagation 
of  the  dominant  scale  of  the  TC.  The  closure  scheme  depends  sensitively,  but  predict¬ 
ably,  on  the  two  adjustable  features  of  the  BSAM  associated  with  the  transition  radius. 
In  a  comparison  with  equivalent  numerical  model  solutions,  the  closure  scheme  error 
was  0.25  m  s  and  12°  without  any  adjustment  of  the  two  available  parameters. 

B.  SUMMARY  OF  RESULTS 
1.  Barotropic  vortex  stability 

The  first  important  result  of  this  research  is  that  the  angular  windshear  of  the 
symmetric  circulation  of  TC-like  vortices  acts  to  make  the  vortex  barotropically  stable 
to  small  asymmetric  perturbations  from  axisymmetry.  The  windshear  in  the  symmetric 
circulation  tilts  the  perturbations  downshear,  which  results  in  a  barotropically  stable 
transfer  of  kinetic  energy  from  the  perturbation  to  the  "basic  state"  represented  by  the 
symmetric  vortex.  The  damping  dependence  is  an  algebraic  continuous  spectrum-type 
response  analogous  to  the  damping  of  perturbations  imposed  on  an  /-plane  Couette 
flow.  The  damping  rate  is  proportional  to  perturbation  wavenumber,  and  is  faster 
(slower)  for  perturbations  with  a  radial  structure  that  decays  rapidly  (slowly)  with  in¬ 
creasing  radius.  A  wavenumber  one  perturbation  associated  with  the  advection  of 
planetary  vorticity  by  the  symmetric  TC  loses  about  80°'o  of  the  initial  kinetic  energy 
by  24  h,  which  agrees  well  with  the  adjustment  times  in  recent  numerical  modeling 
studies  of  TC  adjustment  on  a  /?- plane.  Thus,  the  quasi-steady  propagation  of  TC's  in 
barotropic  models  may  be  regarded  as  a  balance  between  the  transfer  of  kinetic  energy 
from  the  symmetric  TC  to  the  wavenumber  one  interaction  flow  due  to  external  asym- 
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metric  forcing  and  the  transfer  of  kinetic  energy  to  the  symmetric  TC  from  the 
wavenumber  one  asymmetry  due  to  the  symmetric  angular  windshear. 

The  presence  of  significant  tangential  windshear  in  TC's,  and  the  qualitative 
similarity  of  TC  propagation  tracks  in  baroclinic  and  barotropic  models,  strongly  sug¬ 
gests  that  this  barotropic  vortex  stability  mechanism  also  may  play  a  fundamental  role 
in  actual  TC  propagation.  It  should  be  emphasized  that  the  asymmetry-damping  influ¬ 
ence  of  symmetric  angular  windshear  is  not  limited  to  vortices  scaled  to  approximate 
TC's.  Rather,  the  shear-damping  mechanism  is  applicable  to  any  vortical  flow  that  may 
be  reasonably  approximated  by  barotropic  dynamics,  and  in  which  the  strength  of  the 
vortical  flow  is  strong  compared  to  other  influences  over  a  significant  horizontal  area. 
Examples  are  Gulf  Stream  rings  and  other  intense  ocean  eddies. 

2.  Dependence  of  propagation  speed  on  TC  strength 

The  piecewise-analytic  symmetric  TC  wind  profile  in  the  BSA.M  permits  propa¬ 
gation  predictions  for  various  tangential  wind  profiles.  During  testing  of  the  closure 
scheme,  it  has  been  demonstrated  that  the  BSAM-predicted  propagation  speed  depends 
strongly  on  the  magnitude  of  the  symmetric  vorticity  gradient  in  the  vicinity  of  the 
transition  radius.  This  result  complements  the  conclusions  of  Fiorino  and  Elsberry 
(1989)  regarding  strength  changes  quite  well,  since  the  transition  radius  is  typically  lo¬ 
cated  near  the  middle  of  their  300-800  km  critical  annulus.  Specifically,  the  BSAM 
demonstrates  that  the  speed  of  TC  propagation  depends  on  the  second  derivative  of  the 
TC  tangential  winds,  which  explains  why  the  propagation  wind-strength  dependence  is 
so  strong. 

The  model  has  also  been  used  with  composite  TC  tangential  wind  profiles  to 
obtain  a  practical  estimate  of  the  wind  measuring  accuracy  required  to  avoid  significant 
forecast  error  due  to  misrepresenting  TC  strength.  Piecewise-analytic  wind  profiles  have 
been  constructed  that  underestimate,  approximate,  or  overestimate  the  outer  wind 
strength  of  a  composite  pressure-averaged  typhoon.  A  quasi-linear  relationship  is  found 
between  the  symmetric  TC  windspeed  at  550  km  and  the  associated  propagation  speed. 
Specifically,  the  speed  of  propagation  increases  by  1  m  s  for  a  2.4  m  s  increase  in  rs.(550 
km),  which  suggests  that  the  TC  outer  wind  strength  must  be  measured  to  about  ±1 
m  s  to  avoid  having  24  h  forecast  errors  greater  than  85  km.  The  capability  to  alter 
outer  wind  strength  by  changing  only  one  parameter  suggests  that  the  BSAM  might  be 
used  to  estimate  an  "effective"  TC  outer  wind  strength  given  sufficiently  accurate  infor¬ 
mation  about  the  environmental  windficld,  the  past  propagation  of  the  TC  relative  to 
environmental  steering  and  the  intensity  of  the  storm. 
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3.  Dependence  of  propagation  on  environmental  vorticity 

The  influence  of  uniform  environmental  vorticity  in  altering  the  phase  of  the 
wavenumber  one  gyres  in  the  Dispersion  Region  can  be  included  in  the  BSA.M  for  values 
of  Z£.  (x  10-5  s_l)  inside  the  range  -0.5  to  1.5.  An  increased  (decreased)  propagation 
speed  and  a  more  (less)  westward  propagation  direction  is  predicted  when  cyclonic 
(anticyclonic)  uniform  vorticity  is  present  compared  to  /^-induced  propagation  in  a 
quiescent  environment.  Although  the  direction  changes  are  generally  consistent  with 
propagation  vectors  in  an  unpublished  numerical  study  by  Chan  and  Williams  (1989),  a 
verification  for  the  BSA.M  -predicted  speed  changes  is  not  available.  An  unresolved  issue 
here  is  that  a  steadily  increasingly  propagation  speed  is  evident  in  the  Chan  and  Williams 
solutions  for  anticyclonic  shear,  whereas  an  essentially  steady  propagation  speed  is 
found  for  cyclonic  shear.  If  the  steadily  evolving  propagation  speed  and  asymmetry 
situation  is  relevant  to  TC's,  then  the  question  of  how  long  the  TC  has  been  in  a  par¬ 
ticular  environment  may  be  important  to  track  forecasting. 

BSAM  predictions  of  TC  propagation  have  been  made  for  situations  involving 
zonal  and  meridional  gradients  of  environmental  vorticity,  the  influence  of  p,  and  a 
uniform  component  of  environmental  vorticity.  For  example,  a  TC  north  (south)  of  the 
subtropical  ridge  where  a  westward  gradient  of  environmental  vorticity  is  also  present  is 
predicted  to  have  a  more  northward  (westward)  and  faster  (slower)  propagation  in  gen¬ 
eral  agreement  with  composite  observations  in  the  western  North  Pacific  (c.g..  Fig.  1.1). 

4.  Numerical  model  initialization 

An  important  potential  application  of  the  BSAM  is  to  provide  the  initially 
symmetric  bogus  vortex  of  a  barotropic  numerical  model  with  the  correct  wavenumber 
gyre  structure  so  that  quasi-steady  propagation  occurs  immediately,  and  thus  overcome 
a  well-known  slow  bias.  Such  an  application  may  also  be  regarded  as  an  indirect  test 
of  the  validity  of  the  modeling  concepts,  techniques  and  assumptions  employed  in  de¬ 
veloping  the  BSAM.  For  a  quiescent  environment,  the  BSAM-predictcd  wavenumber 
one  gyre  produces  virtually  steady  propagation  regardless  of  whether  an  internally  or 
externally-derived  propagation  velocity  is  used.  Including  the  BSAM -generated  gyres  in 
the  numerical  model  initialization  reduces  the  along-track  forecast  error  from  more  than 
150  km  to  less  than  60  km.  This  preliminary  result  shows  that  the  BSAM  has  significant 
potential  as  an  initialization  tool  for  operational  barotropic  forecast  models  such  as 
SANBAR. 

The  unresolved  issue  regarding  nonsteady  /?-induced  propagation  in  anticyclonic 
shear  mentioned  above  precludes  a  precise  estimation  of  the  potential  4S  h  forecast  error 
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reduction  that  might  be  achieved  by  a  numerical  model  initialization  with  the  BSAM. 
Nevertheless,  the  along-track  (cross-track)  differences  of  up  to  280  km  (80  km)  between 
the  48  h  positions  of  TC  with  and  without  BSAM -initialization  illustrates  the  potential 
sensitivity  of  TC  track  prediction  to  the  initial  asymmetric  structure  of  the  TC  when  a 
spatially  variable  environment  is  present.  Using  the  adjustable  parameters  in  the  BSAM 
to  minimize  barotropic  forecast  model  track  error  in  some  statistical  sense  may  be  a  vi¬ 
able  approach  for  developing  an  initialization  scheme  that  can  improve  barotropic  nu¬ 
merical  model  forecast  skill,  and  provide  additional  insights  into  the  interaction  of  TC's 
with  the  surrounding  environment. 
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APPENDIX  A.  COMPOSITE  DATA  CONVERSION 

Composite  studies  have  generally  characterized  TC  propagation  in  relative  terms 
( e.g .,  speed  and  direction  differences)  using  a  rotated  coordinate  system  aligned  with 
storm  motion.  This  compositing  methodology  tends  to  make  theoretical  interpretation 
of  the  data  difficult  because  the  analytical  and  numerical  studies  predict  that  TC  prop¬ 
agation  will  possess  a  particular  orientation  with  respect  to  the  direction  of  the  large- 
scale  vorticitv  gradient.  In  particular,  a  rotated  storm-relative  coordinate  system  would 
tend  to  obscure  TC  propagation  associated  with  the  gradient  of  the  Coriolis  parameter, 
since  that  gradient  has  a  storm  and  environment-independent  northward  orientation. 
Thus,  part  of  the  difficulty  in  comparing  theory  with  the  composite  observations  may¬ 
be  readily  overcome  by  representing  TC  propagation  as  a  vector  quantity  in  a  north- 
oriented  earth-relative  coordinate  system. 

The  TC  motion  and  steering  flow  composite  data  are  taken  from  the  latitude,  di¬ 
rection,  speed  and  intensity  stratifications  of  Chan  and  Gray  (1982)  and  George  and 
Gray  (1976)  for  the  western  North  Pacific  region,  and  the  direction  and  recurvature 
stratifications  of  Holland  (1984)  for  the  Australian-Southwest  Pacific  region.  Holland 
used  a  single  steering  flow  definition  based  on  a  800  to  300  mb  pressure-weighted  mean 
wind  averaged  over  an  annulus  extending  5°  to  7°  lat.  from  the  TC  center.  Although 
Chan  and  Gray  used  the  same  horizontal  domain,  several  vertical  averaging  schemes 
were  tested.  Only  the  Chan  and  Gray  steering  flow  based  on  a  surface  to  300  mb  vertical 
average  is  used  here,  since  it  most  closely  approximates  the  steering  flow  definition  used 
by  Holland.  Vertically-averaged  steering  flows  have  been  chosen  for  this  analysis  rather 
than  individual  steering  levels  (c.g.,  George  and  Gray  1976)  to  more  appropriately  com¬ 
pare  the  observations  with  the  theoretical  modeling  results  that  are  predominantly  based 
on  barotropic  dynamics. 

Table  A.l  summarizes  the  conversion  process  for  the  western  North  Pacific  com¬ 
posite  data.  The  columns  labelled  VPM,  DD,  Vc  and  Dc  contain  previously  published 
data,  and  the  last  two  columns  have  been  computed  using  the  relationships 


Db  =  Dc  +  DD 


(-T1) 


VB 


^  C  +  ^  PM 

cos  DD 


(.1.2) 
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Table  A.l  Original  and  converted  composite  TC  motion  data  for  the  western  North 
Pacific  region.  Column  heading  meanings:  VPM  is  the  speed  of  the  steering  flow  compo¬ 
nent  parallel  to  the  direction  of  the  composite  TC  minus  the  speed  of  the  TC;  DD  is  the 
difference  between  the  direction  of  TC  motion  and  the  steering  flow;  Vc  and  Dc  are  the 
speed  and  direction  of  motion  of  the  TC  respectively;  and  VB  and  DB  are  the  speed  and 
direction  of  the  steering  flow  respectively.  The  data  in  columns  VPM  and  DD  are  taken 
directly  from  Chan  and  Gray  (1982),  and  the  data  in  columns  Vc  and  Dc  are  taken  di¬ 
rectly  from  George  and  Gray  (1976).  The  data  in  columns  VB  and  DB  have  been  com¬ 
puted  as  described  in  the  text.  Directions  are  measured  clockwise  from  North  and  the 
data  in  the  last  four  columns  are  relative  to  a  reference  frame  fixed  to  the  surface  of  the 
earth. 


Composite 

Stratification 

Vp\l 

(m  s) 

DD 

(deg) 

Vc 

(ms) 

Dc 

(deg) 

v. 

(ms) 

Db 

(deg) 

Latitude: 

>20°N  (GT) 

-1.0 

19 

5.6 

352 

4.9 

Oil 

<20°N(LT) 

-1.6 

6 

5.1 

300 

3.5 

306 

Direction: 

Westward  (W) 

-2.3 

17 

6.2 

285 

4.0 

302 

Northward  (N) 

-1.0 

17 

5.3 

324 

4.5 

341 

Eastward  (E) 

16 

7.1 

027 

6.8 

043 

Speed: 

Slow  (SE) 

27 

2.4 

338 

1.7 

005 

Moderate  (M) 

20 

5.2 

326 

4.4 

346 

East  (I  ) 

-1.3 

14 

10.1 

006 

9.1 

020 

Intensity: 

Weak  (WK) 

-1.1 

14 

4.9 

319 

3.9 

333 

Intense  (I) 

-1.1 

20 

5.0 

326 

4.2 

346 

Very  Intense  (VI ) 

-1.7 

26 

5.2 

319 

3.9 

345 

All  column  labels  arc  defined  in  the  table  caption.  Table  A. 2  is  analogous  to  Table  A.l 
for  the  Australian-Southwest  Pacific  composite  data.  The  columns  labelled  SD.  DD.  VB 


and  Db  contain  previously  published  data,  and  the  columns  labelled  Vc  and  Dc  have 
been  computed  using  (A.l)  and 


Vc  -  VB  -  SD  .  (A. 3) 

The  TC  propagation  vectors  shown  in  Fig.  1.1  were  then  computed  for  each  composite 
stratification  using  the  last  four  entries  in  each  row  of  Tables  A.l  and  A. 2. 


Table  A. 2  Analogous  to  Table  A.l  for  the  Australian-Southwest  Pacific  region.  The 
column  headings  DD,  Vc,  Dc,  VB  and  DB  have  the  same  meanings  as  in  Table  A.l  and 
SD  is  the  speed  difference  between  the  composite  TC  and  steering.  The  data  in  columns 
SD,  DD.  VB  and  DB  are  taken  directly  from  Holland  (1984),  except  that  the  steering  flow 
directions  are  measured  clockwise  from  North.  The  data  in  columns  Vc  and  Dc  have 
been  computed  as  described  in  the  text. 


Composite 

Stratification 

SD 
(m  s) 

DD 

(deg) 

Vc 

(m.s) 

Dc 

(deg) 

vB 

(ms) 

Db 

(deg) 

Direction: 

Westward  (W) 

1.3 

-4 

4.0 

247 

2.7 

243 

Soutlnvestward  (SW) 

1.1 

-18 

3.6 

241 

2.5 

223 

Southward  (S) 

i 

O 

-32 

3.0 

172 

3.5 

140 

Southeastward  (SE) 

-1.1 

-23 

3.8 

150 

4.9 

127 

Eastward  (E) 

-1.2 

-3 

4.2 

106 

5.4 

103 

Rccurvature: 

Before  (B) 

1.1 

-13 

3.5 

238 

2.4 

225 

Near  (N) 

0.5 

-20 

3.7 

193 

3.2 

167 

After  (A) 

-0.5 

-2 

5.0 

162 

5.5 

160 

APPENDIX  B.  PIECEWISE-ANALYTIC  VORTEX  CONSTRUCTION 


The  following  outline  describes  the  procedure  for  constructing  piecewise-analytic 
radial  profiles  of  TC  tangential  wind  and  relative  vorticity.  The  approach  here  is  to  use 
the  linear  degrees  of  freedom  in  (4.6),  (4.16)  and  (5.2)  to  develop  closed-form  expressions 
for  the  coefficients  A„  A2,  A3,  B„  B2,  B3  and  C3  such  that  the  piecewise-analytic  profiles 
must  have  continuous  windspeed  and  relative  vorticity.  The  remaining  nonlinear  degrees 
of  freedom  are  then  used  to  adjust  the  fit  of  the  piecewise-analytic  profiles  to  the  corre¬ 
sponding  analytic  profiles  by  selecting  various  values  for  X,  and  X2  and  assessing  the 
results  via  interactive  computer  graphics. 

1.  Given  the  trial  values  of  R0  and  R,  from  Step  1  of  the  model  solution  procedure 
(Chapter  IV  Section  A.4),  choose  trial  values  for  vs(R0)  and  vs(R,).  Generally  R0 
will  be  in  the  vicinity  of  the  radius  of  maximum  winds,  and  thus  v^RJ^l  will  pro¬ 
vide  a  good  first  estimate.  The  value  of  vs(R,)  for  the  piecewise-analytic  vortex  is 
chosen  to  closely  approximate  the  windspeed  of  the  analytic  profile  at  R,. 

2.  Determine  the  coefficients  A,  and  B,  for  (4.6)  in  the  first  annulus,  by  solving  the 
2x2  linear  system  that  results  from  applying  the  boundary  conditions  of  Step  1 
above  to  (4.6). 

3.  Determine  the  coefficients  A2  and  B2  by  requiring  the  windspeed  and  vorticity  of  the 
piecewise-analytic  TC  structure  to  be  continuous  at  R,.  As  in  Step  2  above,  this 
entails  solving  a  2x2  linear  system. 

4.  Repeat  the  process  in  Step  3  to  determine  the  coefficients  A3  and  B3  in  (4.16)  by 
requiring  the  piecewise-analytic  windspeed  and  vorticity  to  match  at  the  transition 
radius  Rr.  For  the  piecewise-analytic  profiles  used  in  Chapter  V.  this  step  is  mod¬ 
ified  to  include  evaluation  of  C3  in  (5.2)  by  solving  the  linear  3x3  system  that  results 
from  requiring  the  piecewise-analytic  windspeed,  vorticity  and  vorticity  gradient  to 
be  continuous  at  Rr. 

5.  Using  interactive  computer  graphics,  vary  the  parameters  X,  and  X2  to  adjust  the 
piecewise-analytic  windspeed  and  vorticity  profiles  to  approximate  the  analytic 
counterparts  as  described  in  Chapter  IV  Section  B.I.  During  this  step,  also  adjust 
Ro,  R,.  vs(R0)  and  vs(R,)  if  necessary  to  achieve  an  acceptable  fit. 
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